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Abstract 

The  goals  of  this  seedling  program  were  to  analyze  requirements  on  key  devices  and 
algorithms  to  realize  scalable  photonic  quantum  information  processing  (QIP)  systems 
that  can  significantly  outperform  classical  /  conventional  computing  methods. 
Specifically,  we  investigated  resources  requirements  for  linear  optical  quantum 
computing  to  achieve  scalable  quantum  computation  and  a  detailed  resource-cost  vs. 
performance  study  for  using  all-optical  quantum  computing  to  build  quantum  repeaters 
for  long-range  quantum- secured  communications  (viz.,  QKD). 

On  the  experimental  side,  we  investigated  advanced  photonic  integrated  circuit  (PIC) 
technology,  scaling  of  on-chip  superconducting  single-photon  detectors,  and 
requirements  for  feedforward.  This  work  showed  that  silicon  photonic  circuits  can  enable 
new  generations  of  quantum-limited  detectors  within  the  next  few  years,  and  that  an 
all-optical  scalable  architecture  of  quantum  computing  is  a  promising  approach,  but  also 
realized  that  photonic  quantum  computing  will  require  the  development  of  several  new 
theoretical  and  device  concepts.  Our  analysis  showed  that  PIC  technologies  that  could  be 
operated  at  or  near  room  temperature  would  bring  significant  advantages  for  integration 
with  high-speed  digital  electronics  and  room-temperature  single  photon  detector 
technologies,  and  we  propose  a  new  visible  photonics  platform  based  on  Ill-Nitride  that 
is  also  compatible  with  complementary  metal  oxide  semiconductor  (CMOS)  circuits. 

On  the  theory  side,  we  were  able  to  dramatically  reduce  resource  requirements  for  two 
key  QIP  applications:  photonic  cluster-state  quantum  computing  and  photonic  quantum 
repeaters.  Our  resource  analysis  on  ballistic  linear  optical  quantum  computing  (LOQC) 
considered  how  to  build  a  large  continuous  “sheet”  of  optical  cluster  state  that  is 
“connected  enough”  for  universal  quantum  computing.  The  sheet  of  cluster  state  is 
prepared  by  preparing  small  k-photon  cluster  states  (micro-clusters)  and  injecting  them 
into  a  linear  optical  circuit.  At  the  output,  one  detects  a  subset  of  the  modes,  and  the 
remainder  of  the  modes  carry  a  resource  universal  for  quantum  computing  regardless  of 
the  measurement  results  (the  measurement  results  tell  us  how  to  use  the  giant  entangled 
cluster  on  the  remaining  unmeasured  modes  to  map  any  quantum  algorithm  to  it).  The 


special  case  of  k=l  in  the  aforesaid  way  of  looking  at  Ballistic  quantum  computing  is 
exactly  the  Aaronson-Arkipov  model  of  boson  sampling.  Our  analysis  revealed  that  a 
clear  trade-off  between  the  size  of  the  input  clusters  N,  and  loss  tolerance  in  the  circuit, 
for  the  final  sheet  of  cluster  that  comes  out  to  be  universal  for  quantum  computing.  In 
particular,  by  using  new  techniques  from  percolation  theory,  we  showed  that  having  k=3 
size  clusters  as  a  starting  resource  is  a  must  for  loss-tolerant  Ballistic  LOQC  using 
two-qubit  fusion  operations.  This  work  establishes,  for  the  first  time,  a  clear  trade-space 
between  Aaronson-Arkipov  boson  sampling,  i.e.,  k=l  (which  does  not  allow  universal 
quantum  computing)  to  photonic  cluster  states  that  are  universal  for  quantum  computing 
(i.e.,  with  GHZ  states  of  size  k  >  3 ).  k=2  (Bell  states)  happens  to  be  a  border-line  case, 
where  Ballistic  LOQC  is  barely  possible  is  all  the  linear-optic  elements  are  100% 
lossless,  which  of  course  is  impractical.  This  study  tells  us  that,  for  scalable  LOQC,  we 
need  to  consider  scalable  and  efficient  ways  to  directly  produce  k>3  size  photonic 
clusters  (rather  than  first  producing  single  photons  and  probabilistically  stitching  them 
together). 

We  also  investigated  a  special  application  of  all-photonic  quantum  computing,  viz.,  in 
building  all-optical  quantum  repeaters  for  long-distance  QKD  and  entanglement 
distribution.  We  investigated  what  resources  are  required  to  beat  the  best-possible  rate 
achievable  without  repeaters  using  photonic  quantum  repeaters,  even  when  the  photonic 
repeater  is  constructed  using  lossy  devices  (lossy  sources  and  detectors,  in-line  losses, 
and  coupling  losses).  We  implemented  several  protocol  improvements  that  reduced 
physical  resources  for  such  schemes  by  5  orders  of  magnitudes  from  the  state  of  the  art  at 
the  beginning  of  the  seedling  program,  to  a  few  million  single  photon  sources  per 
repeater  node,  or  just  thousands  of  GHZ  states  per  repeater  station.  Towards  the 
conclusion  of  our  seedling  program,  we  further  improved  the  resource  requirement  by 
another  order  of  magnitude  by  developing  a  “one-way”  all-photonic  quantum  repeater 
scheme,  in  collaboration  with  Prof.  Liang  Jiang  of  Yale  U.,  which  will  be  published  soon. 

The  program  revealed  a  highly  promising  QIP  concept  based  on  weak  optical 
nonlinearities.  In  particular,  we  demonstrated,  for  the  first  time,  an  approach  for  a 
quantum  Ising  machine  based  on  weak  optical  nonlinearities.  As  opposed  to  currently 
pursued  approaches  based  on  superconducting  quantum  computing  architectures,  our 
all-optical  approach  could  shrink  devices  by  >  5  orders  of  magnitude,  and  could  operate 
at  room  temperature.  A  similar  architecture  also  holds  promise  to  single-photon 
nonlinearities  at  in  dielectric  cavities  that  promises  circuit-model  all-optical  quantum 
computing  analogous  to  superconducting  (SC)  qubits  based  quantum  computing,  but  at 
room-temperature  and  at  optical  frequencies. 

We  developed  scalable  methods  for  making  and  operating  arrays  of  both 
waveguide-integrated  and  fiber-coupled  superconducting  nanowire  single-photon 
detectors  (SNSPDs).  We  also  developed  high  performance  SNSPDs  on  AIN,  a 
wide-bandgap,  piezoelectric,  and  electro-optic  semiconductor.  SNSPD  on  AIN  can  work 
for  a  broad  spectrum  range  from  UV  to  mid-IR,  and  can  be  potentially  combined  with 
other  on-chip  superconducting  electronics  to  achieve  optical  modulation.  In  conjunction 


with  detector  development,  we  also  built  cryogenic  components  necessary  for  large  array 
operation,  including  cryogenic  amplifiers  and  flexible  RF  cable  bundles. 

1.  Progress  on  Experimental  Devices  and 
Platforms 

1. 1  Demonstration  of  high-fidelity  linear-optics  mode 
transformations 

Realizing  scalable,  high-fidelity  interferometric  networks  is  a  central  challenge  to  be 
addressed  on  the  path  towards  linear  optical  quantum  computation  as  well  as  for 
mediating  optical  interactions  between  nonlinear,  matter-based  qubits.  We  demonstrated 
a  programmable  nanophotonic  processor  composed  of  88  ultra-high  contrast 
Mach-Zehnder  interferometers  each  exhibiting  a  record  extinction  ratio  exceeding  80  dB, 
see  below.  We  benchmarked  the  performance  of  the  processor  on  single-qubit  rotations 
and  showed  the  first  experimental  demonstration  of  a  new,  error-tolerant  universal 
unitary  encoding  protocol  by  implementing  9-dimensional  unitary  transformations 
sampled  from  the  Haar  measure.  We  also  showed  how  the  fidelity  of  these 
transformations  can  be  boosted  using  in-situ  nonlinear  optimization  techniques.  In 
addition,  we  introduce  new  methods  for  characterizing  these  large  interferometric 
networks. 
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Figure  1  Photonic  integrated  circuit.  Left:  programmable  PIC.  Right:  Transmission  at  one 
of  the  output  modes,  as  one  of  the  internal  MZI’s  phase  shifters  is  modulated. 

1.2  On-chip  nanowire  detector  array 


We  developed  high-performance  NbN  nanowire  detectors  on  AIN  PICs.  These  detectors  have 
near-unity  quantum  efficiency,  sub-6-ns  reset  time,  and  ~60  ps  timing  jitter.  Since  they  were 
fabricated  directly  on  the  waveguides,  the  optical  coupling  was  efficient  and  low-loss.  Figure 


2  shows  a  scanning  electron  micrograph  of  AIN  waveguide-integrated  SNSPDs.  The 
detectors  are  arranged  in  such  a  way  that  on-chip  photon  correlation  measurement  is  possible. 


Figure  2:  AIN  waveguide-integrated  SNSPDs. 


To  read  out  a  large  number  of  detectors,  we  proposed  a  time-tagged  readout  scheme.  Figure  3 
shows  a  prototyping  time-tagged  SNSPD  array  fabricated  on  AIN.  This  array  consists  of  4 
chains,  each  connecting  16  detectors,  which  makes  a  total  number  of  64  detectors.  The  time 
tagging  is  realized  by  separating  detectors  using  superconducting  slow-wave  transmission 
lines.  Impedance  matching  taper  is  needed  to  transform  the  high-impednace  SNSPD  to  50  Q 
readout  electronics.  This  readout  scheme  is  scalable,  and  greatly  reduces  the  resources 
needed  for  large-array  cryogenic  readout. 


Figure  3:  Time-tagged  SNSPD  array  fabricated  on  AIN.  The  array  consists  of  4  detector 
chains,  each  connecting  16  detectors  in  series.  Detectors  are  separated  using  superconducting 
slow-wave  transmission  line,  and  terminated  with  impedance-matching  tapers  to  50  Cl 
readout  circuits. 

1.3  Fiber-coupled  detector  array 

Fiber-coupled  detector  array  can  be  used  for  off-chip  single  photon  detection.  We 
developed  a  new  packaging  technique  to  couple  light  from  multiple  optical  fibers  to  a 
detector  array  on  a  single  chip.  Figure  4  shows  an  8-channel  fiber-coupled  SNSPD  array. 
The  detectors  were  fabricated  on  a  silicon  nitride  substrate,  and  arranged  to  match  the 
spacing  of  a  commercial  fiber  array.  The  fiber  array  was  aligned  and  glued  to  the  detector 
chip.  These  detectors  have  a  circular  active  area  with  1 5  pm  in  diameter,  and  the  device 
yield  is  >  80%.  This  packaging  technique  can  be  extended  to  larger  array  size  without 
much  additional  effort. 

To  enable  large  array  operation,  we  also  developed  compact  cryogenic  amplifier  and 
flexible  RF  cable  bundles.  These  cryogenic  components  are  essential  for  both 
waveguide-integrated  and  fiber-integrated  SNSPD  arrays. 


Figure  4:  An  8-channel  fiber-coupled  SNSPD  array. 


1.4  Post-fabrication-tunable  linear  optic  fabrication 

We  have  analyzed  the  performance  of  the  programmable  nanophotonic  processor  (PNP)  that 
is  dynamically  tunable  via  post-fabrication  active  phase  tuning  to  predict  the  scaling  of  such 
PICs  from  tens  to  hundreds  of  optical  channels.  In  particular,  we  have  demonstrated  that  the 
PNP  can  perform  the  mode  transformations  necessary  for  a  8x8  Green  Machine[l], 


Figure  GM1.  Green  Machine  circuit  implemented 
using  Mach-Zehnder  interferometers  in  the  PNP  array 
configuration. 


Building  on  the  development  of  a  circuit  (Figure  GM1)  to  implement  the  Green  Machine 
on  the  programmable  nanophotonic  processor  (PNP),  control  hardware  and  software  has 
been  engineered  to  calibrate  the  PNP  and  manipulate  light  routing  within  the  chip.  These 
features  are  necessary  for  configuring  the  precise  phase  settings  required  to  demonstrate 
the  Green  Machine  circuit  on  a  physical  chip  with  non-idealities.  Setup  and  light 
propagation  for  the  currently  controlled  PNP,  which  has  5  input  modes  and  1 1  output 
modes,  is  shown  in  Figure  GM2. 


Figure  GM2.  (a)  Hardware  control  of  thermal  heater  voltages  and  optical  coupling  setup  for  the  current 
PNP  chip.  Not  pictured:  photodiode  array  setup  used  to  detect  light  coupled  from  the  output  waveguides, 
(b)  Propagation  of  light  through  the  chip,  controlled  by  the  setup  in  (a). 


QPPO: 

The  PNP  circuit  shown  in  Figure  1  requires  12  input  modes,  which  is  larger  than  the 
existing  fabricated  chip,  hereby  referred  to  as  QPPO.  An  alternative  test  circuit  (Figure 
GM3)  was  designed  to  simulate  the  features  of  the  Walsh  transform  on  QPPO,  accounting 
for  the  decreased  number  of  input  modes.  Given  that  every  other  output  mode  on  the  chip 
is  not  routed,  the  power  on  only  four  output  modes  can  be  detected  simultaneously; 
therefore,  the  circuit  needed  to  be  altered  to  read  all  eight  modes  from  the  output.  The 
output  probability  distribution  was  measured  for  this  test  circuit  with  very  high  fidelity, 
as  shown  in  Figure  GM4. 


Figure  GM3.  Test  circuit  for  demonstrating  features  of  the  Walsh  transform  on  the  compressed  QPPO  chip. 
Mach  Zehnder  interferometer  (MZI)  internal  arm  phase  differences  are  denoted  in  the  gray  boxes.  Here,  nil 
performs  a  2x2  Hadamard  operation,  n  performs  an  identity,  and  0  performs  a  perfect  swap. 


_ 


!5 

05 

_Q 

2 

0. 


OJ 

"O 

o 


0.20 

0.15 

0,10 

0,05 

0,00 


□ 


I  I  Experiment 

F  =  99.9949  dt  0.0033%  “  simulation 


Mode  1  Mode  2  Mode  3  Mode  4  Mode  5  Mode  6  Mode  7  Mode  8 


Figure  GM4.  Output  probability  distribution  of  Green  Machine  test  circuit  with  one  excited  input  mode. 
Fidelity  of  the  experimental  data  with  respect  to  software  simulations  of  the  circuit  yield  high  fidelity  of 
operation  of  99.9949±0.0033%. 

PNPlb: 

A  larger  chip,  PNPlb,  consisting  of  26  input  modes,  26  output  modes,  and  176 
individually-controlled  thermal  phase  shifters  (Figure  5)  has  been  fabricated  and 
packaged.  The  software  for  the  QPPO  has  been  designed  to  be  compatible  with  PNPlb  as 
well.  A  free  space  coupling  system  has  been  developed  and  is  being  tested  that,  in 
conjunction  with  a  PID  thermal  control  loop,  could  enable  coupling  to  and  from  the  chip 
that  is  approximately  ldB  less  than  that  which  can  be  obtained  with  standard  edge 
coupling.  The  free  space  system  makes  use  of  a  beam  expansion  system  deployed  using 
two  AR-coated  lenses,  shown  in  Figure  GM6. 
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Figure  GM5.  Layout  of  the  larger  26-mode,  264-heater  PNP  that  is  currently  packaged.  MZI  internal  phase 
settings  corresponding  to  an  implementation  of  the  Green  Machine  circuit  are  shown  in  the  shaded  boxes. 


Figure  GM6.  (a)  Beam  expansion  system  schematic  for  5x  demagnification  of  fiber  array  modes  to  one 
facet  of  the  PNP  chip.  Using  this  lens  system  allows  for  coupling  from  a  fiber  array  with  a  pitch  and  mode 
field  diameter  (MFD)  of  127pm  and  10pm,  respectively,  to  the  chip,  having  waveguide  pitch  and  MFD  of 
25.4pm  and  2pm,  respectively,  (b)  Image  of  the  setup,  using  an  objective-lens  system  to  couple  to  the  left 
and  a  lens  beam  expansion  system  on  the  right. 


Detection: 

We  are  designing  a  homodyne  detector  array  to  realize  shot-noise  limited  detection  of  the  output  modes. 
InGaAs  P-I-N  photodiodes  obtained  from  Beijing  SWT  Optical  Communication  Technology  constitute  the 
detectors  in  this  scheme. 


Task  2:  Resource-performance  scaling  in  optical 
quantum  processing 

2.1  Resource  scaling  for  an  LOQC  implementation  in  PIC-based 
linear  optics 

The  goal  of  this  task  was  to  estimate  resource  scaling  requirements  for  a  useful  all-optical 
quantum  computing  application.  It  is  known  that  general-purpose  quantum  computers  are 
in  principle  physically  possible  [2].  Recently,  resource  estimates  for  linear-optics 
quantum  computing  systems  have  been  produced  [3],  though  using  non-optimal 
encoding.  In  fact,  the  optimal  encoding  depends  critically  on  the  application,  which 
makes  the  problem  statement  poorly  defined.  To  address  a  well-defined,  specific 
problem,  we  proposed  to  instead  analyze  the  recently  proposed  concept  of  all-optical 
quantum  repeaters  [4]  that  uses  photon  sources,  linear  optics,  photon  detectors  and 
classical  feedforward  at  each  repeater  node,  but  no  quantum  memories.  We  considered 
resource  costs  assuming  realistic  (though  aggressive)  device  performance  for  the  photon 
sources,  detectors,  and  PIC-based  mode  transformations.  This  analysis  quickly  revealed 
that  the  original  proposal  would  requirement  astronomically  many  single  photon  sources 
and  detectors  per  repeater  node. 


In  theoretical  work  about  to  be  published  in  Physical  Review  A,  we  showed  that  the 
quantum-secure  key  rate  has  the  form  R(r\)  =  Dr\s  bits  per  mode,  where  r|  is  the 
end-to-end  channel's  transmissivity,  and  the  constants  D  and  s  are  functions  of  various 
device  inefficiencies  and  the  resource  constraints,  such  as  the  number  of  available  photon 
sources  at  each  repeater  node.  Even  with  lossy  devices,  we  show  that  it  is  possible  to 
attain  s<l,  and  in  turn  outperform  the  maximum  key  rate  attainable  without  quantum 
repeaters,  R_direct(r|)=— log2(l-r|)~(l/ln2)r|  bits  per  mode  for  q<  1,  beyond  a  certain 
total  range  L,  where  q~exp(-aL)  in  optical  fiber.  We  also  propose  a  suite  of 
modifications  to  original  all-optical  repeater  protocol[4],  which  lower  the  number  of 
photon  sources  required  to  create  photonic  clusters  at  the  repeaters  so  as  to  outperform 
R_direct(ri),  from  -TO11  to  ~106  photon  sources  per  repeater  node.  We  show  that  the 
optimum  separation  between  repeater  nodes  is  independent  of  the  total  range  L,  and  is 
around  1.5  km  for  assumptions  we  make  on  various  device  losses. 

This  work  is  available  on  ArXiv[5]  and  accepted  for  publication  in  PRA.  The  detailed 
paper  is  also  attached  at  the  end  of  this  report  as  an  addendum. 

2.2  Ballistic  quantum  computing  with  microclusters 

It  has  been  argued  that  BosonSampling  —  sampling  from  the  joint  photon-number 
distribution  at  the  output  of  a  large  multiport  interferometer  when  multiple 
identical  bosons  are  input  into  it  -  is  computationally  hard  for  classical  computers, 
but  by-definition  efficient  on  a  quantum  system  [6],  since  the  aforesaid  (quantum) 
system  samples  efficiently  from  that  distribution  just  by  the  virtue  of  being  itself.  A 
few  practical  uses  of  BosonSampling  have  been  proposed  in  quantum 
simulation[7-9],  but  the  scheme  does  not  permit  general  purpose  quantum 
computation  —  or,  so  it  is  believed.  On  the  other  hand,  it  was  recently  shown  that 
cluster  states  universal  for  quantum  computation  could  be  constructed  by 
probabilistically  fusing  3-photon  GHZ  states,  i.e.,  inputting  3-photon  clusters  at  the 
input  of  a  large  interferometer,  and  without  the  use  of  feed-forward  and  switching, 
the  output  joint  quantum  state  is  universal  for  quantum  computing  [10],  in  the 
sense  that  no  matter  what  is  the  actual  shape  of  the  instance  of  the  sheet  of 
entangled  cluster  that  comes  out  of  the  linear  optical  interferometer,  one  can  always 
detect  a  portion  of  it  (some  subset  of  modes),  and  use  the  measurement  results  to 
determine  how  to  map  any  quantum  algorithm  instance  into  the  cluster  that  lives  on 
the  remainder  of  the  modes.  This  result  by  Terry  Rudolph  and  collaborators, 
published  in  2015,  marked  a  major  advance  that  suggested  that  building  a 
linear-optical  quantum  computer  may  be  far  less  challenging  than  previous 
estimates  [3]. 

During  the  course  of  our  SEQUOIA  program,  members  of  our  team  were  able  to 
make  significant  additional  improvements  on  the  efficiency  of  Ballistic  cluster  state 
generation.  Most  notably, 

(1)  we  were  able  to  establish,  for  the  first  time,  a  fundamental  trade-off  between 
the  inline  loss  tolerance  in  the  source-detection-efficiency  product  of  each 


source  and  every  detector  in  the  system,  and  the  size  k  of  the  initial 
micro-cluster  resource,  so  that  universal  quantum  computing  is  possible,  as 
shown  in  the  figure  below. 

(2]  Rudolph’s  work  showed  that  as  long  as  linear-optic  Fusion  operation  at 
success  probability  >  62.5%  is  possible,  one  can  do  universal  QC  with  k=3 
(GHZ)  micro-clusters  as  the  starting  resource.  We  reduced  that  achievability 
threshold  to  59%.  We  also  showed  sufficient  evidence  that  54%  is  achievable 
(full  proof  is  yet  to  be  done),  but  showed  that  one  has  to  have  a  success 
probability  of  Fusion  greater  than  l/(k-l),  which  for  k  =  3  microclusters  is 
50%.  So,  we  have  essentially  closed  the  efficiency  gap  for  Ballistic  QC. 

(3)  In  Rudolph’s  work  that  demonstrated  the  62.5%  achievable  threshold  for 
universal  Ballistic  QC,  they  used  3D  clusters.  We  show  it  is  possible  with  a  2D 
cluster,  which  reduces  the  device  requirements  tremendously.  The  reason 
this  threshold  was  an  important  result  is  that  Peter  van  Loock  and  others  had 
shown  in  recent  work  that  a  Bell  measurement  (Fusion)  at  75%  success  rate 
is  possible  using  a  linear-optics  based  scheme  with  single  photons  injected 
into  the  passive  linear-optic  circuit,  if  all  the  devices  are  lossless.  Without 
those  injected  photons,  the  threshold  can  be  no  larger  than  50%.  Since  62.5% 
<  75%,  Rudolph  et  al.  had  argued  that  Ballistic  QC  is  possible  as  long  as 
3-photon  GHZ  states  are  available  to  start  with.  We  have  shown  a  particular 
2D  lattice,  that  a  74%  success  probability  suffices  the  percolate  the  output 
cluster,  and  hence  make  universal  Ballistic  QC  possible  with  this  approach, 
and  hence  doable  in  principle  with  the  aforesaid  75%  boosted  Fusion  gate. 

Remarkably,  our  work  shows  that  experimentally  near-term  efficiencies  may  be 
tolerable  for  even  moderately  sized  input-cluster  states.  Experimental  progress 
towards  systems  that  can  produce  such  cluster  states  has  been  surprisingly  fast: 
recently,  it  was  demonstrated  experimentally  that  photonic  cluster  states  can  be 
produced  deterministically  by  photon  scattering  off  a  lambda-level  system  in 
semiconductor  quantum  emitters  [11],  thought  the  efficiency  is  still  low  and  the 
system  requires  cryogenic  cooling.  As  an  alternative  that  may  be  operated  at  room 
temperature  and  fully  lithographically  controlled  (as  opposed  to  quantum  dots),  we 
have  developed  dielectric  nonlinear  sources  of  single  photons  and  small  entangled 
states[12];  see  Section  3.3. 
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Figure  2.2.1  Illustration  of  generation  of  scalable  cluster  state  from  micro-cluster  input,  k  =  1  (Boson 
sampling),  output  cluster  is  disconnected:  not  a  good  enough  resource  for  universal  QC;  k  =  2  (Bell  pair 
initial  resource),  output  barely  percolates  if  all  efficiencies  are  100%  (no  loss  tolerance);  k  =  3  (3-photon 
GFIZ  is  initial  resource),  then  ballistic  QC  is  possible;  3D-diamond  (Rudolph  et  al.);  Xc  =  62.5%;  loss 
tolerance:  r)sr)d  >  96.2%;  (10,3)-b  lattice  (our  method);  he  =  54%;  loss  tolerance:  r)sr)d  >  93.4% 

Best  possible  lattice  (converse  proof  /  Bethe  lattice);  he  =  50%;  max  possible  loss  tolerance:  91.6% 

Sources  of  >=  3-qubit  microclusters  needed  if  we  want  to  avoid  feed- forward,  switching,  and  associated 
losses 

Our  resource  analysis  on  ballistic  linear  optical  quantum  computing  (LOQC)  considered 
how  to  build  a  large  continuous  “sheet”  of  optical  cluster  state  whose  connectivity  is 
pervasive  and  long-ranging  enough  to  renormalize  it  (i.e.,  continuously  carve  out  a 
logical  cluster  state  from  the  random  large  spanning  cluster  that  comes  out  of  the 
interferometer)  for  universal  quantum  computing.  The  sheet  of  cluster  state  is  prepared 
by  preparing  small  k-photon  cluster  states  (micro-clusters)  and  injecting  them  into  a 
linear  optical  circuit.  At  the  output,  one  detects  a  subset  of  the  modes,  and  the  remainder 
of  the  modes  carry  a  resource  universal  for  quantum  computing  regardless  of  the 
measurement  results  (the  measurement  results  tell  us  how  to  use  the  giant  entangled 
cluster  on  the  remaining  unmeasured  modes  to  map  any  quantum  algorithm  to  it).  The 
special  case  of  k=l  in  the  aforesaid  way  of  looking  at  Ballistic  quantum  computing  is 
exactly  the  Aaronson-Arkipov  model  of  boson  sampling.  Our  analysis  revealed  that  a 
clear  trade-off  between  the  size  of  the  input  clusters  k,  and  loss  tolerance  in  the  circuit, 
for  the  final  sheet  of  cluster  that  comes  out  to  be  universal  for  quantum  computing.  In 
particular,  by  using  new  techniques  from  percolation  theory,  we  showed  that  having  N=3 
size  clusters  as  a  starting  resource  is  a  must  for  loss-tolerant  Ballistic  LOQC  using 
two-qubit  fusion  operations.  This  work  establishes,  for  the  first  time,  a  clear  trade-space 
between  Aaronson-Arkipov  boson  sampling,  i.e.,  k=l  (which  does  not  allow  universal 
quantum  computing)  to  photonic  cluster  states  that  are  universal  for  quantum  computing 
(i.e.,  with  GHZ  states  of  size  k  >  3 ).  k=2  (Bell  states)  happens  to  be  a  border-line  case, 
where  Ballistic  LOQC  is  barely  possible  is  all  the  linear-optic  elements  are  100% 
lossless,  which  of  course  is  impractical. 

This  study  told  us  that,  for  scalable  LOQC,  we  need  to  consider  scalable  and  efficient 
ways  to  directly  produce  k>3  size  photonic  clusters  (rather  than  first  producing  single 


photons  and  probabilistically  stitching  them  together).  This  observation  spurred  us  to 
investigate  a  structured  means  to  directly  generate  3 -photon  entangled  GHZ  state 
micro-clusters  by  exploiting  six-wave  mixing  in  a  non-linear-optical  waveguide.  This 
idea,  and  associated  calculations  will  be  described  in  Task  3  below. 

Task  3:  Exploratory  survey  of  key  enabling 
technologies  for  scalable  optical  information 
processing 

We  have  performed  a  thorough  exploratory  study  surveying  various  key  enabling 
technology  components — including  ones  not  being  developed  or  pursued  in  Task  1  of  this 
project  for  optical  quantum  information  processing. 


3.1  Wide-bandgap  photonic  integrated  circuits  for  quantum 
information  processing  in  the  visible  spectrum 

How  discussion  so  far  on  photonic  integrated  circuits  has  focused  on  the  silicon  on 
insulator  (SOI)  platform,  which  is  technologically  the  most  mature.  However,  this 
platform  also  has  several  important  shortcomings:  high-speed  modulators  are  based 
on  free  carriage  dispersion  effects,  which  by  their  nature  introduce  some 
non-negligible  loss  and  do  not  work  well  at  cryogenic  temperatures  because  of 
carrier  freezeout;  in  addition,  the  circuit  requires  infrared  photons  above  about  1.1 
q  in  wavelength  come  out  which  in  turn  requires  superconducting  detectors  that 
add  complications  as  compared  to  silicon  avalanche  photodetectors  that  could  be 
used  at  shorter  wavelengths.  Therefore,  we  investigated  here  the  possibility  of  an 
entirely  different  type  of  photonic  integrated  circuit  based  on  the  white  bandgap 
Ill-nitrite  material  platform.  This  platform  allows  for  electro-optic  modulators 
based  on  chi-2  coefficient,  which  function  well  at  low  temperature  and  do  not  add 
free-carrier  dispersion  loss.  Furthermore,  wavelengths  in  the  visible  and  potentially 
even  the  ultraviolet  are  compatible,  allowing  the  use  of  room-temperature  APDs. 
Without  the  need  for  superconducting  detectors,  the  entire  photonics  platform  can 
operate  at  room  temperature.  Moreover,  the  Ill-nitride  platform  allows  for  Kerr 
nonlinearities  as  well  as  chip-integrated  light  sources  and  detectors. 

Specifically,  we  analyzed  a  PIC  platform  comprised  of  a  crystalline  AlxGal-xN  optical 
guiding  layer  on  an  AIN  substrate  for  the  ultraviolet  to  visible  (UV-vis)  wavelength 
range.  An  A1  composition  of  x~0.65  provides  a  refractive  index  difference  of  ~  0.1 
and  a  small  lattice  mismatch  (<  1%)  that  minimizes  crystal  dislocations  at  the 
AlGaN /AIN  interface.  This  small  refractive  index  difference  is  beneficial  at  shorter 
wavelengths  to  avoid  extra-small  waveguide  dimensions.  The  platform  enables 
compact  waveguides  and  bends  with  high  field  confinement  in  the  wavelength  range 


from  700  nm  down  to  300  nm  [and  potentially  lower]  with  waveguide  cross-section 
dimensions  comparable  to  those  used  for  telecom  PICs  such  as  silicon  and  silicon 
nitride  waveguides,  allowing  for  well-established  optical  lithography.  This  platform 
can  potentially  enable  cost-effective,  manufacturable,  monolithic  UV-vis  photonic 
integrated  circuits.  This  study,  which  was  a  collaboration  between  BBN  (Mohammad 
Soltani),  MIT  (Dirk  Englund  and  Tomas  Palacios]  and  the  University  of 
Massachusetts  at  Boston  (Richard  Soref],  will  appear  in  Optics  Express. 


3.2  Review  of  atom-like  solid-state  on-demand  single  photon 
sources 

Single  photon  sources  represent  an  essential  resource  in  the  applications 
considered  in  this  program.  However,  experimentally,  we  are  still  some  ways  from 
an  'ideal'  on-demand  single-  photon  emitter.  A  wide  range  of  promising  material 
systems  have  been  developed,  and  several  have  transitioned  from  proof-of-  concept 
to  engineering  efforts  with  steadily  improving  performance.  We  reviewed  recent 
progress  in  the  race  towards  the  ideal  single-photon  emitter  required  for  a  range  of 
quantum  information  processing  applications.  We  focused  on  solid-state  systems 
including  quantum  dots,  defects  in  solids,  two-dimensional  hosts  and  carbon 
nanotubes,  as  these  are  well  positioned  to  benefit  from  recent  breakthroughs  in 
nanofabrication  and  materials  growth  techniques. 

The  central  performance  metrics  of  single  photon  sources  are:  source  purity  (lack  of 
multi-photon  emissions];  photon  indistinguishabilitiy  (every  photon  emitted  has  the 
same  mode  shape];  and  efficiency  (product  of  internal  quantum  efficiency  and 
collection  efficiency  into  a  desired  electromagnetic  mode].  Our  review  of  recent 
works  showed  major  deficiencies  in  the  reporting  of  these  source  qualities. 
Reporting  was  most  rigorous  for  InAs/GaAs  self-assembled  quantum  dots,  which 
were  also  the  most  advanced  technologically,  indistinguishability  ~  99%,  efficiency 
~  75%,  and  photon  purity  g(2](0]<l%. 

A  review  article  on  these  sources  is  published  in  Nature  Photonics[13] 


3.3  Proposal  for  on-demand  single  photon  sources  based  on 
quantum  feedback  of  a  nonlinear  dielectric  cavity 

As  mentioned  above,  increased  infidelity  in  the  single  photon  states  produced  by 
sources  sharply  increases  the  resource  overhead  for  quantum  repeaters  and 
ballistic  quantum  computing.  Single  photon  sources  based  on  atomic  emitters  have 
improved  greatly  over  recent  years  —  for  example,  emission  from  InAs  quantum 
dots  can  now  achieve  indistinguishability  between  consecutive  photons  in  excess  of 
99%  at  an  efficiency  greater  than  75%  [13]  —  but  their  performance  still  does  not 
reach  the  desired  >99%  or  so [5].  Moreover,  single  photon  sources  based  on 


quantum  emitters  in  solids  typically  have  considerable  inhomogenous  distribution 
in  emission.  On  the  other  hand,  single  photon  sources  based  on  heralded  entangled 
pair  generation  can  be  very  bright,  operate  at  room  temperature,  have  excellent 
photon  indistinguishability.  Moreover,  many  heralded  sources  based  on 
nanophotonic  circuits  can  be  tuned  perfectly  into  resonance[14,15].  A  heralded 
single  photon  source  relies  on  the  information  obtained  from  measurements  of  the 
device  to  predict  its  quantum  state.  It  is  normally  implemented  using  a  photon  pair 
(signal  and  idler)  source  and  a  single  photon  detector  that  performs  measurements 
on,  e.g.  the  idler  Hilbert  space.  The  success  probability  of  single  photon  generation 
can  be  increased  by  creating  photons  in  multiple  modes  and  multiplexing  them  into 
a  single  output  mode.  All  the  degrees  of  freedom  of  photons  may  be  used  as  modes 
and  multiplexing 

both  spatial,  temporal,  and  frequency  modes  have  been  demonstrated.  Spatial-  and 
frequency  multiplexing  can  produce  modes  in  parallel  from  several  individual 
sources  or  a  frequency  comb,  respectively.  Temporal  multiplexing  is  a  serial  process 
where  a  single  (a)  source  is  triggered  multiple  times. 

In  our  recent  work,  we  considered  temporal  multiplexing,  with  the  option  to  also 
use  spectral  multiplexing.  The  proposed  architecture  is  shown  in  Fig.  la,  consisting 
of  an  ultrahigh  Q  microcavity  with  a  X  (3)  nonlinearity  that  allows  signal  and  idler 
photons  to  be  generated  in  pairs.  The  resonator  is  coupled  to  a  switchable  driving 
laser,  a  frequency-  and  photon  number  resolving  detector  (FPNRD),  and  an  output 
channel  via  frequency  selective  tunable  gates  that  control  the  coupling  rates  to  the 
resonator. 


Fig.  3.3.1:  (a)  Proposed  on-demand  single  photon  source  based  on 
dynamic  cavity  storage.  Subscripts  i,  p,  s  refer  to  idler,  pump,  and 
signal  wavelength.  Solid  lines  are  optical  waveguides,  while  dashed 


lines  represent  electrical  control  signals.  FP-  NRD:  Frequency-  and 
photon  number  resolving  detector,  (b)  Example  of  a  gate 
implementation  that  has  a  large  coupling  for  the  pump  frequency  but 
is  uncoupled  for  the  signal  and  idler  frequencies.  This  is  achieved  by 
matching  the  arm  length  difference  of  the  MZI  to  the  free  spectral 
range  of  the  micro  cavity. 

Figure  3.3.2  plots  the  single-photon  fidelity  Fth  (%)  against  the  probability  of  a 
heralding  event  in  a  given  clock  cycle,  for  a  range  of  detector  efficiencies  r| .  This 

device  uses  temporal  heralding  in  a  single  ring  only.  With  a  heralding  probability  of 
even  just  P(success)=30%  [corresponding  to  r\  =0.99  and  Fth=0.98) ,  spectral 

multiplexing  over  15  multiple  sources  increases  the  probability  of  heralding  at  least 
one  event  increases  to  1-(1-0.3)15>99.5%.  Thus,  such  a  source  could  produce  single 
photons  of  sufficient  efficiency  and  fidelity. 


80 


95  96  97  98  99 


Fth  [*] 


Figure  3.3.2:  Fidelity  vs  heralding  trade-off. 

3.4  On-demand  GHZ  source  based  on  QND  and  dynamic  cavity 
control  (manuscript  in  preparation) 

As  mentioned  above,  3-photon  GHZ  states  form  a  minimal  resource  for  loss-tolerant 
production  of  cluster  states  for  ballistic  quantum  computing.  However,  no  efficient 
3-photon  GHZ  state  source  exists  today.  To  this  end,  we  have  extended  the  concept 
of  the  on-demand  single  photon  source  by  cavity  feedback  (Section  3.3)  to  produce 
3-photon  GHZ  states.  The  central  idea  is  to  monitor  the  cavity  photon  population 
using  quantum  nondemolition  measurements  of  the  cavity  modes,  using  cross-Kerr 


nonlinear  interactions  with  the  probe  beam  and  homodyne  detection,  as  shown  in 
the  figure  below. 
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3.5  3-photon  GHZ  source  based  on  6-wave  mixing  (manuscript  in 
preparation) 


Another  way  to  produce  3-photon  GHZ  states  is  by  6-wave  mixing,  as  illustrated 
below.  We  are  currently  investigating  this  method. 
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3.6  Quantum  Logic  with  Interacting  Bosons  in  ID 

We  also  developed  a  scheme [16]  for  implementing  high-fidelity  quantum  logic  gates 
using  the  quantum  walk  of  a  few  interacting  bosons  on  a  one-dimensional  lattice. 
The  gate  operation  is  carried  out  by  a  single  compact  lattice  described  by  a 
one-dimensional  Bose-Hubbard  model  with  only  nearest-neighbor  hopping  and 
on-site  interactions.  We  find  high-fidelity  deterministic  logic  operations  for  a  gate 
set  [including  the  CNOT  gate]  that  is  universal  for  quantum  information  processing. 
We  analyzed  the  applicability  of  this  scheme  in  light  of  recent  developments  in 
controlling  and  monitoring  cold-atoms  in  optical  lattices,  as  well  as  an 
implementation  with  realistic  nonlinear  quantum  photonic  devices.  This  work  is 
currently  under  review.  Figure  3.6.1  shows  the  evolution  of  a  two  dual-rail  photon 
states  (control  and  signal]  undergoing  a  controlled  logic  gate,  assuming  a 
two-photon  interaction  (such  as  mediated  by  the  Kerr  nonlinearity  or  an  atomic 
nonlinearity].  For  details,  see  Ref.  [16] 
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Figure  3.6.1  An  implementation  of  the  controlled  NOT  (CNOT)  gate  (a)  The  real  part 
and  (b)  the  imaginary  part  of  the  two-particle  unitary  transform,  U.  The  CNOT  gate 
operation  corresponds  to  the  sub-matrix  of  the  logic  states,  shown  in  solid-color 
bars  and  marked  with  red  axis  labels.  Plots  (c)-(f)  show  the  position  (in  terms  of  the 
lattice  sites,  1-4)  of  the  particle  density  as  a  function  of  time,  t,  revealing  the 
operation  principle  of  the  gate  on  each  logical  state  (|00>,  |10>,  |01>,  and  |11> 
respectively).  One  observes  that  the  target  qubit  (in  sites  3  &  4)  performs 
Rabi-oscillations  that  are  perturbed  by  the  state  of  the  control  qubit  (in  sites  1  &  2) 
—  the  target  qubit  performs  one  fewer  Rabi-flip  if  the  control  qubit  is  in  the  |1> 
state. 


3.7  Limitations  of  two-level  emitters  as  non-linearities  in  passive 
two-photon  controlled  phase  gates 

Two-level  atomic  systems  in  cavities  are  often  cited  as  one  possible  mechanism  for 
producing  two-photon  phase  gates  (or  other  two-photon  logic  gates).  We 
investigated  various  architectures  and  found  that  it  is  not  possible  to  reach  unity 
fidelity  for  such  gates  for  a  time-invariant  cavity. 

Using  a  "dual  Hong-Ou-Mandel”  geometry  shown  in  Figure  3.7.1,  we  were  able  to 
ensure  that  the  incoming  and  existing  photon  states  are  in  the  dual-rail  logic 
representation. 
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Figure  3.7.1  Schematic  of  two-photon  phase  gate. 

We  analyzed  such  gates  under  the  action  of  finite-duration  photon  wave  packets  for 
signal  and  control  photons.  The  two-photon  gate  fidelity  is  plotted  in  Figure  3.7.2, 
where  vg  is  the  group  velocity  in  the  waveguides  (for  signal  and  control  photons) 
and  f  the  emitter  decay  rate  into  waveguide  modes. 
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Figure  3.7.2  Gate  fidelity  simulation  results. 

We  plan  to  expand  this  analysis  to  time-dependent  waveguide-cavity  coupling, 
which  promises  to  actually  allow  near-unity  gate  fidelity. 


4.  Future  Directions  in  Optical  Quantum 
Computing  and  Networking 


This  seedling  program  has  improved  devices  and  developed  theoretical  protocols 
that  greatly  lowered  the  resource  requirements  of  linear  optics  quantum  computing 
and  all-optical  quantum  repeaters.  This  section  considers  promising  future 
directions. 

4.1  Photonic  nanocavity  design  to  achieve 
room-temperature  single-photon  strong  coupling 
with  dielectric  materials 

We  developed  a  photonic  crystal  nanocavity  design  with  arbitrarily  small  mode 
volume(Veff),  which  allows  nonlinear  optical  response  at  the  single-photon  level. 
The  design  relies  on  enhancement  based  on  two  electric  field  boundary  conditions. 
These  methods  can  be  concatenated  to  further  reduce  Veff.  Depending  on  the 
termination,  both  dielectric  and  air  cavity  can  be  designed  for  a  target  application. 
We  illustrate  the  nanocavity  concept  in  a  silicon-air  nanobeam  at  1550  nm  that 
reaches  an  ultrasmall  volume  Veff~2*10  /n)3.  With  a  suitable  chi-3  medium,  this 

cavity  allows  single-photon  nonlinearity  —  even  at  room  temperature. 


Fig  4.1.1  (a)  Geometry  of  base  structure  (b)~(e)  Electric  field  distribution,  (inset) 
Electric  energy  density,  (b)  Type  1  enhancement  (slot  cavity),  (c)  Type  2 
enhancement  with  narrow  bridge  (5nm).  (d)  Wide  bridge  (40nm).  (e)  intermediate 
size  bridge  (20nm). 


Optical  nanocavities  with  small  mode  volume  ( V eff)  and  high  quality 
factor  (Q)  can  greatly  increase  the  light-matter  interaction  and  have  a  wide  range 
of  application  including  nanocavity  lasers,  cavity  quantum  electrodynamics 
(cQED),  and  nonlinear  optics.  Photonic  crystal  (PhC)  cavities  with  an  air  slot  can 
achieve  exceptionally  small  V eff  on  the  order  of  0.01  X3 ,  where  X  is  the 
free-space  wavelength[17,18].  This  cavity  relied  on  enhancement  through  a 
boundary  condition  on  normal  electric  displacement.  We  demonstrated  further 
reduction  of  Veff  using  the  second  boundary  condition  on  parallel  electric  field. 


Furthermore,  concatenation  of  both  enhancement  enables  arbitrarily  small 
limited  only  by  practical  considerations  such  as  fabrication  resolution.  The 
extreme  field  confinement  of  this  cavity  design  opens  new  possibilities  in 
nonlinear  optics;  in  particular,  such  cavity  designs  can  enable  Kerr-nonlinear 
optical  phenomena  at  the  single  photon  level. 

Typical  photonic  crystal  cavity  obtains  field  confinement  through 
interference  that  limits  the  size  of  confinement  as  half  wavelength  (  Veff~  QJlrif 

):  diffraction  limit.  However,  V eff,  which  is  local  in  the  sense  that  it  is  determined 
by  field  at  one  position  (maximum  field),  is  not  strictly  limited  by  the  diffraction 
limit  if  field  confinement  does  not  come  from  interference:  strong  inhomogeneity 
of  medium.  Mode  volume  reduction  beyond  the  diffraction  limit  has  been  studied 
with  slot  cavities  by  Lipson  and  coworkers  [19,20].  In  a  slot  cavity,  field  in  low 
index  region  is  enhanced  through  boundary  condition  on  electric  displacement: 

£h^h  =  Dh-L  =  D[j_  =  &lEl ...  (4. 1 .2) 

Ei=%Eh...{AA.2>) 

where  subscript  I  and  h  mean  low  and  high  index  respectively,  and  -1-  means 
normal  component  of  field  with  respect  to  boundary.  Maximum  electric  energy 
density  is  increased  by  a  factor  of: 


Eel=lEl 

^eo  e„\Eh\2 


? -(4.1-4) 


If  we  assume  that  a  slot  is  narrow  that  the  numerator  of  (1 )  is  not  changed  by  a 
slot,  then  Veff  is  reduced  by  a  factor  of  zh/zl .  For  a  silicon-air  nanobeam  cavity, 

this  enhancement  is  eh/et~  13.9  and  mode  volume- 0.0R 3  was  achieved  [17,18]. 
Field  enhancement  through  normal  electric  displacement  continuity  (Type  1 
enhancement)  is  advantageous  in  that  it  is  inherently  wavelength  independent 
endowing  more  tolerance  to  errors  occured  from  simulation  and  fabrication. 
However,  since  a  cavity  has  high  electric  field  in  low  index  region  (air  cavity), 
there  has  been  limited  applications. 

Here,  we  propose  a  method  to  further  reduce  Veff  with  parallel  electric 
field  continuity  across  a  boundary  (Type  2  enhancement).  More  specifically,  this 
is  achieved  by  a  high  index  bridge  across  a  slot  (Fig.  1(a)).  Parallel  electric  field 
across  the  boundary  is  given  by  (5): 


Eh\\  =  Ei\\  .5) 


Through  Type  2  enhancement,  large  electric  field  in  low  index  region  (slot)  is 
transferred  into  high  index  region  (bridge).  As  a  result,  the  bridge  has  the  largest 
electric  energy  density;  cavity  is  a  dielectric  cavity.  This  cavity  is  useful  when  one 
needs  light-matter  interaction  in  high  index  material.  Maximum  electric  energy 
density  enhancement  is  given  by  (4.1.6): 


=  £#/ 
We  1 


6/ 


(4.1.6) 


In  the  same  logic  with  the  Type  1 ,  if  a  bridge  is  narrow  enough  that  adding  a 
bridge  out  of  a  slot  does  not  change  the  numerator  of  (1 ),  then  Veff\s  reduced  by 
a  factor  of  sh/s, . 

To  demonstrate  our  assertion,  3D  finite-difference  time-domain  (FDTD) 
simulation  was  conducted  on  the  nanobeam  cavities  depicted  in  figure  1(a)  with 
variations  on  the  dashed  box.  Simulated  Veff  is  9.48  x  10“21(2.54x  10“3(A,)3), 
which  is  ~10  times  smaller  than  the  V eff  of  slot  cavity 
(9.47  x  10~20,  2.51  x  10~2(A,)3) .  This  enhancement  is  slightly  smaller  than 
eh/sl  =  12.09  because  of  the  finite  size  of  bridge  and  electric  field  components 
perpendicular  to  the  boundary.  "Enhancement  by  boundary  condition"  only 
happens  when  one  boundary  condition  is  dominant  over  the  other.  In  a  slot  cavity 
(figure  4.1.2(b)),  lateral  boundary  between  low  index  and  high  index  should  be 
dominant  over  the  vertical  boundary  between  slot  and  holes.  For  the  field  in  the 
bridge  to  be  forced  same  as  the  field  in  slot,  vertical  boundary  between  high 
index  region  and  low  index  region  should  be  dominant  over  the  lateral  boundary 
between  the  bridge  and  high  index  slab.  Figure  4.1.2(d)  and  4.1.2(e)  illustrate 
simulation  results  when  the  previously  stated  requirement  is  not  satisfied.  For  a 
wide  or  intermediate  width  bridge  case,  field  in  the  bridge  is  not  same  with  that  in 
the  slot  because  effects  from  boundary  between  bridge  and  slab  is  not  negligible. 
However,  there  is  still  enhancement. 


Fig  4.1.2.  Electric  field  distribution,  (a)  2nd  level  slot  (SBS)  (b)  Tip  shape  as  a 
limiting  case  of  concatenation.  Here,  we  show  air  cavity  limit  of  a  tip.  (c)  Quality 
factors  and  mode  volume  in  each  level,  (d)  Two  figure  of  merit:  Q/V  as  a  general 
criteria  and  QVM/V eff  for  strong  coupling. 


We  furthermore  find  that  the  extreme  light  confinement  allows  us  to  reach  the 
strong  coupling  regime:  a  single  photon  injected  into  the  cavity  produces  a 
sufficient  index  shift  (by  Kerr  nonlinearity)  that  the  cavity  frequency  shifts  by  more 
than  one  linewidth,  assuming  a  Q  factor  >  100,000  and  materials  inside  the  cavity 
gap  such  as  Graphene  (0.5  *  10 “13m2/F2  Zhang  et  al.  2012).  Indium  Tin  oxide  ( 

0.5  *  10"16m2/F2[21]),  and  J-aggregate  (0.5  *  10 ~l3m2/V2  Zhuravlev  et  al.  19921. 


4.2  All-Optical  Ising  Machine 

Devices  that  can  solve  the  Ising  model  via  quantum  annealing,  a  restricted  class  of 
adiabatic  quantum  computers  have  recently  gained  popularity,  since  the  Ising  model 
can  encode  many  interesting  NP  hard  optimization  problems.  It  is  however  not  yet 
clear  if  quantum  annealing  provides  speed  up  in  terms  of  how  the  number  of 
computational  steps  scale  as  the  size  of  the  problem  grows  as  compared  to  the  best 
classical  algorithms.  All  quantum  annealing  architectures  realized  to  date  employ 
'matter'  qubits,  such  as  superconducting,  NMR  and  Bose-Einstein  condensates.  Even 
in  the  absence  of  a  computational  complexity  advantage,  an  all-optical  quantum 
annealer  could  be  very  attractive  both  due  to  its  ability  to  gain  a  large 
constant-factor  speedup,  as  well  as  achieve  a  potentially  large  power  saving,  over 


Mode  Volume  (nT 


conventional  electronic-computing  solutions.  We  designed  an  all-optical  realization 
of  quantum  annealing  and  an  explicit  integrated  photonic  architecture  for  its 
implementation.  The  proposed  architecture  employs  an  array  of  spectral-dual-rail 
encoded  optical  qubits  in  a  single  pair  of  coupled  non-linear  cavities,  and 
dynamically-tunable  linear-optic  and  cross-phase  interactions  across  pairs  of  qubits, 
with  very  small  yet  tunable  non-linear  phases.  Our  calculations  show  that  the 
maximum  nonlinearities  required  in  the  protocol  are  of  the  order  of  10'4.  Such 
nonlinearities  would  be  easily  achievable  using  the  fractal  cavities  described  in  the 
previous  section. 
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Rate-distance  tradeoff  and  resource  costs  for  all-optical  quantum  repeaters 
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10  Moulton  Street,  Cambridge,  MA  02138,  USA 

We  present  a  resource-performance  tradeoff  of  an  all-optical  quantum  repeater  that  uses  photon 
sources,  linear  optics,  photon  detectors  and  classical  feedforward  at  each  repeater  node,  but  no 
quantum  memories.  We  show  that  the  quantum-secure  key  rate  has  the  form  R(rj)  =  Dr]s  bits  per 
mode,  where  77  is  the  end-to-end  channel’s  transmissivity,  and  the  constants  D  and  s  are  functions 
of  various  device  inefficiencies  and  the  resource  constraint,  such  as  the  number  of  available  photon 
sources  at  each  repeater  node.  Even  with  lossy  devices,  we  show  that  it  is  possible  to  attain 
s  <  1,  and  in  turn  outperform  the  maximum  key  rate  attainable  without  quantum  repeaters, 
^direct  (tj)  =  —  log2(l  —  rj)  ~  (1/ In  2)77  bits  per  mode  for  77  <C  1,  beyond  a  certain  total  range  L, 
where  77  ~  e~aL  in  optical  fiber.  We  also  propose  a  suite  of  modifications  to  a  recently-proposed  all- 
optical  repeater  protocol  that  ours  builds  upon,  which  lower  the  number  of  photon  sources  required 
to  create  photonic  clusters  at  the  repeaters  so  as  to  outperform  indirect  (77),  from  ~  1011  to  ~  106 
photon  sources  per  repeater  node.  We  show  that  the  optimum  separation  between  repeater  nodes  is 
independent  of  the  total  range  L,  and  is  around  1.5  km  for  assumptions  we  make  on  various  device 
losses. 

PACS  numbers:  42.50.Ex,  03.67.Dd,  03.67.Lx,  42.50.Dv 


I.  INTRODUCTION 

Quantum  key  distribution  (QKD)  enables  two  distant 
authenticated  parties  Alice  and  Bob,  connected  via  a 
quantum  (e.g.,  optical)  channel,  to  generate  information- 
theoretically  secure  shared  secret  bits.  No  knowledge  of 
the  channel  conditions  (noise  model,  or  any  channel  es¬ 
timate)  is  required  a  priori  to  ensure  security.  However, 
the  shared  secret  is  generated  at  a  rate  commensurate 
with  the  worst-case  adversary  physically  consistent  with 
the  channel  conditions  actually  presented  to  Alice  and 
Bob.  The  reason  is  that  all  the  perceived  channel  im¬ 
perfections  (absolutely  anything  that  causes  the  channel 
map  to  deviate  from  a  noiseless  identity  transformation) 
is  attributed  to  the  actions  of  the  most  powerful  adver¬ 
sary  allowed  by  physics — even  though  some  (or  all)  of 
that  deviation  of  the  channel  from  an  identity  map  may 
actually  stem  from  non- adversarial  sources,  such  as  losses 
due  to  free-space  diffraction,  fiber  loss,  detection  ineffi¬ 
ciency,  thermal  noise  from  blackbody  at  the  operating 
temperature  and  wavelength,  and  detector  noise.  An  im¬ 
portant  consequence  of  this  assumption  is  that  all  the 
signal  power  transmitted  by  Alice  that  is  not  collected 
by  Bob  is  made  available  coherently  to  the  eavesdropper, 
Eve.  This  model  for  Eve  is  the  intuition  behind  why  the 
secret  key  rate  for  a  direct-transmission  based  QKD  pro¬ 
tocol  must  decrease  linearly  with  77,  the  Alice-Bob  power 
transmissivity,  in  the  77  1  regime  EEt  For  any  direct- 

transmission  protocol  over  the  pure-loss  optical  channel 
of  transmissivity  77,  and  assuming  unlimited  authenti¬ 
cated  two-way  public  classical  communication,  it  was  re¬ 
cently  shown  that  the  key  rate  cannot  exceed  —  log2(l— 77) 
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bits  per  mode  [2],  which  is  «  I.4477  for  77  1.  For 

a  pure-loss  channel,  the  Pirandola-Laurenza-Ottaviani- 
Banchi  (PL OB)  upper  bound  improves  over  the  Takeoka- 
Guha-Wilde  (TGW)  bound  [QQ  by  a  factor  of  2  in  the 
77  <C  1  regime.  The  TGW  bound  is  an  upper  bound  on 
the  secret-key  agreement  capacity  with  unlimited  two- 
way  classical  communication  P2  ( Af ) ,  applicable  to  a  gen¬ 
eral  quantum  channel  Af.  For  the  pure-loss  channel  Afv, 
the  PLOB  bound  coincides  with  the  best-known  achiev¬ 
able  rate  [3],  thus  establishing  P2  (Af^ )  =  —  log2(l  —  77) 
bits  per  mode.  From  hereon,  we  denote  by  indirect  (v)  = 
—  log2(l  —  77)  the  maximum  bits-per-mode  secret  key 
rate  achievable  by  any  direct-transmission  QKD  protocol, 
i.e.,  without  the  use  of  quantum  repeaters.  The  bits/s 
rate  of  a  QKD  protocol’s  implementation  is  obtained  by 
multiplying  the  bits/mode  rate  by  the  spatio-temporal- 
polarization  bandwidth  (modes/s),  which  is  governed  by 
the  channel  geometry,  and  the  transmitter  and  detector 
bandwidth.  Since  loss  increases  exponentially  with  dis¬ 
tance  L  in  optical  fiber  (i.e.,  77  =  e~aL ),  for  77  1, 

the  key  rate  generated  by  any  direct-transmission  QKD 
protocol  must  decay  exponentially  with  the  range  L.  Ex¬ 
pressed  as  a  function  of  L,  indirect  (A)  =  —  log2(l  —  e~aL) 
bits/mode,  which  is  ~  lA4e~aL  bits/mode,  for  L  large. 

Quantum  repeaters,  proposed  in  [4],  are  devices  which 
when  inserted  along  the  length  of  the  optical  chan¬ 
nel,  can  help  generate  shared  secret  at  a  rate  that  sur¬ 
passes  indirect  (v)  at  any  value  of  Alice-to-Bob  channel 
transmissivity  77  [2].  Quantum  repeaters  need  not  be 
trusted  or  physically  secured  in  order  to  ensure  the  se¬ 
curity  of  the  keys  generated.  If  n  quantum  repeaters 
are  inserted  along  the  length  of  the  channel  connect¬ 
ing  the  communicating  parties  Alice  and  Bob,  and  if 
there  are  absolutely  no  physical  constraints  placed  on 
the  repeater  nodes  (i.e.,  the  repeaters  are  assumed  to  be 
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lossless,  error-corrected,  general  purpose  quantum  com¬ 
puters),  then  the  maximum  key  rate  achievable  by  Alice 
and  Bob  is  given  by  —  log2(l  —  ^min)  bits/mode,  where 
Vmin  =  min  (?7i ,  772, . . . ,  Vn+i),  With  rj  =  r)i...T}n+1,  is 
the  transmissivity  of  the  lossiest  link  between  successive 
repeater  nodes  [5]  (see  [6]  for  a  different  upper  bound 
based  on  squashed  entanglement  HD-  Given  n  ideal  re¬ 
peater  nodes,  their  optimal  placement  is  to  lay  them 
equally-spaced,  in  which  case,  the  maximum  achievable 
rate  is  —  log2(l  —  77 VO+1))  bits/mode.  As  n  00,  the 
rate  is  unbounded.  However,  assuming  repeaters  to  be 
lossless  error-corrected  quantum  computers  is  not  practi¬ 
cal.  A  more  practically  relevant  question  to  ask  is  if  the 
repeater  nodes  have  finite  resources  with  lossy  and  im¬ 
perfect  components  (where  ‘resources’  may  be  different 
physical  entities  depending  upon  the  type  of  quantum  re¬ 
peater  and  the  protocol  employed),  then  what  rate  can 
Alice  and  Bob  achieve,  and  more  importantly  what  would 
it  take  to  build  repeater  nodes  so  as  to  be  able  to  signif¬ 
icantly  outperform  indirect  (jj)  =  ~  l°g2(l  —  v)  bits/mode. 
This  is  the  topic  addressed  in  this  paper,  for  repeaters 
that  are  built  solely  using  photonic  components — single¬ 
photon  sources,  detectors,  electro-optic  feedforward,  but 
no  matter-based  quantum  memories.  As  we  will  see  later 
in  this  paper,  that  given  physical  constraints  on  a  re¬ 
peater  node,  placing  more  repeaters  (higher  n)  between 
Alice  and  Bob  may  not  always  improve  the  rate,  i.e.,  de¬ 
pending  upon  the  total  distance  L  (or  equivalently,  the 
transmissivity  77)  between  Alice  and  Bob,  and  given  the 
physical  device  constraints  in  a  repeater  node,  there  may 
be  an  optimal  number  77*  (77)  of  nodes,  which  achieves  the 
highest  end-to-end  rate. 

Traditional  quantum  repeaters  work  in  the  framework 
of  entanglement-based  QKD.  At  the  end  of  the  trans¬ 
mission  phase  of  an  entanglement-based  QKD  proto¬ 
col,  Alice  and  Bob  share  (noisy,  or  imperfect)  entangled 
pairs  (e.g.,  of  photons  or  matter-based  stationary  qubits) 
which  could  have  been  tampered  with,  or  which  could 
have  deteriorated  due  to  channel  loss  and  noise.  At  that 
point,  if  the  end  goal  of  Alice  and  Bob  is  to  generate 
shared  entanglement  (for  use  in  some  quantum  protocol 
that  consumes  shared  entanglement,  such  as  telepora- 
tion  [7]  or  dense  coding  [8]),  they  would  perform  entan¬ 
glement  distillation  to  sieve  out  a  small  number  of  clean 
maximally  entangled  Bell  pairs  by  performing  local  op¬ 
erations  and  classical  communications  (LOCC).  If  the 
end  goal  of  Alice  and  Bob  is  to  generate  shared  secret 
(a  strictly  less  demanding  goal  than  generating  shared 
entanglement),  they  directly  measure  the  noisy  shared 
entangled  pairs,  and  perform  (classical)  error  correction 
and  privacy  amplification  on  their  correlated  measure¬ 
ment  results  over  an  authenticated  public  channel  to  dis¬ 
till  a  quantum-secure  shared  secret  key. 

Several  different  genres  of  repeater  protocols  have  been 
proposed  [9].  The  two  primary  ingredients  in  any  of  the 
traditional  repeater  architectures  are:  (1)  some  form  of 
a  quantum  memory,  and  (2)  the  ability  to  perform  a  cer¬ 
tain  restricted  class  of  quantum  logic,  i.e.,  gates  and  mea¬ 


surements  on  the  flying  (photonic)  qubits  as  well  as  the 
static  (memory)  qubits.  In  the  most  basic  repeater  pro¬ 
tocol,  the  restricted  quantum  operation  required  is  Bell 
state  measurement  (BSM)  on  pairs  of  qubits.  A  BSM  on 
qubit  b  and  qubit  c  converts  two  independent  Bell  pairs 
|T)a6  and  |T)cd  into  one  Bell  pair  |T)ad,  upto  local  sin¬ 
gle  qubit  operations,  a  process  known  as  entanglement 
swapping. 


A.  Quantum  repeaters  based  on  mode  multiplexing 
and  Bell  state  measurements 

In  the  following  discussion,  we  will  focus  on  a  class 
of  quantum  repeaters  that  rely  solely  on  probabilistic 
BSMs,  quantum  memories,  and  multiplexing,  i.e.,  the 
ability  to  ‘switch’  qubits  across  (spatial,  spectral,  or  tem¬ 
poral)  modes.  The  essence  of  such  a  repeater  protocol 
was  developed  by  Sinclair  et  al.  pm,  which  employed 
spectral  multiplexing  in  multimode  quantum  memories 
across  m  parallel  (spectral)  channels,  and  entanglement 
swapping  using  linear  optics  and  single  photon  detectors 
(the  success  probability  of  which  can  at  most  be  50%). 
Guha  et  al.  analyzed  the  secret  key  rates  achievable  by 
the  above  protocol,  with  a  fixed  m  (memory  size)  and 
found  that  even  when  photon  loss  is  the  only  source  of 
noise,  the  achievable  key  rate  is  of  the  form  R(rj)  =  Dr]s 
bits/mode,  where  D,  and  s  <  1  are  constants  that  are 
functions  of  various  losses  in  the  system  (e.g.,  detection 
efficiencies,  coupling  losses,  memory  loading  and  readout 
efficiencies,  and  BSM  failure  probability)  [Tl] .  Since  the 
exponent  of  77,  i.e.,  s  is  strictly  less  than  1,  the  key  rate 
must  beat  indirect  (??)  (which  scales  as:  oc  77  for  77  <C  1) 
beyond  a  certain  minimum  distance  determined  by  the 
actual  values  of  the  system’s  loss  parameters,  which  is 
around  a  couple  of  hundred  kilometers  for  reasonable  es¬ 
timates  of  the  losses  m-  Since  77  =  e  aL  in  fiber,  the 
rate  achieved  by  this  repeater  protocol  for  a  fixed  mem¬ 
ory  size,  R(L)  =  De~saL  still  scales  exponentially  with 
the  range  L,  albeit  with  a  smaller  exponent  compared  to 
the  best  possible  rate  without  any  repeater,  which  could 
turn  into  a  huge  absolute  improvement  in  the  end-to-end 
secret  key  rate  m- 

Azuma  et  al.  recently  proposed  an  all-photonic  variant 
of  this  protocol  in  which  they  substituted  matter  based 
quantum  memories  with  optical  cluster  states  E2,  based 
on  a  proposal  by  Varnava  et  al.  to  mimic  a  quantum 
memory  (i.e.,  protect  against  photon  losses)  by  append¬ 
ing  each  physical  photonic  qubit  by  an  entangled  ‘tree 
cluster’  state  m ■  As  long  as  the  losses  incurred  by  each 
photon  (i.e.,  photons  being  protected  as  well  as  the  ad¬ 
ditional  photons  in  the  trees  added  for  loss  protection) 
is  less  than  3  dB,  the  effective  loss  of  the  logical  qubit 
can  be  made  to  approach  zero,  by  increasing  the  size  of 
the  tree  cluster,  i.e.,  the  number  of  photons  in  the  logi¬ 
cal  qubit  [14] .  Thus,  Azuma  et  aids  proposal  showed  the 
theoretical  feasibility  of  a  quantum  repeater  architecture 
(i.e.,  one  that  can  beat  the  scaling  of  direct-transmission 
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QKD)  using  only  flying  qubits,  with  the  repeater  nodes 
being  equipped  only  with  single  photon  sources,  passive 
linear-optical  circuits  (beamsplitters  and  phase  shifters), 
single  photon  detectors,  and  classical  feedforward. 

Azuma  et  aids  result  marked  a  promising  conceptual 
leap  towards  all-optical  quantum  repeaters.  However, 
important  unanswered  questions  remained,  including  the 
achievable  secure  key  generation  rate  and  how  it  scales 
with  distance  (or  loss),  as  well  as  the  physical  resource 
requirements:  e.g.,  the  number  of  photon  sources  and  de¬ 
tectors  at  the  repeater  nodes.  As  an  example,  a  calcula¬ 
tion  in  their  paper  shows  that  at  a  range  of  L  =  5000  km, 
an  entanglement-generation  rate  of  69  kHz  is  achievable 
in  a  fiber  based  linear  optic  system  with  100  kHz  repeti¬ 
tion  rate,  150  ns  feed  forward  time  and  a  source-detector 
efficiency  product  of  95%  whereas  sharing  a  single  entan¬ 
gled  photon  pair  via  a  direct  transmission  scheme  with 
the  same  parameters  would  require  1081  years.  The  level 
of  error  protection  required  to  achieve  the  aforesaid  re¬ 
peater  performance  at  L  =  5000  km  would  require  one  to 
build  entangled  clusters  of  ~  104  photons  at  the  100  kHz 
clock  rate  at  each  repeater  node.  Building  such  a  cluster 
using  linear  optics  and  feed-forward  G51QU  would  require 
around  1024  single  photon  sources  at  each  repeater  node. 
Furthermore,  since  every  photon  used  for  error  correc¬ 
tion  is  sent  between  repeater  nodes  in  [12],  their  scheme 
would  require  around  20,  000  parallel  channels  connecting 
the  neighboring  nodes.  Thus,  while  Ref.  [12]  showed  the 
theoretical  possibility  of  all-optical  repeaters,  clearly  fur¬ 
ther  work  is  needed  to  address  their  practical  feasibility. 
These  results  open  up  a  compelling  line  of  research  to  in¬ 
vestigate  improved  all-photonic  repeater  architectures  of 
various  genres  which  could  be  built  with  practically  fea¬ 
sible  resources,  and  also  a  thorough  comparative  study 
of  rates  achievable  with  each  such  all-optical  repeater 
scheme. 


B.  Main  results 

Our  contributions  in  this  paper  are  twofold.  The  first 
is  a  rigorous  analysis  of:  (a)  the  secret  key  rates  achiev¬ 
able  with  the  the  aforesaid  all-photonic  repeater  archi¬ 
tecture  given  the  size  of  the  clusters  generated  at  each 
repeater  station,  and  (b)  the  resources  required  (e.g., 
number  of  single  photon  sources  and  detectors  required 
at  each  repeater  node)  to  build  that  cluster,  while  tak¬ 
ing  into  account  in  explicit  detail  each  step  in  building 
the  required  clusters  using  a  network  of  passive  linear  op¬ 
tics  (i.e.,  beamsplitters  and  phase  shifters),  imperfect  on- 
demand  sources  with  loss  (see  section  [II|  for  a  description 
of  the  source),  single  photon  detectors  (with  some  num¬ 
ber  resolving  capability),  and  feed-forward.  We  find  that 
the  achievable  secret  key  rate  scales  as  Dr]s  bits/mode, 
where  D  and  s  <  1  are  functions  of  the  number  of  photon 
sources  at  each  repeater  node  (the  resource  constraint — 
which  is  parametrically  related  to  the  size  of  the  cluster) , 
all  the  ‘inline’  losses  (e.g.,  losses  in  the  optical  fiber  or 


waveguide  used  while  creating  the  cluster,  independent  of 
the  fiber  loss  between  repeater  stations),  and  the  source 
and  detector  efficiencies.  With  77  ~  e~aL  in  fiber,  the  key 
rate  still  scales  exponentially  with  L,  but  with  a  smaller 
exponent  compared  to  the  best  direct-transmission  pro¬ 
tocol.  This  is  no  surprise  given  the  analysis  of  HD  ,  since 
the  tree-cluster  construction  of  [12]  essentially  mimics 
an  imperfect  quantum  memory,  but  one  whose  efficiency 
cannot  simply  be  modeled  by  a  constant  per  mode  as  in 
Ref.  m-  Using  the  cluster  building  scheme  proposed  by 
Li  et  al.  [16],  we  find  that  to  a  good  approximation,  the 
resource  requirements  are  determined  by  the  number  of 
probabilistic  fusion  steps  k  required  to  build  the  cluster 
starting  from  single  photons,  and  hence,  we  calculate  the 
performance  with  the  best  cluster  that  can  be  built  in  k 
fusion  steps.  We  use  the  scheme  of  Li  et  al.  because  it 
has  been  shown  to  be  more  efficient  than  the  scheme  of 
Varnava  et  al.  [15]  at  building  clusters  [16].  Given  all  the 
inline  and  device  losses,  we  evaluate  the  number  of  pho¬ 
ton  sources  (and  detectors)  needed  at  each  repeater  node 
to  beat  indirect  (A)  at  a  given  total  range  L  between  Al¬ 
ice  and  Bob.  We  also  prove  that  given  the  device  losses, 
there  is  an  optimal  spacing  between  the  repeater  nodes 
(which  evaluates  to  roughly  1.5  km  for  a  set  of  system 
parameters  we  choose),  regardless  of  the  overall  range  L. 

Our  second  major  contribution  in  this  paper  is  a  sig¬ 
nificant  improvement  to  the  all-photonic  repeater  archi¬ 
tecture  in  12 — both  in  terms  of  the  resources  required 
at  each  node  and  the  number  of  parallel  optical  channels 
connecting  the  neighboring  nodes.  We  find  that  barely 
beating  indirect  (A)  using  the  all-optical  scheme  of  [12]  re¬ 
quires  more  than  1011  photon  sources  at  each  repeater 
node  for  realizing  the  required  optical  cluster  states  and 
measurements.  It  also  requires  208  parallel  channels  con¬ 
necting  neighboring  nodes,  even  when  assuming  very  op¬ 
timistic  device-loss  parameters.  Assuming  the  same  de¬ 
vice  losses,  our  improved  repeater  architecture  reduces 
the  number  of  photon  sources  (to  barely  beat  indirect  (A)) 
by  5  orders  of  magnitude,  while  reducing  the  number 
of  parallel  channels  to  8.  In  both  of  these  calculations, 
each  source  is  used  only  once  per  clock  cycle,  i.e.,  they 
are  not  temporally  multiplexed.  We  prove  a  tight  ana¬ 
lytical  lower  bound  for  the  performance  of  our  improved 
scheme.  These  performance  advances  are  enabled  pri¬ 
marily  by  the  following:  (1)  using  boosted  fusion  logic 
that  improves  the  success  probability  of  the  BSM  to  75% 
by  using  four  ancilla  single  photons  [17,,  (2)  employing  a 
more  resource-efficient  scheme  for  creating  tree  clusters, 
building  on  the  work  of  [15j  [16],  (3)  retaining  all  the 
ancilla  photons  used  for  loss  protection  (i.e.,  to  mimic 
a  quantum  memory)  locally  at  the  repeater  nodes  in  a 
lossy  waveguide,  and  (4)  optimizing  the  timing  of  several 
single  qubit  measurements  in  the  entire  protocol. 

We  will  limit  our  analysis  to  include  photon  losses  (dur¬ 
ing  the  entire  ‘lifetime’  of  each  photon,  i.e.,  from  the  time 
of  generation  to  detection)  but  will  not  consider  ‘multi¬ 
photon’  errors  stemming,  for  instance,  from  multi-photon 
emissions  from  the  source,  or  detector  dark  clicks.  We 
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should  note  however  that  the  error  correction  scheme  an¬ 
alyzed  here  also  provides  some  protection  against  depo¬ 
larizing  noise  [12],  a  variant  of  which  arises  when  one 
assumes  multi-photon  errors,  and  errors  stemming  from 
imperfect  mode  matching  within  the  passive  linear  opti¬ 
cal  circuits  at  the  repeater  nodes. 

The  remainder  of  the  paper  is  organized  as  follows. 
Section  [TT|  reviews  preliminaries  and  notation  used  in  the 
paper.  Section  III  describes  our  (improved)  all-photonic 
quantum  repeater  architecture  with  a  detailed  descrip¬ 
tion  of  each  step  starting  from  the  creation  of  the  tree 
clusters  for  error-protection,  photon  transmission,  mea¬ 
surements  at  the  repeater  nodes,  and  the  measurements 
by  Alice  and  Bob,  followed  finally  by  key  generation.  Sec¬ 
tion  |IV]  derives  a  closed  form  expression  for  a  lower  bound 
to  the  rate-distance  envelope  (i.e.,  an  achievable  rate  by 
the  protocol),  which  we  show  (numerically)  to  match  the 
true  rate-distance  envelope  extremely  closely.  Section  [V] 
compares  our  scheme  to  that  of  Ref.  [12|  in  terms  of 
resource  requirements  and  rates,  and  discusses  possible 
avenues  for  further  improvement.  The  concluding  sec¬ 
tion  [VT]  provides  concrete  directions  for  future  research 
in  order  to  further  improve  the  prospects  of  a  quantum 
communications  network  based  solely  on  flying  qubits. 


II.  PRELIMINARIES 

In  this  paper,  we  work  with  dual-rail  photonic  qubits, 
where  the  logical  |0)  and  |1)  are  encoded  by  a  single  pho¬ 
ton  in  one  of  two  orthogonal  (spatial)  modes.  A  photonic 
cluster  state  (or,  graph  state),  on  a  graph  G(V,E)  with 
vertices  in  set  V  and  edges  in  set  E,  can  be  constructed 
by  preparing  each  of  the  \V\  qubits  (one  stationed  at 
each  vertex)  in  the  state  (|0)  +  |1))/a/2,  and  applying 
\E\  controlled-phase  operations  (a  two-qubit  unitary  gate 
that  applies  a  pauli  Z  gate  to  the  second  qubit  if  the 
first  qubit  is  in  the  |1)  state  and  applies  an  identity  oth¬ 
erwise)  on  each  pair  of  vertices  that  share  an  edge  m ■ 
The  (entangled)  quantum  state  of  the  | V|  qubits  thus 
obtained  is  an  eigenstate  of  the  \V\  stabilizer  operators 
Xi  njeAr{i)Zj,  where  the  index  i  runs  over  all  the  ver¬ 
tices,  Xi  and  Zj  are  Pauli  X  and  Z  operators  on  qubit  i 
and  qubit  j  respectively,  and  A f(i)  is  the  set  of  all  near¬ 
est  neighbor  vertices  of  vertex  i.  One  simple  observa¬ 
tion,  given  that  the  cluster  state  is  an  eigenstate  of  the 
aforesaid  stabilizer  operators,  is  that  an  X  measurement 
on  qubit  i,  and  Z  measurements  on  all  but  one  of  the 
qubits  in  would  deterministically  reveal  what  the 

outcome  of  a  Z  measurement  on  that  unmeasured  qubit 
in  A f{i)  would  have  been,  even  if  that  unmeasured  qubit 
had  been  lost.  This  realization  is  at  the  heart  of  the 
tree-based  counterfactual  error  correction  for  protection 
against  photon  losses,  developed  by  Varnava  et  al.  EH- 
The  idea  is  to  attach  a  tree  cluster  to  each  physical  pho¬ 
tonic  qubit  in  the  graph  state  that  needs  to  be  protected 
against  qubit  loss.  One  can  then  deduce  the  result  of  any 
measurement  on  that  qubit  via  an  appropriate  sequence 


of  measurements  on  the  qubits  of  the  attached  tree.  The 
physical  qubit  and  the  qubits  of  the  tree  together  form 
a  protected  (logical)  qubit.  We  consider  regular  trees 
described  by  the  branching  vector  b  =  {60,  b\ , . . . ,  6m}, 
which  signifies  that  the  root  of  the  tree  has  bo  children 
nodes,  and  each  of  those  nodes  have  b\  children  nodes, 
and  so  on  until  bobi  . .  .bm  nodes  at  depth  m.  For  such 
regular  trees  used  for  loss-error  protection,  one  can  write 
an  explicit,  yet  recursive,  expression  for  the  success  prob¬ 
ability  P  of  performing  an  arbitrary  single-qubit  mea¬ 
surement  on  the  protected  qubit  m  It  was  shown  that 
one  can  push  P  arbitrarily  close  to  1  as  long  as  the  prob¬ 
ability  of  losing  each  photon  is  less  than  1/2.  Fig.  [l] 
illustrates  how  to  attach  a  {3,  2,  2}  tree,  shown  by  the 
dark  (purple)  shaded  nodes,  to  a  physical  qubit  of  a  clus¬ 
ter,  shown  by  light  (green)  shaded  nodes.  Note  that  after 
the  tree  cluster  is  attached  to  the  physical  qubit,  X  basis 
measurements  must  be  performed  on  the  physical  qubit 
itself  and  the  root  node  of  the  tree.  These  X  basis  mea¬ 
surements,  if  successful,  create  additional  edges  (shown 
in  dashed  blue  in  Fig.  [l])  between  each  neighboring  qubit 
of  the  root  node  and  each  neighboring  qubit  of  the  phys¬ 
ical  qubit,  after  which  the  tree-protected  logical  qubit  is 
ready  to  use. 


FIG.  1.  Attaching  a  {3,2,2}  tree  to  a  node  of  a  photonic 
cluster. 

An  ideal  loss-less  photonic  cluster  state  on  graph  G 
is  a  pure  state,  \^)G-  A  lossy  cluster  state  on  G  is  ob¬ 
tained  when  all  the  photonic  qubits  of  \if)G  are  transmit¬ 
ted  through  independent  pure-loss  beamsplitters  each  of 
transmissivity  r\.  We  call  1  —  77  the  loss  rate  of  such  a 
lossy  cluster  state.  Clearly,  the  loss  rate  of  1^)^  itself  is 
0. 

Arbitrary  photonic  cluster  states  can  be  prepared — 
with  non- unity  probability — using  ideal  single  photons, 
passive  linear  optics  (i.e.,  beamsplitters  and  phase 
shifters)  and  single  photon  detectors  [19].  As  examples, 
in  the  absence  of  losses,  a  two-photon  maximally  entan¬ 
gled  (Bell)  state  can  be  prepared  with  success  probability 
3/16  [20].  whereas  a  three-photon  maximally  entangled 
(GHZ)  state  can  be  prepared  with  success  probability 
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1/32  H5].  Browne  and  Rudolph  introduced  linear-optical 
Type  I  and  Type  II  two-qubit  fusion  gates,  which  if  suc¬ 
cessful  (with  probability  1/2),  can  fuse  two  cluster  frag¬ 
ments  into  one,  according  to  specific  rules  [19  .  These 
fusion  gates,  in  conjunction  with  Bell  states  and  GHZ 
states,  can  be  used  to  construct  arbitrary  cluster  states. 
The  success  probability  of  the  fusion  gates  can  be  im¬ 
proved  to  3/4  if  additional  (ancilla)  single  photons  are 
available  to  be  injected  on-demand  into  an  otherwise- 
passive  linear  optical  circuit,  and  if  the  detectors  have 
up  to  two-photon  number  resolution  m •  We  assume 
such  boosted  fusion  gates  in  our  all-optical  repeater  con¬ 
struction  described  in  this  paper. 

We  model  a  lossy  single  photon  source  of  efficiency 
T]s  as  one  that  emits,  on  demand,  the  mixed  state 
Tjs  1 1)  (1 1  T  (l  —  T]s)  |0)  (0|.  We  use  rjd  for  the  efficiency 
of  all  detectors  in  the  system.  We  will  assume  that 
the  cluster  is  created  on  a  photonic  chip  to  allow  for 
easier  scalability  after  which  the  photons  are  coupled, 
with  efficiency  ?7C,  into  fiber  with  loss  coefficient  a  and 
speed  of  light  Cf.  Pc hip  =  e_^TsCch  denotes  the  sur¬ 
vival  probability  of  a  photon  on-chip  during  one  feed¬ 
forward  step,  where  /3  is  the  loss  coefficient,  cc h  is  the 
speed  of  light  and  rs  is  the  feed-forward  time,  all  on- 
chip.  tjquz  =  Vs'nd/i^  —  VsVd)  is  the  survival  rate  of  the 
photons  that  are  input  into  a  linear-optical  circuit  in¬ 
tended  to  produce  3-photon  maximally-entangled  GHZ 
states  m-  The  final  measurement  step  requires  feed¬ 
forward  in  fiber.  The  survival  probability,  P^b,  during 
feed-forward  time  in  fiber,  t/,  is  P kb  =  e~arfCf.  The 
values  for  device  performance  assumed  for  the  plots  that 
appear  later  in  the  paper,  are  summarized  in  Table  [T| 


III.  REPEATER  ARCHITECTURE 

Before  we  discuss  the  all-photonic  repeater  architec¬ 
ture,  it  is  instructive  to  review  a  generic  quantum  re¬ 
peater  architecture  based  on  multimode  quantum  mem¬ 
ories,  probabilistic  BSMs,  and  multiplexing  over  m  par¬ 
allel  channels  depicted  in  Figs.  [2](a)  and  (b),  which  was 
proposed  by  m,  and  analyzed  in  m-  The  parallel  chan¬ 
nels  can  be  a  combination  of  mutually-orthogonal  spec¬ 
tral,  spatial,  and  polarization  modes,  over  each  of  which 
dual-rail  photonic  qubits  can  be  transmitted  simultane¬ 
ously  at  the  clock  rate  (determined  by  the  source  and 
detector  bandwidth).  Alice  and  Bob  are  separated  by 
optical  fiber  of  length  L  (i.e.,  end-to-end  transmissivity, 
r]  =  e-aL),  interspersed  with  n  repeater  stations  spaced 
Lq  =  L/n  apart,  with  Alice  and  Bob  Lo/2  away  from  the 
terminal  repeaters  in  the  chain. 

Each  of  the  n  repeater  nodes  (or,  ‘major  nodes’), 
shown  by  a  gray  box,  consists  of  a  multimode  quantum 
memory  straddled  between  sources  of  m  Bell  pairs  on 
its  left  and  another  m  on  its  right.  Each  major  node 
loads  one  half  of  an  entangled  Bell  state  onto  the  mem¬ 
ory,  while  transmitting  the  other  half  towards  the  middle 
of  the  adjoining  elementary  link.  Each  major  node  does 


the  above  synchronously  on  every  clock  cycle.  At  the 
center  of  each  elementary  link  is  a  ‘minor  node’,  shown 
as  dark-blue-shaded  boxes  in  Fig.  |2][b).  After  the  qubits 
from  the  major  nodes  reach  the  minor  nodes  (i.e.,  after 
propagation  through  a  distance  Lq/ 2),  each  minor  node, 
simultaneously,  performs  BSMs  on  each  of  the  m  pairs  of 
qubits  received  from  the  repeater  nodes  on  its  either  side. 
The  successful  BSMs  within  each  elementary  link  are 
shown  by  thick  (green)  line  segments.  Immediately  after 
the  minor  node  BSMs,  each  minor  node  sends  back  the 
information — about  which  of  the  m  channels  were  suc¬ 
cessfully  measured — to  its  two  neighboring  major  nodes, 
on  an  authenticated  classical  channel.  Upon  receipt  of 
that  information,  each  major  node  performs  a  BSM  on 
two  qubits  held  in  its  memory  that  had  been  entangled 
halves  of  qubits  that  participated  in  successful  BSMs  at 
the  minor  node  to  the  left  of  that  major  node,  and  the 
minor  node  to  its  right,  respectively.  Simultaneous  with 
the  minor-node  BSMs,  Alice  and  Bob  measure,  in  one  of 
the  two  randomly-chosen  mutually-unbiased  bases,  the 
m  qubits  they  receive  at  their  respective  ends  of  the  ter¬ 
minal  half-elementary-link  segments  (see  Fig.  |2jb)),  and 
send  the  information  about  which  channels  generated  a 
‘click’  on  their  detectors,  back  to  their  respective  neigh¬ 
boring  major  nodes.  Finally,  each  major  node  sends  the 
information  on  whether  its  BSM  succeeded,  to  Alice  and 
Bob.  Hence,  at  every  clock  cycle,  with  some  probability 
(i.e.,  if  all  the  minor  nodes  heralded  at  least  one  suc¬ 
cess  each,  all  major  node  BSMs  were  successful,  and  Al¬ 
ice  and  Bob  both  detected  a  photon  on  at  least  one  of 
the  m  channels  each  while  using  the  same  measurement 
bases),  Alice  and  Bob  obtain  a  shared  (raw,  sifted)  bit. 
A  long  sequence  of  sifted  bits  is  thereafter  used  to  distill 
a  quantum- secure  shared  secret  via  error  correction  and 
privacy  amplification. 

The  all-optical  repeater  architecture  we  now  discuss 
builds  upon  a  recent  proposal  by  Azuma  et  al.  fT2],  al¬ 
though  there  are  some  important  differences,  which  we 
will  point  out  later  in  Section |V|  The  key  idea  is  to  mimic 
a  quantum  memory  (whose  goal  is  essentially  to  protect 
photonic  qubits  against  loss  for  a  certain  time  duration) 
by  using  the  tree  cluster  approach  described  in  Section  [II] 
The  authors  of  [12]  went  one  step  further  and  subsumed 
the  functionalities  of  all  the  subcomponents  of  the  ma¬ 
jor  node  (the  quantum  memory  as  well  as  the  2 m  Bell 
pair  sources)  into  one  single  giant  optical  cluster  state, 
which  we  describe  next.  Fig. [2](c)  illustrates  the  construc¬ 
tion  of  this  cluster.  We  start  with  a  depth-2  star  cluster 
with  a  degree-2m  root  node,  and  4m  +  1  total  qubits. 
The  ‘outer’  qubits,  shown  as  white  circles,  play  a  role 
analogous  to  the  white  qubits  in  Fig.  [2](a)  that  are  trans¬ 
mitted  to  the  minor  nodes  on  fiber  channels.  The  2 m 
‘inner’  qubits,  shown  as  gray  circles,  are  each  attached 
with  a  tree  cluster  of  an  appropriately-chosen  branching 
vector  6,  thereby  creating  a  giant  tree  cluster.  The  loss- 
protected  (logical)  inner  qubits  play  a  dual  role,  that  of 
the  black  qubits  in  Fig.  [2](a)  that  are  held  in  the  quan¬ 
tum  memories  locally  at  the  major  nodes,  and  that  of 
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FIG.  2.  (a)  and  (b)  show  schematics  of  one  elementary  link, 
and  a  chain  of  them  connecting  Alice  and  Bob,  respectively, 
for  a  repeater  architecture  that  employs  quantum  memories, 
Bell  pair  sources,  probabilistic  BSMs,  and  multiplexing  over 
mn  orthogonal  parallel  channels,  (c)  depicts  the  construction 
of  a  photonic  cluster  state  that  can  subsume  the  roles  of  the 
quantum  memory  and  the  Bell  pair  sources,  thereby  resulting 
in  a  quantum  repeater  architecture  based  solely  on  ‘flying’ 
qubits.  The  outer  (white)  photonic  qubits  are  transmitted 
on  the  fiber  channels,  and  the  inner  (black)  qubits  are  held 
locally  in  a  (lossy)  waveguide  at  the  repeater  node.  See  text 
for  a  detailed  description. 


the  memories  themselves.  We  make  the  two  X  measure¬ 
ments  corresponding  to  each  tree  appended  to  the  star, 
as  described  in  the  previous  section  (i.e.,  a  total  of  4m  X 
measurements).  Finally,  we  make  a  Y  measurement  on 
the  root  node  of  the  star,  which  has  an  effect  of  creating 
a  clique  among  all  the  (logical)  inner  qubits,  shown  by 
black  circles  in  Fig.[2|c).  The  clique  of  the  2 m  logical  in¬ 
ner  qubits,  connected  to  the  2 m  outer  qubits,  forms  the 
full  photonic  cluster  state  that  each  major  node  creates 
every  clock  cycle,  and  sends  out  the  2 m  outer  qubits  (the 
white  circles)  towards  the  neighboring  minor  nodes  (m  to 
the  left  and  m  to  the  right)  on  fiber  channels.  Note  that 
the  final  cluster  state  (after  the  X  and  Y  measurements) 
is  not  a  tree. 

Each  major  node  is  equipped  with  single  photon 
sources,  reconfigurable  passive  linear  optics,  and  single 
photon  detectors.  The  clusters  are  created  using  lin¬ 
ear  optics  and  feed-forward  [15,  16  .  Since  the  cluster 
creation  process  is  probabilistic,  the  resources  (number 


of  photon  sources,  detectors,  size  of  linear  optic  circuit) 
must  be  chosen  to  ensure  a  near-unity  success  probability 
of  creating  the  cluster  in  every  clock  cycle  (see  Fig.  [5]). 

The  minor  nodes  are  identical  to  what  was  described 
earlier.  The  remainder  of  the  protocol  proceeds  exactly 
as  described  at  the  beginning  of  this  Section  in  the  con¬ 
text  of  the  memory-based  architecture,  except  for  the 
following  difference  of  the  action  at  the  major  nodes. 
When  the  information  about  which  modes  were  success¬ 
ful  comes  back  at  a  major  node  (from  the  two  neighboring 
minor  nodes),  instead  of  doing  a  BSM  between  a  pair  of 
qubits  held  in  a  memory,  the  major  node  applies  X  mea¬ 
surements  on  the  two  logical  inner  qubits  corresponding 
to  the  successful  modes  on  either  side  of  the  clique,  and 
makes  Z  measurements  on  the  remaining  2 m  —  2  logical 
inner  qubits  (see  Fig.  [2jc)).  The  X  measurements  have 
the  effect  of  fusing  the  successful  outer  qubits  into  an  en¬ 
tangled  chain,  and  the  Z  measurements  have  the  effect 
of  removing  the  extraneous  qubits  from  the  cluster. 

So,  in  any  given  clock  cycle,  if  the  photonic  clusters 
at  each  major  node  are  successfully  created  (which  in¬ 
cludes  success  in  performing  the  4m  X  measurements 
and  one  Y  measurement),  if  all  the  minor  nodes  herald 
at  least  one  BSM  success,  if  the  logical  (inner)  qubits  sur¬ 
vive  the  local  storage  at  the  major  nodes  while  the  outer 
qubits  fly  to  the  minor  nodes  and  the  classical  informa¬ 
tion  (about  which  modes  were  successful)  arrives  back,  if 
the  two  X  measurements  and  2m  —  2  Z  measurements 
done  to  prune  the  clusters  at  the  major  nodes  using  that 
classical  information  are  successful,  and  if  Alice  and  Bob 
get  at  least  one  click  each  while  using  same  measurement 
bases,  then  Alicte  and  Bob  obtain  a  raw  sifted  shared  bit. 
In  Section  |TV|  we  explicitly  calculate  this  overall  success 
probability,  and  the  resulting  secret-key  generation  rate. 
As  we  will  see,  larger  error-protection  trees  afford  better 
rate  performance  (up  to  a  limit  governed  by  the  device 
loss  rates),  but  creating  larger  clusters  at  the  major  nodes 
requires  more  resources  (sources  and  detectors). 

In  Section  [Til  A[  we  describe  in  detail  the  construction 
of  the  clusters  at  the  major  nodes  using  linear  optics, 
and  calculate  the  success  probability.  In  Section  III  B ,  we 


will  describe  how  the  measurements  on  the  major-node 
clusters  are  done,  after  the  BSMs  at  the  minor  nodes, 
to  stitch  together  an  end-to-end  entangled  state  between 
Alice  and  Bob. 


A.  Constructing  the  clusters  at  the  major  nodes 

The  cluster  as  described  above,  prepared  at  each  ma¬ 
jor  node  in  every  clock  cycle,  is  pieced  together  by  fusing 
single  photons  into  progressively  larger  cluster  fragments, 
probabilistically,  using  linear-optical  circuits  and  photon 
detectors.  The  optimal  algorithm  for  creating  photonic 
cluster  states  using  linear  optics — in  terms  of  minimizing 
the  total  number  of  photons  consumed  and  maximizing 
the  eventual  probability  of  success — is  not  known  even  for 
a  general  A-node  line  cluster.  With  losses  from  sources 
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detectors  and  waveguides  during  cluster  construction, 
finding  the  optimal  recipe  becomes  even  harder.  One  de¬ 
sign  knob  is  the  number  of  redundant  cluster  fragments 
attempted  at  each  step.  A  higher  number  of  attempts 
improves  the  probability  of  successfully  creating  the  fi¬ 
nal  cluster,  but  with  a  higher  number  of  required  photon 
sources  and  detectors.  We  refer  to  this  trick  of  attempt¬ 
ing  the  creation  of  multiple  identical  cluster  fragments  at 
each  step  of  the  process  as  multiplexing. 

We  now  describe  the  resource  counts  and  success- 
probability  calculations  for  two  methods  to  create  the 
cluster  at  the  major  node.  The  first  one  is  a  method 
implied  by  previous  rough  estimates  of  the  resource  re¬ 
quirements  pans].  We  then  discuss  an  improved  scheme 
that  decreases  the  resource  requirements  during  the  cre¬ 
ation  process.  Fig.  [4]  provides  a  schematic  for  these  two 
schemes,  which  we  refer  to  in  the  discussion  below. 


Number  of  potentially- 
distinct  clusters  at 
each  stage 


(2°) 


(21) 


C 


km 


t  (4m  X  measurements  and 
I  one  Y  measurement) 

ck 


ct 1 


•Sf  If 


(2fc) 


(GHZ  state) 


ro 

(GHZ  state) 


(GHZ  state) 


FIG.  3.  The  tree  cluster  Ck  (and  the  final  cluster  Ckrn  after 
the  X  and  Y  measurements),  shown  in  Fig. [2j  are  created  by 
a  sequence  of  probabilistic  linear-optical  fusion-II  operations, 
starting  from  3-photon  maximally-entangled  (GHZ)  states. 


Let  us  label  the  final  cluster  Ckrn  (see  Fig.  [2])  where 
the  letter  m  signifies  that  the  Y  measurement  required 
to  turn  the  inner  qubits  of  the  star  into  a  clique  (a  fully 
interconnected  graph)  and  the  X  measurements  required 
to  connect  the  error  protection  trees  to  the  inner  qubits 
have  already  been  applied.  Before  these  measurements, 
the  (tree)  cluster  is  labelled  as  Ck .  We  label  the  daughter 
clusters  that  are  fused  together  to  create  Ck  as  Ck _1  and 
Ck~x .  The  daughter  clusters  that  are  fused  together  to 
create  C\~x  are:  CkJ2  and  Ck^2.  The  clusters  that  are 
fused  together  to  create  Ck 22  are:  Ck 2\  and  Cf  2  2,  and 
so  on  (See  Fig.  [3|.  At  the  bottom  of  the  stack  are  3- 
photon  GHZ  states,  (7?  with  i  =  hd2  •  •  •  G/c,  which  are 
in  turn  created  by  groups  of  6  photons  fed  into  linear- 
optical  circuits  that  generate  the  3-photon  GHZ  states 
with  probability  Pghz  =  [^^(2  -  r]sVd)f  /32  [15] .  The 
loss  rate  of  the  heralded  GHZ  states  is,  1  —  ??ghz  where 

Vguz  =  VsVd/iZ  ~  VsVd )  US- 

We  assume  that  the  cluster  Ck  can  be  prepared  in 
a  series  of  k  fusion  steps,  where  at  each  step,  clusters 
of  roughly  equal  sizes  are  fused  together,  thus  roughly 
doubling  the  cluster  size  in  each  step  m-  This  as¬ 


sumption  becomes  accurate  in  the  limit  of  large  clusters. 
This  method  ties  the  final  size  of  the  intended  cluster 
(TVciuster  =  2k  +  2  photons)  to  the  number  of  fusion  steps 
(fc),  and  this  relationship  becomes  increasingly  exact  as  k 
becomes  large.  In  other  words,  we  assume  that  Cl^\  and 
C \~2  are  two  clusters  each  of  p  photons,  which  when  fused 
successfully  using  a  fusion-II  gate  (applied  to  one  photon 
each  of  the  above  two  clusters)  creates  the  2p  —  2  photon 
cluster  C],  i  =  i\,  12  . . . ,  ik-i-  Starting  with  the  3-photon 
GHZ  states  C^i  i<2  ikl  the  size  of  Ck  is  2k  +  2  photons. 
Hence,  the  minimum  number  of  fusion  steps  required  to 
build  a  A^duster  photon  cluster  is  k  =  [log2(7Vciuster  -  2)]. 
The  label  fc,  the  number  of  fusion-II  steps  used  to  arrive 
at  Ck:  also  translates  to  the  resource  requirements,  and 
the  loss  rate  of  each  photon  in  the  final  cluster,  as  we 
show  below.  Note  that  k  is  a  function  of  the  branching 
vector  b  of  the  error-correction  trees  used.  The  larger 
the  error-correction  trees,  the  larger  is  the  final  cluster 
Ck,  and  the  larger  is  the  number  of  steps  k  required  to 
prepare  that  cluster. 


1.  The  naive  multiplexing  scheme 


Let  us  now  examine  the  cluster  creation  process  (de¬ 
picted  for  k  =  2  in  Fig.  [4ja)).  At  every  point  we  need 
the  cluster  fragment  C],  we  attempt  to  create  ub  copies 
of  that  identical  cluster  (rts  =  3  shown  in  Fig.  Eh)), 
of  which  hopefully  one  is  successfully  created  and  her¬ 
alded  for  further  use.  Therefore,  creating  one  usable  copy 
of  Ck  requires  (2? %s)k  GHZ  states  C^i  i2  ik  at  the  bot¬ 
tom  of  the  stack.  Each  GHZ  state  is  picked  from  ughz 
parallel- attempted  GHZ  states  (tighz  =  4  shown  in  Fig. 
|4ja)),  and  creating  each  GHZ  state  requires  6  single  pho¬ 
tons.  Therefore,  creating  one  usable  copy  of  Ck  requires 
(2nB)k  x  6riG hz  single  photons.  Finally,  at  the  top  of  the 
chain,  we  create  nmeas  copies  of  Ck  in  parallel  (nmeas  =  4 
shown),  on  each  of  which  the  4m  X  measurements  and 
one  Y  measurement  are  performed,  to  prepare  copies  of 
the  final  required  cluster  Ckrn.  We  choose  nmeas  such 
that  we  obtain  with  high  probability  one  successfully- 
created  copy  of  Ckm.  Therefore,  the  total  number  of  sin¬ 
gle  photon  sources  (shown  by  black  dots  at  the  bottom  of 
Fig.Qa))  that  need  to  simultaneously  fire  on  every  clock 
cycle,  Ns  =  6nGHZ  nmeas(2nB)fc. 

The  probability  of  successfully  creating  a  GHZ  state 
is  Po  =  1  -  (1  -  Pg Hz)nGHZ-  The  success  prob¬ 
ability  of  fusion  at  the  Z-th  step — i.e.,  that  of  combining 
CX  and  77  into  Gi  is  given  by  Qt  =  (rjG}izP^hip)2 /2. 
The  success  probability  of  heralding  one  cluster  C\  (from 
the  ub  parallel  copies  attempted)  is  given  by  the  recur¬ 
sive  formula,  Pi  =  1  —  (1  —  P2_1Qi)TlB ^  with  P0  given 
as  above.  The  4m  X  measurements  and  one  Y  mea¬ 
surement  required  to  convert  Ck  to  the  final  cluster  Ckm 

/  \4m+l 

succeed  with  probability  P'  —  IvghzPck[p  )  .  Since 

this  step  is  multiplexed  over  nmeas  parallel  attempts,  the 
success  probability  of  heralding  one  copy  of  the  final  clus- 


ter  at  a  major  node  is  given  by,  Pc\  =  1  —  (QkP')nmea-s  - 
The  success  probability  of  all  n  repeater  nodes  creating 
the  clusters  Ckm  locally  during  any  given  clock  period, 
is  Pcn  =  P^.  The  blue  (dashed)  plot  in  Fig. [5] shows  Pcn 
as  a  function  of  Ns  for  n  =  250  repeater  stations  (ma¬ 
jor  nodes),  k  =  7,  and  for  device  parameters  as  given  in 
Tabled 

2.  The  improved  multiplexing  scheme 

The  improved  multiplexing  scheme  we  now  describe 
addresses  the  following  deficiencies  of  the  scheme  de¬ 
scribed  above. 

•  The  protocol  presented  above  does  not  make  the 
most  optimal  use  of  the  multiple  copies  of  identi¬ 
cal  clusters  that  are  successfully  created  at  a  given 
step.  To  illustrate  this  point,  let  us  consider  the 
tib  =  3  copies  of  (attempted)  C 2  clusters  that 
are  shown  in  Fig.  [4ja),  of  which  one  successfully 
created  C 2  is  picked.  The  first  of  those  three  at¬ 
tempted  C 2  clusters  is  shown  to  be  created  by  fus¬ 
ing  a  C\  cluster  and  a  C\  cluster.  The  C\  is  cho¬ 
sen  out  of  tib  =3  copies  of  (attempted)  C\  clus¬ 
ters,  as  shown.  If  two  of  those  three  copies  of  C\ 
are  actually  successfully  created,  the  second  suc¬ 
cess  goes  waste.  Note  however  that  the  second  and 
the  third  (of  the  three  attempted)  C 2  clusters  also 
each  need  to  be  created  by  fusing  a  C\  and  a  C\. 
Those  two  C\  clusters  are  also  picked  from  ub  —  3 
copies  each  of  (attempted)  C\  clusters  (not  shown 
in  the  figure).  It  is  thus  simple  to  see  that  at  each 
time  step,  a  total  of  ( ubY  =  9  copies  of  C\  are  at¬ 
tempted,  but  the  selection  of  successes  only  happen 
within  groups  of  three,  which  is  clearly  inefficient. 
A  far  more  efficient  approach  is  to  maintain  one 
single  “bank”  of  copies  of  C\  and  similarly  one  sin¬ 
gle  bank  for  copies  of  C\ ,  and  attempt  fusions  on 
clusters  from  these  two  banks  pairwise  (and  throw 
away  the  excess  clusters  in  the  bank  that  has  more 
copies),  to  produce  a  single  bank  of  C 2  clusters. 
This  way,  one  does  not  have  to  choose  the  multi¬ 
plexing  numbers  ncHZ  and  nmeas,  and  the  total 
number  of  single  photons  Ns  directly  translates  to 
an  overall  probability  of  success  Pc\  of  creating  the 
final  cluster  Ckrn.  In  general,  we  maintain  single 
banks  of  each  distinct  cluster  fragment  consumed 
in  the  entire  stack  shown  in  Fig.  [3J  and  for  each  fu¬ 
sion  step  shown  in  Fig.  [3j  we  apply  pairwise  fusion 
to  all  cluster  copies  from  the  two  banks  correspond¬ 
ing  to  the  two  daughter  clusters  (and  throw  away 
the  excess  clusters  from  the  bank  that  has  more). 

•  The  X  and  Y  measurements  that  were  performed 
at  the  very  end  (on  4m  + 1  nodes  of  the  tree  cluster 
Ck ,  to  convert  it  to  the  required  final  cluster  Ckm ) 
can  be  performed  at  the  very  beginning — on  the  ap¬ 
propriate  photons  (which  would  eventually  become 


those  4m  +  1  photons  in  Ck) — while  they  are  still 
part  of  the  3-photon  GHZ  states,  i.e.,  before  any  of 
the  fusion-II  operations  begin.  Making  these  mea¬ 
surements  at  the  bottom  of  the  stack  makes  failures 
much  less  costly,  which  in  turn  significantly  reduces 
the  resource  requirements  (i.e.,  the  Ns  required  to 
achieve  a  given  final  success  probability  Pcn).  Ap¬ 
pendix  [A]  rigorously  explains  why  these  measure¬ 
ments  can  be  done  on  the  photons  while  they  are 
still  parts  of  the  GHZ  states. 

•  The  success  probability  of  each  of  the  fusion-II  op¬ 
erations  (at  all  k  steps  in  the  cluster  creation  pro¬ 
cess)  can  be  improved  from  1/2  to  3/4  by  injecting 
ancilla  single  photons  im.  These  success  proba¬ 
bility  numbers  diminish  with  source  and  detection 
inefficiencies.  But,  the  cost  of  using  additional  pho¬ 
tons  needed  (as  ancillas)  to  realize  these  boosted 
fusion  gates  is  far  outweighed  by  the  effect  of  the 
success-probability  improvement,  thereby  improv¬ 
ing  the  effective  tradeoff  between  Ns  and  Pcn. 

We  start  with  Ns  photons  and  send  them  all  through 
GHZ  factories,  hence  attempting  the  creation  of  [Ns/6\ 
3-photon  GHZ  states.  The  number  of  GHZ  states 
x  successfully  created  follows  a  binomial  distribution 
B(x,  |AV6J,Pghz)  where  B(x,n,p)  =  (”)px(  1  -p)n~x. 
Hereonafter,  let  us  follow  an  illustrative  set  of  numbers 
for  a  k  =  2  cluster,  which  is  depicted  schematically  in 
Fig.  gb).  Suppose  we  get  x  =  18  successfully-created 
GHZ  states.  These  GHZ  states  are  now  split  into  4  banks 
corresponding  to  C®  1?  ^21  and  ^22-  Out  of  these, 

let  us  say  C®  1  and  C °  2  consist  of  photons  that  would  be 
eventually  measured  in  Ck .  As  discussed  in  Appendix  |a| 
these  qubits  can  be  measured  now.  Since  the  measure¬ 
ment  of  photons  has  a  success  probability  Pchip^GHZ,  the 
number  of  C®™  cluster  states  (x)  created  as  a  result  of 
making  measurements  on  y  C®  1  states  follows  a  binomial 
distribution  P(x,  y,  PchiP^GHz)-  The  banks  correspond¬ 
ing  to  Ci  1  and  C22  are  given  a  fraction  l/(Pchip^GHz) 
more  GHZ  states.  Hence,  these  banks  have  5  GHZ  states 
each  whereas  the  other  two  have  4  each.  Suppose  that 
measuring  the  5  copies  of  Ci  1  results  in  4  copies  of  C®™, 
and  measuring  the  5  copies  of  C °  2  results  in  4  copies  of 
C2 *2  .  The  first  fusion  step  is  now  attempted  (i.e.,  fusing 
Ci™  with  Ci  2^  and  fusing  C21  with  C2™)  resulting  in  2 
successfully  created  copies  of  C\ m  and  3  copies  of  C 2m 
(the  maximum  possible  number  of  successes  in  both  cases 
was  4).  In  the  final  step,  there  are  2  fusion  attempts  from 
which  we  get  one  copy  of  the  final  cluster  state  C2rn. 

In  general,  in  a  level-/  fusion  step  in  Fig.  [3j  and 
with  yi  and  y2  copies  in  the  respective  banks  of  the 
two  daughter  clusters,  the  distribution  of  the  number  x 
of  fused  states  C\  is,  B  (x,  min{?/i,  t/2},.Pz),  where  pi  = 
A  {li.VsVd)2  +  \ (VsVd)4)  HZ]  and  m  =  rteHZ^hip  is  the 
survival  rate  of  photons  up  to  before  the  Ith  fusion  step. 
The  success  probabilities  of  this  scheme,  Pc\  (and  Pcn) 
are  calculated  using  Monte  Carlo  simulations. 
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FIG.  4.  (a)  the  naive  multiplexing  scheme.  A  dashed  rectangle  represents  a  cluster  that  has  some  probability  of  having  been 
been  created  after  a  probabilistic  fusion  step  (red  circle)  or  at  the  output  of  creating  GHZ  states  using  linear  optics  starting 
from  six  single  photons  (labeled  ‘GHZ  Factory’).  A  solid  rectangle  represents  a  cluster  state  that  is  successfully  created  with 
high  probability  by  choosing  a  successful  outcome  (blue  square)  out  of  several  identical  copies  attempted  (dashed  boxes),  (b) 
the  improved  multiplexing  scheme.  A  box  surrounding  clusters  of  the  same  type  represents  a  bank  of  clusters  and  any  operation 
applied  to  the  bank  is  applied  to  all  the  clusters  in  it. 


In  Fig.  [5j  we  plot  the  probability  Pcn  of  successfully 
building  clusters  Ck  (with  k  =  7),  simultaneously  at 
n  =  250  major  nodes,  for  both  schemes,  ughz  and 
nmeas  are  optimized  for  the  naive  scheme  to  maximize 
Pcn  for  any  given  Ns.  The  plot  clearly  shows  that  the 
improved  scheme  leads  to  resource  savings  by  a  factor  of 
^  104.  We  further  observe  that,  for  both  schemes,  Pcn 
undergoes  a  rapid  percolation-like  transition  from  zero  to 
one  as  Ns  is  increased  beyond  a  certain  threshold  value. 
Pcn  is  only  a  function  of  fc,  n,  and  Ns.  We  fix  Pcn  =  0.9 
and  calculate  the  corresponding  minimum  Ns  required, 
for  every  value  of  k  and  n.  This  sharp-transition  behav¬ 
ior  of  Pcn  allows  us  to  conveniently  split  the  problem  of 
designing  the  repeater  architecture  into  two  parts: 

(1)  choosing  an  error-protection  level  by  choosing  m 
(number  of  parallel  qubit  channels)  and  b  (the  branching 


vector  of  the  error  protection  trees),  which  gives  us  k 
(indicative  of  the  total  cluster  size),  and  using  this  to 
calculate  the  key  rate  vs.  distance  achieved — both  with 
n  repeater  stations,  and  also  the  resulting  envelope  over 
all  n;  and 


(2)  given  the  design  choices  (m  and  5),  calculating  the 
number  of  photon  sources  Ns  so  as  to  achieve  a  close- 
to-unity  Pcn  (probability  that  all  n  nodes  create  the  re¬ 
quired  clusters  on  every  clock  cycle),  for  a  given  value  of 
k  (cluster  size  at  each  repeater  node),  and  n  (the  number 
of  repeater  nodes). 
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FIG.  5.  The  probability  that  all  n  —  250  major  nodes  are 
simultaneously  successful  in  creating  clusters  of  size  k  —  7 
fusion  steps  (i.e.,  2k  +  2  =  130  photon  clusters),  using  the 
naive  and  the  improved  multiplexing  schemes. 


B.  Measuring  the  clusters  and  connecting  the  chain 

Once  the  clusters  are  created,  the  outer  qubits  are  sent 
to  minor  nodes  at  the  middle  of  the  elementary  links,  as 
shown  by  the  arrows  in  Fig.  [2jc).  The  outer  qubits 
are  measured  in  the  Bell  basis  at  the  minor  nodes  using 
ancilla-assisted  boosted  fusion  gates  PH-  The  loss  rate 
seen  by  the  outer  qubits  is  etrav  =  1  —  V2nPd tepVGHZVc 
where  77 1/2n  is  the  transmissivity  of  half  of  an  elemen¬ 
tary  link  (of  range  L/2n).  All  the  physical  qubits  cor¬ 
responding  to  the  inner  (logical)  qubits  are  stored  lo¬ 
cally  in  a  fiber  bundle  with  the  same  attenuation  as  the 
communication  fiber  between  the  repeater  stations.  Due 
to  the  classical-communication  delay,  the  core  qubits  see 
more  loss  than  the  outer  qubits  do,  which  we  define  as 
fistat  =  1  -  ??"-Pchtp2-Pfib^GHZ??c-  However,  it  is  important 
to  note  that,  just  like  in  the  architecture  of  mm,  this 
delay  only  leads  to  a  latency  in  the  scheme  and  does  not 
affect  the  clock  rate  of  the  system. 

When  the  result  of  the  BSMs  on  the  m  qubit  channels 
at  the  two  neighboring  minor  nodes  arrive  back  at  a  ma¬ 
jor  node,  the  major  node  picks  one  successful  qubit  chan¬ 
nel  on  either  side  (if  none  of  the  m  BSMs  were  a  success 
on  any  one  of  the  sides,  then  that  time  period  is  an  overall 
failure) .  The  logical  inner  qubits  corresponding  to  all  the 
outer  qubits  that  are  not  deemed  part  of  the  successful 
BSMs  are  removed  from  the  cluster  by  measuring  them 
in  the  Z  basis  [T4]  (note  that  this  Z  measurement  is  a  log¬ 
ical  one,  which  benefits  from  the  loss-protection  trees). 
On  the  two  logical  qubits  (one  on  either  side)  correspond¬ 
ing  to  the  successful  channels,  X  basis  measurements  are 
performed,  which  has  an  effect  of  extending  the  entangle¬ 
ment.  Alice  and  Bob,  simultaneous  with  the  minor  node 
BSMs,  detect  the  m  outer  photons  sent  to  them  by  the 
first  and  the  last  major  node  in  the  repeater  chain,  over 
links  of  length  Tq/2,  using  one  of  two  randomly-chosen 


mutually- unbiased  bases.  Assuming  the  clusters  at  all  n 
repeater  nodes  were  successfully  created  (which  happens 
with  probability  Pcn),  the  conditional  probability  of  gen¬ 
erating  an  end-to-end  entangled  pair  between  Alice  and 
Bob,  in  one  clock  cycle,  is  given  by  the  probability  that 
all  n—  1  minor  nodes  herald  at  least  one  successful  BSM, 
and  all  the  pruning  logical  X  and  Z  measurements  on 
the  clusters  at  all  n  major  nodes  are  successful,  and  Al¬ 
ice  and  Bob  both  obtain  successful  detects  on  at  least 
one  of  the  m  qubit  channels: 

Pmeas  =  p2n  [1  _  (i  _  P^,,  (1) 

where  Px  and  Pz  are  the  probabilities  of  successful  X 
and  Z  basis  measurements  on  the  logical  inner  qubits, 
respectively.  Pend  is  the  probability  that  Alice  (resp., 
Bob)  obtains  at  least  one  successful  detection  in  one  of 
the  m  qubit  channels. 

We  quantify  the  performance  of  the  repeater  architec¬ 
ture  in  terms  of  the  number  of  shared  secret  bits  gener¬ 
ated  per  mode  (i.e.,  per  clock  cycle  per  spatial  channel, 
where  m  is  the  number  of  spatial  channels  employed). 
Since,  the  channel  noise  comprises  of  only  photon  loss, 
the  success  probability  divided  by  the  number  of  spatial 
channels  per  attempt  is  the  secret  key  rate  (in  bits  per 
mode)  generated  by  this  scheme,  i.e.,  R  =  PcnPmeas/2m 
bits/mode.  Note  that  the  bits  per  mode  is  obtained  by 
dividing  by  the  number  of  spatial  channels  that  is  twice 
the  number  of  qubit  channels  (2m).  This  is  because  we 
assume  single-polarization  dual-rail  encoding  where  each 
qubit  on  any  given  spatial  channel  occupies  two  succes¬ 
sive  temporal  modes. 


IV.  RATE  CALCULATIONS 

In  this  Section,  we  evaluate  the  secret  key  rate  achiev¬ 
able  using  the  all-optical  repeater  architecture  described 
above,  while  accounting  for  all  the  device  and  channel 

losses.  We  first  evaluate  an  expression  for  Pnm’^(P),  the 
bits-per-mode  rate  for  a  given  choice  of  design  parame¬ 
ters:  m  (the  number  of  parallel  channels)  and  b  (branch¬ 
ing  vector  of  the  error-protection  trees).  L  is  the  Alice- 
to-Bob  range  and  n  is  the  number  of  equally-spaced  re¬ 
peater  nodes  that  are  deployed  between  Alice  and  Bob. 
We  evaluate  the  rate- vs. -distance  envelope  R^m,b\L) — 
the  maximum  of  R^^^L)  at  any  L  over  the  choice  of 
n  G  {1,2,...} — and  we  show  explicitly  for  when  b  is 
a  depth-2  tree,  that  P^m,^(P)  >  Dr]s ,  with  D  a  con¬ 
stant,  r]  =  e~aL  and  s  strictly  less  than  1.  We  find 
by  numerical  evaluation  that  this  lower  bound  is  tight. 
We  compare  this  rate-distance  envelope  with  the  best 
rate  achievable  without  the  use  of  quantum  repeaters 
^direct  (L)  =  —  log2(l  —  77),  for  some  (m,  b)  pairs. 

A  given  choice  of  m  and  b  determines  fc,  the  num¬ 
ber  of  fusion  steps  required  to  prepare  the  final  cluster 
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Ck  prepared  by  each  repeater  node  at  every  clock  cy¬ 
cle,  which  in  turn  quantifies  the  size  (iVciuster  =  2^  +  2 
photons)  of  Ck.  Next,  we  choose  a  value  of  k — a  single 
parameter  that  quantifies  the  amount  of  resources  we  are 
willing  to  dedicate  to  each  repeater  node — ,  and  numer¬ 
ically  optimize  the  choice  of  m  and  b  that  is  consistent 
with  the  chosen  and  which  maximizes  the  rate.  We 
denote  the  rate  attainable  with  n  repeater  nodes  con¬ 
ditioned  on  the  per-node-resource-constraint  parameter 

as  Rn\L)  and  calculate  the  optimal  rate- vs. -distance 
envelope  R^k\L)  by  taking  an  envelope  over  the  choice 
of  n.  Finally,  we  compare  the  rate-distance  envelopes 
for  increasing  values  of  k  and  translate  the  values  of  k 
to  the  number  of  single  photon  sources  required  at  each 
repeater  node. 

The  probabilities  of  fault-tolerant  X  and  Z  measure¬ 
ments  on  one  of  the  (logical)  inner  qubits  of  a  major 
node  cluster,  Px  and  Pz,  can  be  expressed  in  terms  of 
the  probabilities  ^  of  a  successful  ‘indirect’  Z  measure¬ 
ment  (as  described  in  Section  0  on  a  qubit  at  the  i-th 
level  of  the  error-protection  tree  imm- 


where  C  =  P^VghzVc- 

We  now  have  the  bits-per-mode  rate  achievable  with 
an  n-repeater-node  chain, 

R^\L)  =  ^p2ndp2(— l)«P2n  [1  _  (1  _  PbH™-1  , 

(9) 

with  Px,  Pz,  Pb  and  Pen d  as  given  in  Eqs.  0,  0,  0 
and  0,  with  r]  —  e  aL  the  transmissivity  of  the  end-of- 
end  channel  (of  range  L).  See  the  dotted  magenta  curves 

in  Fig.  ^  for  the  plots  of  Pim’^(P)  as  a  function  of  L  for 
a  few  chosen  values  of  n. 

One  way  to  obtain  a  lower  bound  of  the  envelope  over 
the  plots  Pim,^(P)  over  all  choices  of  n  (see  black  plot 
in  Fig.  |6]),  is  to  pick  one  point  (Ln,  Pim,^(Pn))  on  each 

of  the  rate-distance  functions  p!m’^(P),  n  =  0, 1,  2, . . ., 
and  connect  them.  Let  us  choose  Ln  as: 

Ln  =  nz  In (AB2)/a,  (10) 


Px  =  £o,  and 

Pz  ~  (1  estat  T  ^stat^l)  °5 


with  z  being  a  constant  that  is  yet  to  be  chosen.  The 
(2)  Alice-to-Bob  channel  transmissivity  at  these  range  values 
{£\  are  therefore  given  by: 


where, 


r]n  =  e~aLn  =  e~nzln(AB2)  m 


(11) 


&  =  1-  [1  -  (1  -  estat)(l  -  Cetat  +  estat&+2)6*+1]**  ,  (4) 

and  i  <  /,  f/+i  =  0,  bi+ j  =  0. 

Let  us  assume  a  tree  depth  of  d  =  2,  i.e.,  b  =  [bo  &i], 
which  is  consistent  with  our  numerical  findings  on  the 
optimal  branching  vector  as  described  later  in  the  pa¬ 
per  (see  table  0-  For  a  depth-2  branching  vector,  using 

Eq.  0,  we  find  that  £0  =  1  -  [l  -  (l  -  estat)6l+1)]  b°  and 

£l  =  1  —  4taf  Thus> 


We  now  evaluate  a  locus  of  the  (range,  rate)  pairs 

(Ln,  R(hn,b\Ln))  over  n  E  {0,1,2,...}  and  choose  the 
parameter  z  we  left  undetermined  in  Eq.  ©  so  as  to 
maximize  the  rate-distance  envelope.  We  call  this  rate- 

distance  envelope  R^b\L)  since  this  is  by  construction 
a  lower  bound  on  the  true  envelope  p(m’6)(P). 

Let  us  evaluate  Px,  Pz,  Pb  and  Pen d  at  r]  =  r]n  (i.e., 
substitute  r\xln  =  (AB2)  2  in  the  respective  expressions) 
and  define  the  following  quantities: 


Px  =  1- 


1  - 


(^)6l+V+*f 


Pz  = 


1  —  ^1  —  Tjn  B^j 


bi  +  l' 


-i  b0 


and  (5) 

(6) 


and  the  Bell  measurement  success  probability  becomes 


„  AB2  i 
Pb  = - 


(7) 


where  A  =  m(±(i]si]d)2  +  ife^d)4)/P|b,  B  = 


PcMv  PfibVGYLZVc- 


hip  iiu'/vjrazj'/c 

The  probability  of  at  least  one  successful  detection  at 
Alice’s  (or  Bob’s)  end  is  given  by 


Pend  — 


/  i  \  rn 

=  l-(l , 


(8) 


Px  =  1  - 


1  -  (AB2) 


2\~z(b  i  +  l) 


Bbl+1 


ib0 


PZ  = 


b  i  +  l' 


1  -  (l  -  (AB2)  z  #) 

Pb  =  —  (^4-B2)1  ^  ,  and 
rn  ' 

Pend  =  l-(l  ~  (AB2)~Z/2  a)"  , 


-i  b0 


and  (12) 

(13) 

(14) 

(15) 


using  which  let  us  define  the  following:  q\  —  p ^ p2x , 
q2  =  1  -  (1  -  Pb )m,  and  q3  =  p2nd,  and  obtain: 


Rirrb\Ln)  =  (qiq2r  (16) 


To  obtain  the  envelope  Rj^b\L),  we  need  to  calcu¬ 
late  the  locus  of  the  distance-rate  pairs  (Ln,  Plm,6^(Pn)) 
over  n  E  {1,2,...}.  We  do  this  by  eliminating  n  from 
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Eqs.  ([10])  and  ([16]).  With  a  little  algebra,  and  expressing 
the  envelope  in  terms  of  77  =  e_aL,  we  get  the  following: 


4b  %)  =  DVS, 


(17) 


where  D  =  and  the  exponent  s  =  — 


Note  that  R^b\L)  in  (fl7|)  is  a  lower  bound  on  the 

actual  rate-distance  function  p(m’6)(P)  for  any  value  of 
the  parameter  2  that  we  left  undetermined  in  our  choice 

of  the  range  values  Ln  we  used  to  evaluate  ( L ).  We 

numerically  optimize  the  choice  of  z  such  that  the  value 
of  the  exponent  s  is  minimized  (note  that  <4,  and  q% 
are  all  functions  of  z). 

In  Fig.pl  we  plot  R^'b\L)  (bits  per  mode)  as  a  func¬ 
tion  of  L  (km)  for  n  =  1, 10,  24,  56, 133,  and  314  (magenta 


dotted  plots),  with  b  =  {7,  3}  and  m  =  4,  and  other  de¬ 
vice  parameters  as  summarized  in  Table  [I]  These  values 
of  m  and  b  translate  to  k  =  8,  i.e.,  28  +  2  =  258  pho¬ 
ton  clusters  created  at  each  node  at  every  clock  cycle. 
We  also  plot  the  analytical  rate-envelope  lower  bound  in 

Eq.  (|T7|) ,  R^h\L)  (black  solid  line),  with  the  optimal 
z  computed  numerically.  For  the  chosen  parameters,  we 
get  D  =  0.11  and  s  =  0.37.  The  analytical  lower  bound 

R^b\L)  is  visually  indistinguishable  at  the  scale  of  the 
plot  from  the  numerically-obtained  rate-distance  enve¬ 
lope  R(m,b\L).  This  excellent  agreement  persists  for  all 
values  of  m  and  b  we  have  have  tried. 

One  interesting  implication  of  the  range  values  Ln  in 
Eq.  (10)  lying  on  the  rate-distance  envelope  is  that  the 
distance  between  each  repeater  (major)  node, 


ln(AB2) 

a 


(18) 


is  a  constant  and  independent  of  the  total  range  L.  In 
other  words,  given  the  device  parameters  and  the  choice 
of  the  major-node  cluster  size  (i.e.,  m  and  5),  there  is  an 
optimal  gap  with  which  repeaters  should  be  placed — no 
more,  and  no  less.  For  the  numbers  used  for  the  plots 
in  Fig.  [6|  Lo  =  1.49  km.  Fig.  [6]  also  shows  Pdirect(A) 
for  comparison  (blue  dashed  plot),  which  the  repeater 
scheme  is  seen  to  outperform  beyond  a  range  of  87  km. 

As  shown  by  the  above  example,  our  repeater  scheme, 
even  when  built  with  lossy  components,  can  achieve 
s  <  1  i.e.  it  outperforms  the  optimum  repeater-less 
rate  indirect  (A).  The  value  of  the  exponent  s  achiev¬ 
able  by  the  repeater  scheme  can  be  improved  (lowered) 
by  enhancing  the  level  of  error  correction  (i.e.,  choos¬ 
ing  a  larger  0).  Doing  so  increases  the  size  of  the  clus¬ 
ters  ( 2k  +  2  photons)  needed  at  each  repeater  nodes, 
and  hence  increases  the  number  of  photon  sources  Ns 
required  locally  at  each  node.  In  Fig.  [7]  we  plot  the 
R^k\L),  numerically-evaluated  envelopes  of  the  rate- 
distance  functions  Ri?\L),  parametrized  by  the  single 
parameter  k  that  quantifies  the  size  of  the  clusters  pre¬ 
pared  by  the  repeaters  at  each  clock  cycle.  It  is  seen  that 


the  rate-distance  exponent  s  improves  (decreases)  as  k 
increases. 


<D 

O 

B 

QJ 
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FIG.  6.  The  key  rate  (in  bits  per  mode)  R ( L )  achieved  by 
an  n-node  repeater  chain  shown  as  a  function  of  range  L,  for 
n  s=  1,10,24,56,133,  and  314  (magenta  dotted  plots),  with 
m  —  4  parallel  channels  and  b  =  {7,  3}  trees.  The  analytical 

lower  bound  to  the  rate-distance  envelope  R^,b\L)  (black 
solid  plot)  is  seen  to  surpass  the  best-possible  repeaterless- 
QKD  rate  indirect  (L)  (blue  dashed  plot)  at  L  =  87  km. 


Device  parameter 

symbol 

value 

fiber  loss  coefficient 

a 

0.046  km"1 
(0.2  dB/km) 

on-chip  loss  coefficient 

p 

0.62  m_I  (2.7 
dB/m) 

feed-forward  time  in  fiber 

Tf 

102.85  ns 

feed-forward  time  on-chip 

Ts 

20  ps 

chip  to  fiber  coupling  efficiency 

Vc 

0.99 

source  detector  efficiency  product 

VsVd 

0.99 

speed  of  light  in  fiber 

cf 

2  x  10 *m/s 

speed  of  light  on  chip 

Cch 

7.6  x  107m/s 

TABLE  I.  Assumed  values  for  device  performance  parameters. 
The  source  detector  efficiency  product  r/srid  is  sufficient  for 
the  purposes  of  the  calculations  in  this  paper,  and  need  not 
be  specified  separately.  Recall  that  Pchip  =  e_/3rsCch,  Pfib  = 
e~ OLT f ,  anc[  rjGuz  —  rjsrjd/(2  —  VsVd)-  Tf  has  been  chosen  to 
make  Pchip  —  Pfib* 


V.  DISCUSSION 

In  this  Section,  we  go  back  to  the  all-photonic  repeater 
architecture  proposed  by  Azuma  et  al.  m,  and  discuss 
the  main  modifications  (improvements)  we  considered  in 
the  architecture  we  described  and  analyzed  above.  We 
also  show  a  comparative  study  of  the  resource  require¬ 
ments  and  rate  performance  of  the  naive  scheme  and  our 
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modified  scheme.  Following  are  the  salient  differences 
between  the  architecture  we  analyzed  above,  and  the  one 
proposed  in  12  . 

Retaining  vs.  transmitting  the  clusters — In  the  pro¬ 
posal  of  [12],  all  the  logical  inner  qubits,  along  with  the 
outer  qubits  (i.e. ,  all  the  N  photons  of  the  cluster  at  a 
major  node)  are  sent  to  the  minor  node,  whereas  we  store 
the  inner  qubit  photons  in  a  fiber  spool  locally  at  the  ma¬ 
jor  nodes.  The  former  has  an  advantage  that  no  classical 
communication  needs  to  happen  from  minor  nodes  back 
to  major  nodes  before  the  logical  X  and  logical  Z  mea¬ 
surements  are  done  to  the  logical  inner  qubits,  since  all 
those  qubits  are  present  locally  at  the  minor  nodes  when 
the  BSMs  are  performed  there  on  outer-qubit  pairs  from 
neighboring  major  node  clusters.  The  advantage  of  our 
(latter)  scheme  is  that  the  number  of  parallel  physical 
channels  needed  (2m)  is  much  smaller  as  compared  to 
the  number  needed  (TV)  for  the  scheme  in  [12].  For  the 
numbers  in  Fig.  [6]  that  is  8  as  opposed  to  208  parallel 
fiber  channels  connecting  successive  repeater  nodes. 

Difference  in  the  bits-per-mode  rate — Further,  the 
bits  per  mode  achieved  by  the  architecture  in  [12] 
would  be  given  by  PcnPmeas/N ,  whereas  the  bits  per 
mode  achieved  by  our  modified  architecture  would  be 
PcnPmeas/2m.  The  Pmeas  of  the  former  is  higher  (due 
to  lower  loss  incurred  by  the  photons  of  the  logical  inner 
qubits  of  the  clusters  as  they  do  not  need  to  wait  in  a 
lossy  fiber  spool  while  waiting  for  the  classical  informa¬ 
tion  to  fly  back  from  the  minor  nodes).  However,  the 
other  improvements  described  below  more  than  compen¬ 
sate  for  the  better  Pme as,  and  the  latter  scheme  achieves 
a  far  better  bits-per-mode  performance  (see  Fig.  [7|. 

Linear  optic  vs.  boosted  linear  optic  fusion  gates — We 
propose  the  use  of  the  improved  Bell-state  measurement 
scheme  of  Ewert  et  al.  m  that  inject  four  single  pho¬ 
tons  to  boost  the  success  probability  of  the  fusion-II  gate. 
Our  calculations  show  that  the  cost  of  using  these  addi¬ 
tional  ancilla  photons  is  far  outweighed  by  the  effect  of 
the  improved  success  probability,  in  the  performance  of 
the  repeater  architecture,  despite  assuming  lossy  sources 
and  detectors. 

Improved  multiplexing  scheme  for  cluster  generation — 
We  use  an  improved  multiplexing  scheme  to  create  the 


clusters  at  the  major  nodes,  as  described  in  Section  III  A 
and  depicted  in  Fig.  0b).  Previous  studies  have  esti¬ 
mated  the  resource  requirements  for  cluster  generation 
based  on  the  average  number  of  attempts  required  for 
each  probabilistic  steps  pa  EES].  However,  in  order  to 
generate  the  required  cluster  at  every  repeater  station 
on  every  clock  cycle  with  high  probability,  the  resources 
required  at  each  repeater  station  need  to  be  greater  than 
the  number  that  would  allow  for  cluster  creation  “on  av¬ 
erage”.  To  our  knowledge,  this  is  the  first  study  that 
explicitly  looks  at  how  probabilistic  operations  need  to 
be  multiplexed  in  a  real  system. 

Pushing  the  measurements  ahead  during  cluster 
creation — The  single  qubits  measurements  that  do  not 
depend  on  the  outcomes  of  Bell  measurements  at  the  mi¬ 


nor  nodes,  are  performed  before  the  fusion  operations, 
directly  on  the  photons  of  the  GHZ  states,  very  early 
during  the  cluster  creation  process. 

Let  us  now  see  what  the  above  modifications  to  the 
architecture  does  to  the  rate  performance.  The  bits- 
per-mode  rates  for  the  naive  and  the  improved  schemes 
are  plotted  in  Fig.  [7ja)  and  (b),  respectively.  We  as¬ 
sume  device  loss  parameters  as  listed  in  Table  [T]  for  both 
sets  of  plots.  In  each  plot,  we  compute  the  rate-distance 
performance  (envelopes  taken  over  n,  the  number  of  re¬ 
peater  nodes)  for  four  different  error-protection  levels 
(i.e.,  k  =  7,8,9,  and  10).  For  every  point  on  each  rate- 
distance  envelope,  m  and  b  are  optimally  chosen  (consis¬ 
tent  with  the  given  k).  Each  rate-distance  plot  exhibits 
the  Drjs  =  De~saL  behavior,  and  the  exponent  s  dimin¬ 
ishes  as  a  higher  k  is  chosen.  For  the  naive  scheme,  the 
minimum  k  for  which  the  repeater  can  beat  indirect  (T) 
(pink-dashed  line)  is  k  =  8  and  the  optimized  clusters 
at  the  major  nodes  have  192  photons  each.  Hence,  the 
scheme  would  require  208  parallel  fiber  links  connecting 
successive  nodes.  In  comparison,  in  the  improved  scheme, 
k  =  7  is  sufficient  to  beat  indirect  (A),  and  requires  2m  =  8 
parallel  fiber  links.  The  optimal  tree  depth,  for  this  k  =  7 
rate  plot  is  found  to  be  d  =  2,  which  is  consistent  with 


the  analytical  development  in  Section  IV 


Table  |TT|  lists,  at  a  range  of  L  —  300  km,  and  for  each 
of  the  cases  [k  —  7,  8,  9, 10),  the  optimal  values  of  m  for 
the  naive  (mnaive)  and  new  schemes  (mnew),  the  optimal 
branching  vector  for  the  naive  (6na ive)  and  new  schemes 
( bnew )  •>  and  the  number  of  parallel  fiber  links  needed  in 
the  naive  scheme  (TVnaive).  In  the  case  of  the  new  scheme, 
the  number  of  parallel  fiber  links  needed  is  simply  2mnew. 


k 

m  naive 

X naive 

^naive 

ni  new 

bnew 

7 

5 

100 

{3,2} 

4 

{4,2} 

8 

8 

208 

{4,2} 

5 

{5,3} 

9 

11 

462 

{5,3} 

6 

{7,4} 

10 

12 

864 

{7,4} 

8 

{10,5} 

TABLE  II.  For  k  =  7,8,9,  and  10,  at  L  =  300  km  range, 
UTnaive  and  77inew  are  the  optimal  values  of  m  for  the  naive 
and  new  schemes  respectively.  6naive  and  bne w  are  the  optimal 
values  of  b  for  the  naive  and  new  schemes  respectively.  Ana ive 
is  the  corresponding  number  of  parallel  fiber  links  needed  be¬ 
tween  successive  repeater  nodes  in  the  naive  scheme.  For  the 
new  scheme,  the  number  of  parallel  links  is  2mnew 


Let  us  now  compare  the  resources  (number  of  pho¬ 
tons,  Ns)  required  to  build  the  major  node  clusters,  for 
the  respective  cases  that  can  (barely)  beat  indirect  (7/). 
The  naive  scheme  requires  1.9  x  1011  photon  sources  at 
each  major  node,  while  the  new  scheme  requires  3.3  x  106 
sources,  an  improvement  of  5  orders  of  magnitude  (see 
Fig-  0-  It  is  also  interesting  to  note  that  if  the  prim¬ 
itive  resources  were  3-photon  GHZ  sources  rather  than 
single  photon  sources,  15  thousand  GHZ  sources  would 
be  required,  a  relatively  smaller  number. 

Given  the  size  of  the  earth,  for  terrestrial  long  dis- 
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FIG.  7.  The  bits  per  mode  rates  R^k\L)  plotted  for  different 
values  of  k ,  the  numbers  of  fusion  steps,  for  the  (a)  naive 
scheme  and  (b)  with  the  improvements  of  this  paper.  The 
repeater- less  rate  bound  indirect  (L)  is  the  pink  dashed  line, 
^cluster  =  2fe  +  2  is  the  total  number  of  photons  in  the  cluster 
generated  at  each  repeater  in  every  clock  cycle. 


tance  communications,  it  is  useful  to  quantify  the  per¬ 
formance  of  our  (improved)  all-optical  repeater  scheme 
at  say  5000  km.  Without  quantum  repeaters,  the  best 
QKD  protocol  realized  with  ideal  devices  cannot  exceed 
a  key  rate  of  2.9  x  10-99  bits  per  mode  at  this  distance. 
Our  all-optical  repeater  scheme,  with  954-photon  clus¬ 
ters  (fc  =  10)  at  each  repeater  node  can  attain  a  key  rate 
of  8  x  10-3  bits  per  mode  using  2 m  =  18  parallel  chan¬ 
nels  and  n  =  12411  repeater  nodes,  which  translates  to  a 
144  kHz  key  generation  rate  assuming  a  1  MHz  repetition 
rate.  If  we  employed  518-photon  clusters  (k  =  9)  instead, 
the  rate  achieved  would  only  be  4  x  10-8  bits  per  mode 
using  2m  =  14  parallel  channels  and  n  =  12255  repeater 
nodes.  The  number  of  photon  sources  required  at  a  re¬ 
peater  node  to  create  the  required  clusters  (using  linear 
optics)  for  the  above  two  example  cluster-size  constraints 
are  1.2  x  108  and  3.6  x  107,  respectively. 

In  the  presence  of  losses  in  the  waveguide,  there  is  a 
maximum  sustainable  size  of  the  clusters  at  the  major 
nodes,  at  least  for  the  error  protection  methods  described 


in  this  paper.  A  larger  cluster  requires  a  greater  creation 
time  and  hence,  each  photon  in  the  cluster  sees  a  larger 
effective  loss  rate  (stemming  from  the  Pkh ip  term  in  etrav 
and  estat)-  Since  the  error  correction  scheme  has  a  maxi¬ 
mum  loss  tolerance  of  50%,  there  is  a  maximum  size  of  the 
clusters  that  can  be  created  and  thus  a  maximum  level  of 
error  protection  that  a  qubit  can  have.  So,  given  a  set  of 
device  losses,  increasing  the  error  protection  level  (viz., 
k )  cannot  indefinitely  improve  the  rate  performance. 

The  aforesaid  detrimental  effect  of  loss  with  an  increas¬ 
ing  cluster  size  has  more  serious  implications  for  cluster- 
state  linear  optical  quantum  computing  (LOQC)  in  gen¬ 
eral,  using  the  tree-based  counterfactual  error  correction 
technique  M-  This  is  because  a  polynomial  scaling  of 
the  number  of  photon  sources  (with  the  size  of  the  clus¬ 
ter)  is  required  in  the  asymptotic  limit  for  the  LOQC 
scheme  to  be  scalable.  The  failure  probability  of  ev¬ 
ery  qubit  needs  to  decrease  exponentially  with  the  size 
of  the  computation.  Hence,  the  level  of  protection  of 
each  qubit  must  increase  with  the  size  of  the  problem, 
which  implies  a  greater  cluster  creation  time  and  hence  a 
greater  loss  rate.  Since  there  is  a  50%  ceiling  on  the  tol¬ 
erable  photon  loss  with  the  tree  code,  it  is  not  possible 
to  achieve  the  required  level  of  protection  for  arbitrar¬ 
ily  large  computations,  as  discussed  above  for  the  case 
of  an  all-photon  quantum  repeater.  Developing  a  scal¬ 
able  method  for  creation  of  arbitrarily  large  clusters  in 
constant  time  would  solve  this  problem  and  will  also  al¬ 
low  for  a  polynomial  scaling  of  the  number  of  photons 
with  computation  size.  A  recent  paper  proposes  using 
counterfactual  error  correction  to  fault-tolerantly  create 
surface  code  data  qubits  m ■  However,  the  resource  re¬ 
quirements  for  this  scheme  are  extremely  high. 


VI.  CONCLUSIONS 

In  conclusion,  we  have  performed  a  rigorous  analysis 
of  the  resource  requirements,  and  the  achievable  secret 
key  rates  of  an  all-optical  repeater  scheme  that  improves 
upon  a  recent  proposal  [12],  while  taking  into  account 
all  the  losses  in  the  system.  While  the  all-optical  re¬ 
peater  proposal  of  [12]  presents  an  important  concep¬ 
tual  advancement,  we  show  that  it  may  not  be  prac¬ 
tically  feasible  given  its  astronomical  resource  require¬ 
ments,  both  in  terms  of  the  number  of  photon  sources 
and  detectors  needed  at  each  repeater  node,  as  well  as 
the  number  of  parallel  optical  fiber  channels  that  must 
connect  successive  repeater  nodes.  Our  scheme  improves 
the  practicality  immensely  in  both  of  the  aforementioned 
metrics,  as  well  as  the  actual  rate- vs. -distance  perfor¬ 
mance  achieved.  In  particular,  the  number  of  photon 
sources  required  at  each  node  is  reduced  by  5  orders  of 
magnitude,  and  the  number  of  parallel  channels  between 
repeater  nodes  required  to  beat  the  performance  of  a 
direct-transmission  QKD  scheme  is  brought  down  from 
more  than  two  hundred,  to  8.  These  results  suggest  that 
further  theoretical  improvements  on  quantum  photonic 
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fault  tolerant  schemes  may  further  improve  the  perfor¬ 
mance  of  all-optical  quantum  repeaters,  as  well  as  other 
applications  of  all-optical  quantum  processing.  One  of 
our  major  contributions  in  this  paper  was  to  rigorously 
prove  that  the  rate-loss  scaling  by  the  aforementioned 
genre  of  all-optical  quantum  repeaters  with  a  fixed  clus¬ 
ter  size  is  given  by  R  =  Dgs  bits  per  mode,  where  D  and 
s  are  constants  that  are  functions  of  various  device  loss 
parameters,  and  that  of  design  choices  made  (to  choose 
the  level  of  error  protection).  The  fact  that  it  is  pos¬ 
sible  to  achieve  a  value  as  the  exponent  s  <  1  proves 
the  fact  that  this  scheme  can  outperform  the  key  rates 
attainable  by  any  QKD  protocol  that  does  not  employ 
quantum  repeaters,  the  rate  performance  of  which  are 
upper  bounded  by  indirect  (v)  ~  1.44  77  for  77  1,  whose 

linear  rate-transmittance  decay  implies  s  =  1. 

In  future  work,  it  will  be  interesting  to  incorporate 
more  realistic  effects  into  the  resource-performance  trade¬ 
off  calculations  of  all-optical  repeaters,  in  particular 
mode-mismatch  errors  in  the  passive  interferometric  ma¬ 
nipulations  on  the  photons  held  locally  at  the  repeaters, 
and  multi-photon  errors  arising  from  imperfect  sources 
and  noisy  detectors.  Finally,  it  would  be  instructive  to 
analyze  and  compare  other  forms  of  quantum  repeater 
architectures,  especially  forward-error-corrected  one-way 
transmission  schemes  0,  realized  only  with  flying  pho¬ 
tons,  linear  optics  and  detectors,  but  no  quantum  mem¬ 
ories. 
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Appendix  A:  Re-ordering  measurements  in  the 
cluster-creation  process 

In  this  Section,  we  explain  why  the  A  measurements 
required  to  attach  trees  for  counterfactual  error  cor¬ 
rection  and  the  Y  measurement  required  to  create  the 
“clique”  from  the  “star”  cluster  can  be  applied  before 
the  fusion  operations.  This  makes  the  cluster  creation 
process  more  efficient.  The  reordering  of  the  operations 
is  depicted  in  Fig.  [8j  Thin  lines  here  represent  pho¬ 
tonic  qubits,  thick  lines  represent  feed- forward  opera¬ 
tions,  boxes  labelled  A,  Y,  Z,  and  H  represent  single 
qubit  A,  y,  Z  rotations,  and  Hadamard  gates  respec¬ 
tively,  and  boxes  labelled  Mx,  My,  and  Mz  represent 
measurement  in  the  A,  Y,  and  Z  bases,  respectively. 

First,  we  show  some  results  regarding  re-ordering  of 


(c)-fiTH^n=  -j  mz  | 

~j  H  | — |  Mz  |=  ~j  Mx  | 
~j~H~| — |  My  |=  -j  M-y  I 


FIG.  8.  Single  qubit  measurements  can  be  applied  before 
fusion  operations,  (a)  A  and  Y  basis  measurements  can  be 
moved  before  conditional  Z  operators,  (b)  Z  operators  before 
Z  basis  measurements  can  be  removed,  (c)  Hadamard  gates 
followed  by  measurement  in  the  A,  Y  or  Z  basis  is  equivalent 
to  direct  measurement  in  a  different  pauli  basis,  (d)  Single 
qubit  measurements  on  the  final  cluster  can  be  moved  before 
fusion  operations. 


single  qubit  measurements  and  rotations.  In  the  left  side 
of  Fig.  Ufa),  the  unitary  operation  U  on  qubit  c  is  con¬ 
ditioned  on  the  result  of  an  A  or  Y  basis  measurement 
on  qubit  b  (that  is  determined  beforehand).  In  addition, 
there  is  a  conditional  operation  Zl  on  the  qubit  b  which 
depends  on  a  feed-forward  signal  from  a  different  part  of 
the  circuit,  which  in  this  case  is  the  result  of  measure¬ 
ment  Ma  on  qubit  a.  The  application  of  a  Z  gate  before 
A  or  Y  measurement  simply  has  the  effect  of  flipping 
the  result  of  the  measurement.  Hence,  the  measurement 
Mx  (resp.  My)  can  be  performed  before  Ma  and  the 
feed-forward  result  of  Ma  can  simply  be  used  to  flip  the 
result  of  Mx  (resp.  My)  as  shown  on  the  right  side  of 
Fig.  J8^a).  The  system  in  Fig.  §b)  is  identical  to  the 
system  in  Fig.  [8](a)  except  for  the  fact  that  measurement 
in  the  A  (resp.  Y)  basis  is  replaced  by  measurement  in 
the  Z  basis.  Since  application  of  a  Z  rotation  does  not 
influence  the  outcome  of  the  Z  measurement,  the  Z  gate 
and  the  associated  feed-forward  can  be  removed  entirely. 
In  Fig.  [8^c),  we  depict  that  a  Hadamard  gate  followed 
by  an  A  basis  measurement  is  equivalent  to  a  Z  basis 
measurement,  a  Hadamard  gate  followed  by  a  Z  basis 
measurement  is  equivalent  to  an  A  basis  measurement, 
and  a  Hadamard  gate  followed  by  a  Y  basis  measurement 
is  equivalent  to  a  Y  basis  measurement  with  the  result 
flipped. 

We  now  use  these  results  to  show  how  measurements 
can  be  pushed  earlier  in  the  cluster  creation  process  at 
the  major  nodes.  The  left  side  of  Fig.  |8jd)  shows  the 
system  with  measurements  applied  after  the  fusion  oper¬ 
ations.  Single  photons  that  are  sent  through  GHZ  facto- 
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ries  to  create  3-photon  GHZ  states,  which  are  then  fused 
using  Bell  measurements  using  ancilla  photons.  The  sur¬ 
viving  photons  require  some  Hadamard  and  conditional 
Z  rotations  as  part  of  the  controlled-phase  and  parity- 
projection  operations  m ■  Finally,  some  of  the  surviving 
photons  require  X  and  Y  basis  measurements,  the  results 
of  which  are  fed  forward  to  photons  in  the  final  “clique” 


cluster.  As  shown  in  Fig.[8ja),  (b)  and  (c),  measurements 
in  the  Pauli  basis  can  be  pushed  in  front  of  Hadamard 
and  conditional  Z  rotations  by  simply  moving  to  a  differ¬ 
ent  Pauli  basis  or  flipping  the  result  of  the  measurement 
result.  Hence,  the  system  is  equivalent  to  the  right  side 
ofFig.gd)  in  which  single  qubit  Pauli  measurements  are 
applied  before  the  fusion  operation. 
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Fundamental  limits  and  improved  percolation  thresholds  for  linear  optical  quantum 
computing  using  feedback-free  probabilistic  fusion  of  photonic  micro  clusters 
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Entangled  &;-qubit  micro  clusters  can  be  stitched  together  using  probabilistic  m-qubit  linear-optic 
fusion  operations  that  succeed  with  probability  A — without  any  measurement-induced  feedforward — 
into  a  giant  connected  component  (GCC)  of  a  random  percolated  instance  of  a  target  lattice  G ,  if 
A  >  \ik,rn\  The  threshold  X<<ck,rn'>  depends  upon  the  choice  of  G  and  the  sequence  of  fusion  steps. 

This  GCC  can  then  be  renormalized  into  a  square-grid  cluster  universal  for  measurement  based 
quantum  computing.  With  two-qubit  fusion  operations  ( m  —  2),  we  give  an  explicit  construction 
that  shows  Ac3,2^  <  0.54.  This  improves  over  a  recently  published  threshold  of  0.625  when  3-qubit 
(GHZ)  states  are  used  as  the  initial  resource;  a  smaller  threshold  on  A  implies  a  higher  resilience 
to  device  losses.  We  also  show  that  A^’2^  >  l/(k  —  1)  for  all  k  >  2.  With  3-qubit  (GHZ)  states  as 
the  initial  resource,  Ac3,2'*  >  0.5,  which  implies  that  a  fusion  gate  based  on  passive  linear  optics  will 
not  work  (due  to  its  0.5  success-probability  limit)  and  that  one  must  use  fusion  gates  boosted  with 
injected  single  photons  or  2-qubit  entangled  Bell  states.  Further,  with  2-qubit  (Bell)  states  as  the 
initial  resource,  ballistic  linear-optical  quantum  computing  is  practically  impossible  using  pairwise 
fusion  operations,  since  Ac2,2'*  =  1.  We  quantify  the  device  loss  thresholds  for  ballistic  LOQC  with 
/c-qubit  micro  clusters  as  the  initial  resource,  and  pairwise  fusion  operations.  Our  results  stress  the 
importance  of  investigating  efficient  ways  to  directly  generate  entangled  clusters  of  3  or  more  qubits; 
as  well  as  the  need  to  investigate  linear  optical  methods  for  simultaneous  fusion  of  3  or  more  qubits. 

PACS  numbers:  42.50.Ex,  03.67.Dd,  03.67.Lx,  42.50.Dv 


The  Knill-Laflamme-Milburn  (KLM)  model  of  quan¬ 
tum  computing,  popularly  known  as  Linear- optical  quan¬ 
tum  computing  (LOQC),  uses  a  single  photon  in  one  of 
two  orthogonal  (spatial,  temporal  or  polarization)  modes, 
i.e.,  1 10)  and  |01)  to  encode  a  qubit,  and  uses  passive  lin¬ 
ear  optical  interferometers  and  single-photon  detectors  to 
implement  gates  and  measurements.  Gates  and  measure¬ 
ments  in  LOQC  are  inherently  probabilistic  even  if  all  the 
single-photon  sources  are  ideal,  and  all  the  linear  optics 
and  detectors  are  lossless.  Component  losses  further  re¬ 
duce  the  success  probabilities  of  gates  and  measurements, 
which  translate  into  astronomical  requirements  on  de¬ 
vices  for  problems  of  practically-relevant  size.  Since  the 
original  KLM  proposal — which  was  largely  deemed  un¬ 
scalable  due  to  the  aforesaid  reason,  several  improved 
proposals  for  LOQC  have  been  proposed  that  use  a  com¬ 
bination  of  injecting  separately-prepared  ancilla  photons, 
and  using  photon  number  resolving  detectors. 

Entangled  fc-qubit  micro  clusters  can  be  stitched  to¬ 
gether  using  probabilistic  ra-qubit  linear-optic  fusion 
operations  that  succeed  with  probability  A — without 
any  measurement-induced  feedforward — into  a  giant  con¬ 
nected  component  (GCC)  of  a  random  percolated  in¬ 
stance  of  a  target  lattice  G,  if  A  >  X^,m^ .  The  threshold 
A (k , ttl )  ^epen(^s  Up0n  the  choice  of  G  and  the  sequence  of 
fusion  steps.  This  GCC  can  then  be  renormalized  into 
a  square-grid  cluster  universal  for  measurement  based 
quantum  computing.  With  two-qubit  fusion  operations 


(m  =  2),  we  give  an  explicit  construction  that  shows 
Ac  ;  <  0.54.  This  improves  over  a  recently  published 
threshold  of  0.625  when  3-qubit  (GHZ)  states  are  used 
as  the  initial  resource;  a  smaller  threshold  on  A  implies 
a  higher  resilience  to  device  losses.  We  also  show  that 
x(k,2)  >  \/(k-  1)  for  all  k  >  2.  With  3-qubit  (GHZ) 
states  as  the  initial  resource,  Ac3’2  >  0.5,  which  im- 
plies  that  a  fusion  gate  based  on  passive  linear  optics 
will  not  work  (due  to  its  0.5  success-probability  limit) 
and  that  one  must  use  fusion  gates  boosted  with  injected 
single  photons  or  2-qubit  entangled  Bell  states.  Further, 
with  2-qubit  (Bell)  states  as  the  initial  resource,  ballistic 
linear-optical  quantum  computing  is  practically  impossi- 
ble  using  pairwise  fusion  operations,  since  Ac  ’  ;  =  1.  We 
quantify  the  device  loss  thresholds  for  ballistic  LOQC 
with  /c-qubit  micro  clusters  as  the  initial  resource,  and 
pairwise  fusion  operations.  Our  results  stress  the  impor¬ 
tance  of  investigating  efficient  ways  to  directly  generate 
entangled  clusters  of  3  or  more  qubits;  as  well  as  the  need 
to  investigate  linear  optical  methods  for  simultaneous  fu¬ 
sion  of  3  or  more  qubits. 
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•  Theory  summary 

-  Ballistic  LOQC 

•  trade-off  between  size  of  input  clusters  and  loss  tolerance 

•  Minimum  resource  state:  GHZ  states 

-  All-optical  quantum  repeater 

•  Second-generation 

-  Tree  code 

-  GHZ  sources 

•  Third-generation  scheme 

-  Parity  code 

-  GHZ  source 

-  In-line  NL  and  Ising 

•  Quantum  device  engineering 

-  How  to  generate  GHZ  resource  states? 

•  Heralded  6-wave  mixing 

•  On-demand  generation  with  QND  measurements 

-  Quantum  programmable  processor 

•  Green  machine  implementation 

•  Prospect  of  realizing  “boosted”  gates:  single-photon  injection  /  mode-matching 

-  Optical  quantum  computing  with  in-line  non-linearities 

•  Spectral  encoding  of  dual  rail  qubits  in  coupled  cavities 

•  Ising  model  implementation 

-  Comparative  device  study 

•  SNSPD 

-  Scalable  readout 

-  Thermal  compatibility  with  QPP 

•  LOQC  program  vision 


LOQC:  some  building  blocks 
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single  (dual  rail)  qubit 


beam-splitter 


Vvei6\l0)  +  y/i 


single  photon 


P(succ)  =  A 


P(fail)  =  1  -  A 


2-qubit  entangled  3-qubit  entangled 

Bell  state  G>\-\Z  state 


A  =  0.5  (max  with  linear  optics) 

A  =  0.75  (linear  optics  with  two  injected  single  photons) 
A  =  1-1 /2N  (2N-2  injected  2-qubit  Bell  pairs) 


(detector  and/or  sources  losses  reduce  all  the  above 
numbers) 


Any  graph  state  can  be  put  together  with  GHZ  states 
and  2-qubit  Type-ll  Fusion 


LOQC  (KLM  and  beyond):  single  photons,  linear  optics, 
feedforward,  and  single  photon  detection 


Ballistic  QC  with  microclusters 
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k 


Goal:  create  a  resource  that  is  sufficient  for  universal  QC  in 
a  loss  tolerant  way  using  minimum  possible  feedforward 

•  k  =  1  (Boson  sampling),  output  cluster  is  disconnected:  not  a  good  enough  resource  for  universal  QC 

•  k  =  2  (Bell  pair  initial  resource),  output  barely  percolates  if  all  efficiencies  are  100%  (no  loss  tolerance) 

•  k  =  3  (3-photon  GHZ  is  initial  resource),  then  ballistic  QC  is  possible 

•  3D-diamond  (Rudolph  et  al.Y,  A  =  62.5%;  loss  tolerance:  n  n .  >  96.2% 

U  o  (J 

•  (10,3)-b  lattice  (our  method);  \  =  54%;  loss  tolerance:  n  n .  >  93.4% 

U  o  U 

•  Best  possible  lattice  (converse  proof  /  Bethe  lattice);  Ac  =  50%;  max  possible  loss  tolerance:  91.6% 

•  Sources  of  >=  3-qubit  microclusters  needed  if  we  want  to  avoid  feed-forward,  switching,  and  associated  losses 


Bits  per  mode  S 


All  optical  repeaters 
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size  of  #  of  # 

state  single-p  on-dema 

hoton-so  nd  GHZ 

urces  states 


(111,11)  size  of  #  of  # 


state 

single-p 

on-dema 

hoton-so 

nd  GHZ 
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states 

(8,3) 

48 

200k 

lk 

(9,3) 

54 

700k 

3.5k 

(12,4) 

96 

2M 

10  k 

(18,5) 

180 

4.4M 

22k 

•  3rd  gen  +  QPC  >  2nd  gen  +  tree  code  by  lOx  in  #sources/node.  (105  x  from  Lo  et  al.) 

•  Availability  of  3-photon  GHZ  sources  reduces  the  #sources  by  >  2  orders  of  magnitude 

•  Both  can  beat  repeaterless  ideal  QKD  protocol,  even  with  losses  (in  all  devices,  coupling,  detectors). 
Mode  mismatch,  g(2)  and  excess  noise  not  included;  initial  network  results 

•  Systematic  study  of  QEC  for  optics-errors:  loss,  indistinguishability,  noise,  g(2),  spectral  purity,  etc. 
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Mlii 


In-line  non-linearities  -  Optical  Ising  solver 
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Krovi,  Guha,  Pant,  Englund  (unpublished,  2016) 
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SEQUOIA  theory  work:  summary 
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•  General  facts  we  have  learnt  about  LOQC 

-  Tradeoff:  offline  (photons/ent.  states)  vs.  inline  (NL)  resource.  Inline  often  hard 

-  LOQC  with  offline  photons/entanglement  (heralded  /  probabilistic  is  OK): 

•  Offline  states  >=  3-qubit  microclusters,  feedforward  need  drastically  reduces: 
“Ballistic”  LOQC  /  relationship  to  BosonSamp 

•  Larger  microclusters,  higher  loss  tolerance;  fewer  sources  (of  such  microclusters) 

-  Hybrid  methods: 

•  Offline  and  inline  can  be  used  simultaneously  to  trade  some  benefits  of  both 

•  CV  cluster  states  can  be  used  for  universal  QC  with  PNR  and  homodyne 

-  LOQC  with  inline:  Special  purpose  processor  /  quantum  annealing,  Ising  model 

•  All-optical  repeaters:  special-purpose  LOQC  (offline  resources) 

-  3rd  gen  +  QPC  >  2nd  gen  +  tree  code  by  lOx  in  #sources/node.  (105x  from  Lo  et  al.) 

-  Availability  of  3-photon  GHZ  sources  reduces  the  #sources  by  >  2  OoM 

-  Both  can  beat  repeaterless  ideal  QKD  protocol,  even  with  losses  (in  all  devices, 
coupling,  detectors).  Mode  mismatch,  g(2)  and  excess  noise  not  included;  initial  network 
results 

•  Challenges  stemming  by  theory  work  in  SEQUOIA 

-  Direct  generation  of  microclusters  (>=  3  qubits)  with  high  fidelity  very  important 

-  In-line  (e.g.,  cross-Kerr)  non-linearities,  even  if  small,  can  be  very  useful 

-  Svstematic  studv  of  QEC  for  ODtics-errors:  loss,  indist.  noise.  a(2L  SDectral  Duritv.  etc. 


Mode  Probability  |  Pi 
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Programmable  linear  optics  transformations 


•  Single-qubit  gate  infidelities  <  10"7 

•  PIC  demonstrated  with  up  to  26  individually  connected  output  modes 
Ill-Nitride  platform  [Soltani,  Soref,  Palacios,  Englund] 


0.20  | - 

8x8  Green  Experiment 

Machine  mode  p  —  99.9949  ±  0.0033%  ™  Simulation 

0.15  transformations 


Mode  1  Mode  2  Mode  3  Mode  4  Mode  5  Mode  6  Mode  7  Mode  8 


+  +  +  +  A  + 

++++A- 

++--++ 


+  +  - 

+  -  + 
+  -  + 

+  -  - 
+  -  - 


-  +  - 
“  -  + 
+  +  - 
+  H  + 


+ 

+ 


+ 


+ 

0 

(1 

0 

,0. 

o 

0 

0 

0 

+ 

0 

0  ! 

foj 

0 

0 

0 

0 

0 

+ 

0 1 

0 

0 

0 

0 

0 

0 

0 

+ 

0 

0 

0 

0 

0 

0 

0 

0 

+ 

0 

0 

0 

0 

0 

0 

0 

0 

+ 

0 

0 

fl 

0 

0 

0 

|0 

0 

+ 

0 

0 

0 

0 

0 

V 

0 

0 

+ 

Raytheon  im- 

BBN  Technologies  III 

Programmable  Nanophotonic  Processor 


•  Reck  and  Zeilinger,  PRL  73  (1994):  any  linear  optics  unitary  can  be  produced 
using  4-port  beamsplitters 


T  ~  lQflS 


sin(0)  e*®  cos{0) 
cos(0)  —  sin(0) 
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Perfect  optics  from  imperfect  components 
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dual  rail  encoding 
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c* 

Cl 
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Assume 
r=50%+-  2.1  % 


CPHASE 


Fidelity 


Fidelity 


J.  Mower*,  N.  C.  Harris*,  G.  R.  Steinbrecher*,  Y.  Lahini,  D.  Englund,  to  appear  in  Phys  Rev  A.  *equal  contributors 


Normalized  T 


Excellent  single-qubit  gate  fidelity 
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Voltage 


Extracted  <f>{v2) 


Voltage 


3rd  Generation  PNP 
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88  Interferometers,  26  individually  connected 
output  modes 


Mode  0 


Mode  25 


MOd*0 


Mode  25 


Column  0 


Column  10 
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PNP  Application:  8x8  Green  machine 


Operation  for  all  8  Hadamard  code  words 
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Light  sources 


Photon  Purity  gl2,(0) 

fndistingulshability 

Efficiency  t| 

Repetition  rate 

Quantum  key 
distribution 

<0,1 

Not  critical, ,  but 
consecutive 
photons  must  be 
uncorrelated 

>0.5* 

>GHz 

Cluster  state 

quantum 

computing 

<0.01  (more  study 
needed  on  many- 
photon  errors) 

>0.99 

>0.9  for  reasonable 

resources 

Ideally  GHz  to 

avoid  long 

buffers, 

maximize 

experiment 

frequency 

All-optical  quantum 
repeater 

<0.001 

>  0.99 

>0.9 

>GHz 

Bell  state  sources 
for  memory-based 
repeaters 

<0,01 

>0.9 

>0.9 

Ideally  GHz 

•  On-Demand  Single  Photon  emission  based  on  quantum  feedback  control  of  a  microcavity  (manuscript  to  be 
submitted  Sept  2016)  [Heuck,  Pant,  Englund] 

•  3-photon  GHZ  source  based  on  6-wave  mixing  (manuscript  in  preparation)  [Pant,  Soltani,  Englund,  Guha] 

•  On-demand  GHZ  source  based  on  QND  and  dynamic  cavity  control  (manuscript  in  preparation)  [Pant,  Guha, 
Englund] 
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Light  sources 


Maximum  count- 
rate  (continuous 
wave  pumping) 
[source  brightness 
[photons/pulse)] 

Lifetime 

Homogeneo 

us 

Linewidth 

at  4K 

Indistinguisha 
ble  photons 
(IP)  & 

entanglement 

(E) 

Spatial 
ta  rgeted 
fabrication 
of  single 
emitters 

Opera 

tion 

Temp. 

Integration  of 
SPEs  with 
dielectric 
cavities  or 
plasmonic 
resonators 

Color 

centers  in 

Diamond 

SiV:  “  3xlOfe  c/s  5"is 

NV:  ~  3x10s  c/s 
[0.2] 

For  other  sources 

i  27 

see  ref 

~  1  ns 

12-22  ns 

NV,  SiV  “ 

lifetime 

limited28'29 

Cr  related  - 

4  GH230. 

NV  -  (IP,  E) 

SiV  —  (IP) 

Only  for  NV 
and  SiV27 

RT 

Dielectric  -  NV, 
SiV  only 
Plasmonics  - 
NV  only 

Defects  in 
SjC,  ZnQ, 
and  BN, 

Rare 

earths  in 
YAG/YSO. 

YAG:  ~  60x10s  c/s31 
ZnO:~l*10sc/s32 
SiC:  “  2x10s  c/s” 
BN~3x106c/s3^ 

19  ns103* 

1-4  ns” 

1-4  ns33 
"3  ns3* 

N/A 

No 

No 

RT 

No 

Arsenide 

QDs 

~  lx  10 '  c/s  ^  [0.7] 

- 1  ns67* 

lifetime 

limited 

Yes 

Yes 

4K 

Yes 

Nitride 

QDs 

Not  available 

~  0.3  ns37 

~  1.5  meV37 

No 

Yes 

RT 

Dielectric -yes 
Plasmonics  -  no 

CNTs 

~3xl03c/s38 

<M  0.4  nsdB 

N/A 

No 

No 

RT 

Dielectric -yes 
Plasmonics  -  no 

2DTMDCS 

~  3.7x10s  c/s3S 

~  1-3  ns 

N/A 

No 

No 

4K 

No 
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A  source  for  scalable  photonic  QIP? 


Can  we  develop  a  source  satisfying  the  requirements  of  LOQC? 
We  investigated  three  promising  architectures: 


eff>  98%  with  spectral  multiplexing,  assuming  good 
enough  detectors  (>99.5%) 


Pump  laser  at 
frequency  2uj 


I 


Generated  colors 
at  frequency  -  uj 


Feedforward  i 


a  b  c 


•  On-Demand  Single  Photon  emission  based  on  quantum  feedback  control  of  a  microcavity  (manuscript  to  be 
submitted  Sept  2016)  [Heuck,  Pant,  Englund] 

•  3-photon  GHZ  source  based  on  6-wave  mixing  (manuscript  in  preparation)  [Pant,  Soltani,  Englund,  Guha] 

•  On-demand  GHZ  source  based  on  QND  and  dynamic  cavity  control  (manuscript  in  preparation)  [Pant,  Guha, 
Englund] 


How  to  build  a  sufficient  source? 
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r. 


Feedback 


PNRD  "Gate  D 


^iD  1  K>i 


Output 


Gate  S 


Protocol 

Generation 

•  Pump  laser  is  turned  off  at  the  first  idler  detection. 

•  A  total  of  «  pairs  are  created  and  L  idler  photons  are  detected. 
Release 

•  Coupling,  7c  is  turned  off  at  the  (Z-7)th  signal  detection. 

•  A  total  of  N  signal  photons  are  detected. 
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Numerical  Example: 


Ql 

QiD 

On-off  time 
Latency 


2T0  7 
2T04 
60  ps 
30  ps 


Detector  efficiency  :  99% 
Release  fidelity  :  99% 
Generation  fidelity  :  98% 

Success  Probability:  72% 


Currently  on-going:  use 
frequency  multiplexing  to 
reach  >  99%  overall 
efficiency 
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ear-^  review  - 

Frequency  Multiplexing 


OQM 

optimal  Quantum 


Measurements 


Time  bins 

•  •••• 

Time  bins 

i  i  i  i  Frequency  Bins 

i  i  i 

I  I  I  I 

time  time 


single  frequency  operation  frequency  multiplexing 

•  Modified  protocol 

•  In  each  time  bin,  pairs  at  multiple  frequencies  are  generated 

•  There  are  many  low-Q  idler  outputs  corresponding  to  different  frequencies 

•  Generation  of  one  photon  at  any  frequency  in  any  of  the  time  bins  heralds  a  success 

•  Efficient  frequency  conversion  at  the  end  of  the  clock  cycle  (Li  et.  al.  CLEO  201 5) 

•  Increases  the  number  of  attempts  at  pair  generation  by  the  number  of  frequencies 

•  Success  probability  >90%  can  be  achieved  by  multiplexing  over  5  frequencies  with  the  same 
parameters  used  in  the  single  frequency  operation 

•  Success  probability  >  99%  appears  possible  assuming  detection  efficiency  >99.5% 
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In-line  nonlinearities? 


Two-photon  interactions  with  in-line  nonlinearities 

-  Proposal  for  a  passive  two-photon  controlled-phase  gate  using  intrinsic  nonlinearities  of  two-level 
emitters:  manuscript  to  be  submitted  Sept  2016  [Lahini  et  al,  ArXiv:  1501. 04349] 

-  Quantum  Logic  with  Interacting  Bosons  in  ID:  manuscript  under  review  [Nysteen,  ..  Englund] 
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How  to  produce  large  cluster  states  on  PICs? 


Browne,  D,  E,,  &  Rudolph,  T,  {2005).  Resource-efficient  linear  optical  quantum  computation.  Physical  Review  Letters,  95(1  )f  010501 . 

Zhang,  Q.,  Bao}  X.  H.,  Lu,  C.  Y.h  Zhou,X.  Q.,  Yang,  T.,  Rudolph,  T.,  &  Pan,  J.  W.  (2003).  Demonstration  of  a  scheme  for  the  generation  of  "even t-ready*  entangled  photon  pairs 
from  a  single-photon  source.  Physical  Review  A,  77(G),  062316. 


Cluster  state  generation 
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fusion  1  — 


Prepared 
cluster  state 


feed -forward 
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Browne,  D.  E.p  &  Rudolph,  T.  {2005).  Resource-efficient  linear  optical  quantum  computation.  Physical  Review  Letters,  05(1 010501 . 
Raussendorf,  R.,  &  Briegel,  H.  J.  (2001 }.  A  one-way  quantum  computer.  Physical  Review  Letters,  86(22),  5188. 
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A  Cluster  State  Architecture  with  photon  storage 


dual  rail  qubit  -.. 

classical 
channel 


quantum 

channel 


switch 


-D-  detector 

delay  line 

reflector 

(►Heralded 
single 
photon 
source 


SPS 


Classical  Computer 

■ 

Fast  feed-forward 
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Nanophotonics  Outlook 


•  Challenges  and  future  directions 

-  Devices  for  microcluster  state  generation  :  off-line  nonlinearity  resource 

-  Strong  coupling  at  room  temperature  :  in-line  nonlinearity  resource 

-  Scaling  full-system  architectures 


lO  o  o  o 

o  o  o  o 


Bus 

waveguide 
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SEQUOIA  nanophotonics  summary 


•  Programmable  linear  optics  transformations 

-  Single-qubit  gate  infidelities  <  10"7 

-  PIC  demonstrated  with  up  to  26  individually  connected  output  modes 
”  Ill-Nitride  photonics  platform 

•  High-performance  light  sources 

-  Review  of  single  photon  sources:  manuscript  accepted  for  publication  in  Nature  Photonics 

-  On-Demand  Single  Photon  emission  based  on  quantum  feedback  control  of  a  microcavity:  manuscript 
to  be  submitted  Sept  2016 

-  3-photon  GHZ  source  based  on  6-wave  mixing 

-  On-demand  GHZ  source  based  on  QND  and  dynamic  cavity  control 

•  Two-photon  interactions  with  in-line  nonlinearities 

-  Proposal  for  a  passive  two-photon  controlled-phase  gate  using  intrinsic  non-linearities  of  two-level 
emitters:  to  be  submitted  Sept  2016 

-  Quantum  Logic  with  Interacting  Bosons  in  ID:  in  review 

•  Challenges  and  future  directions 

-  Devices  for  microcluster  state  generation  :  off-line  nonlinearity  resource 

-  Strong  coupling  at  room  temperature  :  in-line  nonlinearity  resource 

-  Programmable  PICs 

-  System-level  picture  of  cluster  state  generator  and  general-purpose  quantum  computer  using 
networked  nonlinear  qubit  sys 


SEQUOIA  detectors  summary 
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1. 

2. 


Developed  method  to  read  from  64  detectors  simultaneously  using  differential 
readout 


Proof-of-concept  demonstration  in  AIN 


lb 


lb 
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HH>^ 


TL  taper 


active 


^reg 
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delay 


Impedance  matched  readout  of  4  detector  chains 

Each  chain  consists  of  16  detectors,  all  integrated  on  AIN  waveguide 


16  detectors/chain 


Detectors  summary 
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Developed  compact  SNSPD  packaging  and  readout 


System  detection  efficiency  and  K“  Mlii 


reset  time  of  eight  SNSPDs 


Overview 
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Photonic  intearated  circuits 

Sinale  ohoton  detectors 

Photon  Sources 

Programmable  linear  optics 
transformations 

- 

On-Demand  Single  Photon 

Sources 

Demonstrated  for: 

Micro-cluster  sources: 

Single-qubit  rotations 

Green  Machine 

Bell  pair  sources 
-  GHZ  state  sources 

AND. 


LOQC  Program  Vision 
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•  Theory  and  protocols 

-  Novel  ECC  methods:  graph  codes,  finite-length  percolation,  ECC  for  mode 
mismatch  and  excess  noise 

-  Hybrid  qubit  modes:  CV  clusters  for  universal  QC  with  PNR 

-  Special-purpose  QC:  annealing,  repeaters,  quantum  nodes 

-  Novel  ways  to  incorporate  weak  non-linearities  for  universal  QC 

•  Nano-photonics 

-  High-rate  high-quality  sources  of  3-qubit  (or  larger)  microclusters 

-  Exploiting  non-linear  quantum  optics  in  novels  ways  (e.g.,  high-order  NL,  complex 
photonic  molecule  structures)  for  qubit  encoding  and  interactions 

-  On-chip  detectors  &  feed-forward 

•  Single  photon  detectors 

-  Demonstration  of  scalable  readout  using  time  tagging 

-  Electro-optic  feed-forward 

-  New  room-temperature  detector  concepts:  semiconductor  absorber  and 
Yamomoto-Haus  QND  approach 

•  Possible  program  goal 

-  5  yr:  Demo  of  3-qubit  GHZ  state  generation  with  99.99%  Fidelity  at  100  MHz 

-  5  yr:  Translating  SC  circuit  QED  to  optical:  »  higher  density,  RT,  networks 

-  5  yr:  10-qubit  cluster  produced  with  off-line  nonlinearities 

-  10  yr:  Demo  1  logical  qubit  (1  Rausendorf  lattice  block)  w/  Ballistic  QC,  F  =  99% 

-  15  yr:  One  all-optical  quantum  repeater  that  beats  repeaterless  performance 


Back  up  slides  (details) 


Raytheon 

BBN  Technologies 


Theory:  detailed  slides 
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General  purpose  LOQC 
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OVERVIEW  OF  LINEAR  OPTICS  QUANTUM  COMPUTING 

1.  KLM  2001  SCHEME:  SPS,  FF,  dynamic-reconfigurable  PNP 

-  Original  scheme:  Feedforward;  entangled  ancilla  clusters  for  efficient  operation 

-  Improved  schemes:  Boosted  fusion  gates 

2.  BALLISTIC  CLUSTER  STATE:  SPS,  FF  only  at  state  prep  and  cluster  measurement 

-  Sources  of  entangled  3-5  photon  cluster  fragments  improve  resource  requirements 

-  Small  multi-photon  entangled  measurements  help 

3.  Quantum  computing  using  stored  photons 

-  No  nonlinearities:  scalable  architecture  for  LOQC  and  cluster  state  quantum 
computing  (Ballistic  or  non-ballistic) 

-  With  parametric  nonlinearities  (Kerr,  chi-2,  atomic,  etc) 

-  Ising  model 

-  Parallels  to  superconducting  quantum  computing 

-  Why  is  optics  better  than  SCQC?  Room  temperature,  104  higher  density; 
interfacing  with  telecom  quantum  networks/qubit  distribution 

4.  DV  quantum  computing  using  CV  cluster 

Embed  circuit  model  quantum  computing  into  CV-cluster  states  by  Olivier 
Use  homodyne  and  direct  detection 


Nanophotonics 
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detailed  slides 


Design  for  SNSPD-integrated  PICs 
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60  K 

3  K 

Optical  +  Electrical 

Modulators 

Connections 

SNSPDs 

1  cm 


How  powerful  does  the  cryostat  need  to  be? 

Optical  connections:  220  nm  X  500  nm  silicon  waveguides  surrounded  by  lOum  X 
10  um  X  10  um  silica:  1.7  uW/waveauide 

Electrical  connections:  15  um  diameter  gold  :  1 .3  mW/wire 

Radiation:  18uW  assuming  5mm  X  5mm  radiative  surface 


Takeaway:  the  major  contribution  comes  from  electrical  connections. 


Design  for  SNSPD-integrated  PICs 
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•  PNPIb:  ~8.5dB  measured  coupling  loss  per  free-space-coupled  facet 

•  Losses  predominantly  due  to  aperture  clipping  in  objectives  and 
surface  irregularity  in  lenses 
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Low-temp,  high-speed  PICs  concept 


Design  for  electro-optic  version  of  PNP  using 
rings 

NEMS  approach 
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How  good  of  a  single  photon  source  needed? 


Mlii 


Photon  Purity  gl2,[0) 

Indistinguishability 

Efficiency  q 

Repetition  rate 

Quantum  key 
distribution 

<0,1 

Not  critical, ,  but 
consecutive 
photons  must  be 
uncorrelated 

>0.5* 

>GHz 

Cluster  state 

quantum 

computing 

<0.001  (more  study 
needed  on  many- 
photort  errors) 

>0.99 

n>o.9  for  reasonable 

resources 

Ideally  GHz  to 

avoid  long 

buffers, 

maximize 

experiment 

frequency 

All-optical  quantum 
repeater 

<0.001 

>0.99 

>0.9 

>GHz 

Bell  state  sources 
for  memory-based 
repeaters 

<0,01 

>0.9 

>0.9 

Ideally  GHz 

*  -to  be  competitive  against  attenuated  laser  QKD  with  decoy  state 
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Mlir 

Atom-like  single  photon  sources 


Maximum  count- 
rate  (continuous 
wave  pumping) 
[source  brightness 
(photons/pulse)] 

Lifetime 

Homogeneo 

us 

Linewidth 

at  4K 

Indistinguisha 
ble  photons 
(IP)& 

entanglement 

(E) 

Spatial 
targeted 
fabrication 
of  single 
emitters 

Opera 

tion 

Temp, 

Integration  of 

SPEs  with 

dielectric 

cavities  or 
plasmonic 
resonators 

Color 

centers  in 

Diamond 

SfV;  ~  SxlO6  c/s  525 
NV:  ~  3*10sc/s 
[0.2] 

Diamond  Ge-V:  >  1 06, 
[>0.9] 

~  1  ns 

12-22  ns 

NV,  SiV  “ 
lifetime 
limited28, 29 

Cr  related  - 

4  GHzi0. 

NV  -  (IP,  E) 

SiV -(IP) 

Only  for  NV 
and  SiV27 

RT 

Dielectric  -  NV, 

SiV  only 
Plasmonics  - 
NV  only 

Defects  in 
SiC,  ZnOf 
and  BN* 

Rare 

earths  in 
YAG/YSO. 

YAG:  ~  60xlO3c/s31 

ZnO:~  1x10s  c/s w 
SiC:  ~  2x10s  c/s” 

BN  ~  3xlOS c/s” 

19  ns®35 

1-4  ns32 

1-4  ns33 
"3  ns34 

N/A 

No 

No 

RT 

No 

Arsenide 

QDs 

-  l*10'c/s*yB  [0.7] 

~  1  ns6' 36 

lifetime 

limited 

Yes 

Yes 

4K 

Yes 

Nitride 

QDs 

Not  available 

-  0.3  ns3' 

~  1.5  meV37 

No 

Yes 

RT 

Dielectric -yes 
Plasmonics  -  no 

CNTs 

-3x10"  c/s38 

-  0.4  ns3B 

N/A 

No 

No 

RT 

Dielectric -yes 
Plasmonics  -  no 

2D  TMDCs 

~  3.7xio3  c/s3B 

- 1-3  ns 

N/A 

No 

No 

4K 

No 

Summary  of  photophysical  properties  of  solid  state  SPEs.  c/s  -  counts  per  second;  N/A  -  not  available.  The  reported  count  rates  for  each  system  can  be  potentially  optimized  by  integrating  with  cavities ,  or  improving  collection 
optics.  *  -  count  rate  at  the  objective ,  that  is  directly  comparable  to  other  systems.  A  brighter  count  rate  of  ~  lxlO7  c/s  was  reported  for  photons  arriving  at  the  first  collection  lens.  @  -  realized  by  optical  upconversion  to  a 
short-lived  excited  state.  $  -  reported  from  a  nanodiamond  on  iridium.  #  -  recorded  from  a  nanodiamond  positioned  onto  a  solid  immersion  lens.  In  both  cases  emitters  in  bulk  diamond  are  dimmer. 
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Single  photon  sources  as  required  for  LOQC? 


P(n=1)>0.9 

g(2>(0)<0.01  (more  theory  needed) 
HOM  visibility  >  0.99 


Suppose  we  use  spontaneous  pair  generation.  No 
published  scheme  reaches  this  performance. 
Mikkel  Heuck,  Mihir  Pant,  and  DE  proposed  a  new 
scheme  that  could  work  with  reasonable 
technology. 


Mikkel 


Mihir 


How  to  build  a  sufficient  source? 
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r. 


Feedback 


PNRD  "Gate  D 


^iD  1  K>i 


Output 


Gate  S 


Protocol 

Generation 

•  Pump  laser  is  turned  off  at  the  first  idler  detection. 

•  A  total  of  «  pairs  are  created  and  L  idler  photons  are  detected. 
Release 

•  Coupling,  7c  is  turned  off  at  the  (Z-7)th  signal  detection. 

•  A  total  of  N  signal  photons  are  detected. 


2 


<D 

Pi 

Qh 

a 

Ph  2 


10 

8 

6 

4 


01 


0 


Generation 


500  1000  1500 

Time  [ps] 


P(n  = 


Generation  Fidelity 

L  n,4 


Release  Fidelity 

i - 1 - 1 


^  P(n  =  l|<Sjv, 


iV 


First  idler 

detection 

L*  **'  ; 

A  (L-l)th  signal  ■ 

^detection 

y 

( 

1  ■ 

2.5 
2 

Q 

1.5 

i  Q 

1  CO 

£ 


0 


Numerical  Example: 


Ql 

QiD 

On-off  time 
Latency 


2T0  7 
2T04 
60  ps 
30  ps 


Detector  efficiency  :  99% 
Release  fidelity  :  99% 
Generation  fidelity  :  98% 

Success  Probability:  72% 


Currently  on-going:  use 
frequency  multiplexing  to 
reach  >  99%  overall 
efficiency 
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The  Device  Structure  -  Protocol 
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Clock-'*- 
Signal 


Clock- 


Signal 


Laser  j  j 
► - Switch 


43 


The  Device  Structure 


storage 


Darmawan  et  al.  Optics  Express 
(2007) 


Transmission  s'i 


Raytheon 

Gate  Structure  B,NT«hnolo9i“ 


■  Closed  state:  ip  =  n 

■  Open  state:  /  -k 


Wavelength  [nm] 
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ear-^  review  - 

Frequency  Multiplexing 


OQM 

optimal  Quantum 


Measurements 


Time  bins 

•  •••• 

Time  bins 

i  i  i  i  Frequency  Bins 

i  i  i 

I  I  I  I 

time  time 


single  frequency  operation  frequency  multiplexing 

•  Modified  protocol 

•  In  each  time  bin,  pairs  at  multiple  frequencies  are  generated 

•  There  are  many  low-Q  idler  outputs  corresponding  to  different  frequencies 

•  Generation  of  one  photon  at  any  frequency  in  any  of  the  time  bins  heralds  a  success 

•  Efficient  frequency  conversion  at  the  end  of  the  clock  cycle  (Li  et.  al.  CLEO  201 5) 

•  Increases  the  number  of  attempts  at  pair  generation  by  the  number  of  frequencies 

•  Success  probability  >90%  can  be  achieved  by  multiplexing  over  5  frequencies  with  the  same 
parameters  used  in  the  single  frequency  operation 

•  Success  probability  >  99%  appears  possible  assuming  detection  efficiency  >99.5% 
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Producing  entangled  states  on  demand? 


OQM 

optimal  Quantum 


Measurements 


•  3-photon  GHZ  source  can  reduce  resource  requirements 
for  cluster  state  quantum  computing  by  a  factor  of  200 

•  Use  an  n  photon  creation  process: 

pump  — >  po|0)®n  +Pi|l)®n  H-  _P2 12)^^ 

-  n  =  2:  SPDC/sFWM 

-  n  =  3:  3  photon  down  conversion  with  x(3) 

•  Scramble  the  heralding  signal  from  two  such  sources  to 
herald  without  destroying  entanglement 

-  Method  1:QND 


•  n  photon  GHZ  from  an  n  photon  process 
-  Method  2:  Direct  detection 

•  n-1  photon  GHZ  from  an  n  photon  process 


Pant  et.  al.  (in 
preparation) 


•  Use  feedback  from  heralding  to  shut  off  the  pump:  make 
the  source  deterministic 


Method  1 :  QND 
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•  Two  identical  sFWM  sources 

•  The  pumps  go  through  a  beam  splitter  before  homodyne  detection 

•  Weak  nonlinearity:  a0  is  0(1) 

-  a  is  the  number  of  photons  in  the  strong  pump  field 

-  0  is  the  phase  shift  due  to  a  single  photon 

•  n  photon  GHZ  from  an  n  photon  process 


Feedforward 


i  i 


signal 


Pump 


idler 


Clock 


i 


i 


1 

2 


1P1 

P2 


.  4 


L*L 


Phase 

shifter 


|YxY| 


Homodyne 


|XxX| 


4Yx 


Design  of  a  bell  pair  source 
using  sFWM 

Pant  et.  al.  (in  preparation) 


a 


b 


c 


Method  2:  direct  detection 
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•  n-1  photon  GHZ  from  an  n  photon  process 

-  Requires  high  n 

•  Example:  downconversion  from  515  nm  to  1550  band 
using  material  x(3)  for  the  generation  of  dual  rail  bell 
states 

•  Direct  generation  of  3  qubit  dual  rail  GHZ  states  requires 
a  4  photon  process 
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Nonlinear  Six-Wave  Mixing  Process 


Pump  laser  at 
frequency  2u> 


t 


Generated  colors 
at  frequency  ~  w 


Possible  cases: 


□  Conversion  of  two  photon  near  780  nm  into  4  photons  in  the 
1550  nm  band  using  material  x(5) 

□  Efficiency  of  the  process  is  less  than  SPDC  or  sFWM  BUT  The 
important  parameter  for  scalable  quantum  computing  is  heralding 
efficiency  which  can  be  high. 

□  Efficiency  is  proportional  to  Q5:  High  Q  resonators  can 
dramatically  boost  efficiency 

□  Our  initial  results  show  that  the  phase  matching  and  anomalous 
dispersion  required  can  be  achieved 
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Case  2 
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Case  3 
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Pant  et.  al.  (in 
preparation) 


Case  1  is  the  desired  case.  Two  ways  to  filter  out  the  other  cases: 

•  Use  photon  number  resolved  heralding.  Only  pick  cases  with  one  photon  in  the  heralding  mode 

•  The  process  requires  phase  matching  and  anomalous  dispersion.  The  other  processes  will  not  meet  these 
conditions  and  will  be  naturally  suppressed. 


- - 

Clp 


Casel 


m+/?  a>+a  co  co— a  co—fi 


I  mmmn 


Generation  of  3  photon  GHZ  sources  from  six 
wave  mixing  sources _ 
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Pump  laser  at 
frequency  2 w 


Single  photon  detectors 
Detector  1 
Detector  2 

Phase  shift  of  tt  if 
photon  detected 


qubit  1 


qubit  2 


/  qubit  3 


^-uj—a^-(jU—/3  lcu+a  lcu+/3  )  A  |  0<u  — q;0c<j  —  B  |0a;— oS^uj— A  |  leu— a  laj  —  (3  l<u+a  l(u+/?)  B 


V2 

/3^-uj +(3) a\Q 00-/3^00+ -\-/3) B  H“  |0cj— /^Ocj+aOcu-b/?)  A  leu— /3lcu+a;la;+/3).B 


\/2 

Raytheon  BBN  Proprietary 
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Six  wave  mixing  vs  GHZ  creation 
from  single  photons - 


Six  wave  mixing  will  results  in  probabilistic,  heralded  creation  of  GHZ  states. 
Feedforward  on  the  heralding  signal  is  required  to  make  the  source  deterministic 


Why  is  this  better  than  probabilistic  generation  of  single  photons  from  SPDC 
followed  by  postselected  linear  optics? 


SPDC/sFWM->Single 
photons  ->postselected 
linear  optics  ->  GHZ 

Six  wave  mixing  ->  GHZ 

Average  number  of 
attempts 

128 

10.3 

Number  of  Feedforward 
steps 

2 

1 

Six  wave  reduces  resource  requirements  and 
feedforward! 
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Deterministic  Gates  with  Interacting  Bosons 


(a)  \Q) 


(c)  Phase  gate 


M  ID 


(a)  Real  (U) 


(b)  Imaginary  {U) 


Site  number 


Lahini,  Y.,  Steinbrecher,  G.  R.,  Bookatz,  A.  D.,  &  Englund,  D.  (2015).  arXiv:1 501. 04349. 
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Proposal  for  2-photon  nonlinear  gate 


OO 


d k£(k)t{k)  aJc(Ar)|0) 


—oo 


o 


<p/2 


1 

QtunEutu 

IrnMcr 

Dmlioul 

Coupler 

£B 

Phii e 
Shuflcr 

1 

*— 

two-way 


Scattered  two-photon  state  fidelity 
compEircd  la  [he  ideal  stale  -I "leJ’l  ls?1 


t-'=£ttkli1yJ=phuKC  measure  ideal) 


Example:  For  ^mii'onductnr  quiuituni  del  in  PC  Waveguide , 
r/(Y+D=0.9B  =>F~9Q%. 


■  Anders  Mystfien.  MfKkel  H£U£l&»  Dara  P  S,  McCutcheon*  Jesper  MjjrJt,  and  Dirk  R,  Englund,  Proposal  for  a 
passive  two-photon  control  led-phase  gate  using  intrinsic  non-linearities  of  two-level  emitters,  to  be  submitted 
(3016) 

■  Full  set  of  quantum  logic  gates  in  coupled  resonators  under  Bose -Hub  band  Hamiltonian; 


Hjgh-Wiilily  Cranium  Logic Gnho^wilh  Inlw^cng  Sokhu-mi  #  ID  Laflioq,  ^nvl^nhiw,  Gtfjggfy  ft  Slflinhmchnr  Adnm  D.  Boafcl Cirtt E ngluid,  Jffflfoi  1501.114349 


Design  study 

Other  Photonics  Platforms 
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Diamond 

GaN 

AIN 

Si02 

Si3N4 

refractive  index 
(5>  A  =  500  nm 

2.4 

2.42 

2.13 

1.45 

2.04 

bandgap*  (nm) 

230 

365 

200 

140 

250 

Crystalline 

Yes 

Yes 

Yes 

No 

No 

electro-optic  coeff. 

NA 

~  1  pm/V 

~  1  pm/V 

NA 

NA 

thermo-optic  coeff.  (K 1) 

lO6 

1.6xl04 

3.6  x10  s 

10  s 

2.5x10  s 

thermal  conductivity  (W.m^.K 1 } 

2200 

130 

285 

1.4 

30 

active  integration 

No 

Yes 

Yes 

No 

No 

Design  study 


AIGaN  Wide-Bandgap  Photonics 
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Bandgap  (nm) 


w 


AIN 


Mohammad 
Soltani,  Richard 
Soref,  Dirk 
Englund 


(a)  Cross  section  of  the  AlxGal-xN  waveguide,  (b)  Variation  of  AlxGal-  xN/AMsl  lattice-mismatch  vs.  x.  The  top  horizontal  axis 
shows  the  bandgap  of  AlxGal  -xN  in  the  wavelength  unit  for  the  each  x  value. 


Refractive  index  difference  between  the  AlxGal -xN  and  AIN  for  (a)  ordinary  and  (b)  extraordinary 
indices  with  the  x  values  and  lattice  mismatches  shown  in  the  inset  of  (a). 


Design  study 

AIGaN  Photonics 
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ring  radius  (micron) 


10  15 

ring  radius  (micron) 


20 


(a)  Simulated  radiation  Q  for  the  TE  mode  of  an  AI0.65Ga0.35N  ring 
resonator  vs.  its  radius  for  a  resonance  wavelength  ~  300  nm  and  for  h  =  0 
and  h  =  H/2.  The  inset  shows  the  ring  cross  section.  For  this  simulation  W 
=  700  nm,  H  =  350  nm.  (b)-(c)  Cross  section  of  the  radial  electric  field  mode 
profiles  for  a  ring  radius  of  15  microns  for  h  =  H/2  and  h  =  0,  respectively, 
(d)  Roundtrip  loss  of  the  rings  in  (a). 


Monolithic  Visible-to-IR  Quantum  Photonic 
Fabrication  on  a  Foundry  Si  platform 
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Coupling  of  light  Electro-optic  (Pockels) 

between  AIN  and  SiN  , - switching - 

/ 


A  path  to  all  the  required  operations  on  the  same  integrated  photonics 
platform  (through  AIM  Photonics.  SUNY  Fab) 

•  Fast  low-loss  modulation  for  feed-forward  in  AIN  using  pockels  effect 

•  Efficient  detection  in  Si 

•  Broadband  nonlinear  optics  for  generation  of  a  variety  of  heralded  sources  in  SiN  (x<3) 
and  x(5>)  and  AIN  (x(2)) 

•  Compatibility  with  SNSPD  integration. 
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Monolithic  Visible-to-IR  Quantum  Photonic 
Fabrication  on  a  Foundry  Si  platform 
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Coupling  of  light  Electro-optic  (Pockels) 

between  AIN  and  SiN  , - switching - 

/ 


A  path  to  all  the  required  operations  on  the  same  integrated  photonics 
platform  (through  AIM  Photonics.  SUNY  Fab) 

•  Fast  low-loss  modulation  for  feed-forward  in  AIN  using  pockels  effect 

•  Efficient  detection  in  Si 

•  Broadband  nonlinear  optics  for  generation  of  a  variety  of  heralded  sources  in  SiN  (x<3) 
and  x(5>)  and  AIN  (x(2)) 

•  Compatibility  with  SNSPD  integration. 
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Getting  photons  under  control 
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Nir 


Photonic  qubits  are  out  of  control! 

Capture  them  in  resonators 

•  Interactions  by  detection 

•  Interactions  by  nonlinearity 

— ►  Closely  analogous  to  SCQC  in  the  optical 

domain! 

Intermediate  devices/spin-offs  are  already 

important! 

-  Example:  on-demand  entangled  pair  sources 
would  solve  the  arguably  biggest  hurdle  in 
quantum  networking  with  rare  earth  ions  and 
atomic  gases 
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Ising  model  solver  based  on 
Quantum  Annealing - 

i?>  i?> 


Photonic  Molecule 


qubit  1  qubit  2  qubit  3  qubit  4  (jJ, 

Frequency  encoded  dual  rail 
qubits 

Photon  creation:  sFWM  in  the  ring.  Temporal  multiplexing  in  high  Q 
ring  for  near  deterministic  creation. 


Single  qubit  beam-splitter:  Frequency  conversion  by  ring  modulation 
or  Bragg  scattering  FWM 


•  Single  Qubit  Phase:  Phase  on  a  single  qubit  can  be  imparted  by 
interaction  with  a  classical  pump  mode  with  cross-kerr  interaction. 


Hundreds  of  qubits  in  a 
ring.  Minimal  interaction 
of  light  outside  the  ring 
allows  for  higher  Q. 


Two  qubit  phase  gates:  The  modes  will  interact  due  to  cross  kerr 
nonlinearity  in  the  ring  which  will  be  enhanced  due  to  the  high  Q  of  the 
rings. 

The  single  and  two  qubit  phase  gates  can  be  tuned  up  linearly  in  time 
by  increasing  the  interaction  time  between  single  qubit  beam-splitter 
operations  i.e.  by  changing  the  time  spent  in  each  step. 


A  nonlinear  quantum  computing  Raytheon 

1  1  ^  BBN  Technologies 

architecture  inspired  by  superconducting 


9Fico/2  - 
7fi<»/2  - 
5fico/2  - 
3fuo/2  - 
fito/2  - 


Nonlinearity 


9Fico/2  - 
7tico/2  - 
5fico/2  - 
3fico/2  - 
fito/2  - 


•  High  Q  ring  resonators  can  allow  splitting  of  the  energy  levels  similar  to 
superconducting  systems 

•  Optical  qubits  offer  several  advantages  over  superconducting  qubits 
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Microwave  vs  optical  photon  QC 


Superconducting  qubits 

Optical  qubits 

Coherence  time  (T1  &  T2) 

~  100s  of  ps 

~  1 00s  of  ps 

Carrier  frequency  (Hz) 

~  109 

~  1014 

Qubit  area  (pm2) 

104 

10 

Operating  temperature  for 
~  zero  thermal  photons 

10s  of  mK 

Room  temp 

Nonlinearity:  Acosi 
„ .  Ik 

photon 

~102 

In-line:  Assume 
n2~1 0'12cm2/W  for  GaAs. 
For  one  photon  in  PhC 
cav  (V  -0.1  (A/n)3), 

An~10  ,  for  polymer 
(assume  n2~10"7cm2/W) 
and  fractal  cavity 

Strong  coupling  cavity 
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SNSPD:  detailed  slides 
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SNSPD-related  tasks 
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Mlir 


1.  Fabricate  64  functioning  single-photon 
detectors  on  PIC  with  >  90  %  efficiency 

2.  Investigate  on-chip  schematic  for  64 
detectors  parallel  read-out . 

3.  Implement  64-channel  fiber-coupled 
off-chip  detector  system. 


Aluminum  Nitride  for  integrated 
photonics 
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■  Large  bandgap  (6.23  eV  @  77K) 

■  Transparent  from  200  nm  to  10  pm 

■  Low  fluorescence 

■  Piezoelectric  &  electro-optic  effect 

■  Optical  modulator  for  feedforward 
operation 

Require  compatible  fabrication 
processes  for  SNSPDs  ! 

Xiong  et  al.,  New  Journal  of  Physics  (2012) 
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High  quality  NbN  films  on  AIN 
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200  nm  AIN  on  sapphire  substrate 
Reactive  DC  magnetron  sputtering 
@  840  DegC 


RMS  roughness:  0.456  nm 
Thickness:  ~  4.8  nm 


Sheet  Resistance  (Ohm/sq) 
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y  (nm)  y  (nm) 


Absorbance  of  waveguide  coupled 
SNSPD 
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Single  mode  ridge  waveguide  (450  nm  x  200  nm) 


x  (nm) 


TE:  n  =1.83  +  7.52x10 '3/ 

eff 

Absorption:  0.64  dB/pm 


x  (nm) 

TM:  n  =1.779  +  6.31x10 '3i 

eff 

Absorption:  0.55  dB/pm 


Tuning  detection  efficiency  to  achieve 
correlation  measurements  on  chip 


50%  absorption:  5  pm,  95%  absorption:  21.7  pm 
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SNSPD1 


SNSPD2 


N> 

N> 


10  |jm 

EHT  =  3.00  kV 

Signal  A  =  SE2 

Date 

:1  May  201 6 

1  1 

WD  =  3.3  mm 

Mag  =  1.28  KX 

Time 

:1 5:59:59 

Normalized  count  rate 


Characterizations  of  the  detectors 
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Saturated  detection  efficiency 


1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 

0.5  0.6  0.7  0.8  0.9  1 

Normalized  bias  current  (l„/lsw) 


•  2.5  K  cryogenic  probe  station 

•  Back  illumination 

•  Saturation  indicates  high 
quantum  efficiency 
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Detector  timing  performance 


T - 1 - T 


0  5  10  15 


time  (ns) 


6  ns  reset  time 


-200  -100  0  100  200 

time  delay  (ps) 

52  ps  timing  jitter 
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Scalable  integration  of  64  detectors  on  PIC 


Impedance  matched  readout  of  4  detector  chains 

Each  chain  consists  of  16  detectors,  all  integrated  on  AIN  waveguide 


Key  components  in  a  detector  chain 
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(D y™ 


4 


dela 


X 


•  • 


lb 


detecto 


Microstrip  line 


Detector 


Delay 


ine 


A 

&i02 

r - |N 

A 

■ai 

NO 

2 


Styong  mode  confinement 
Au  layer  (grounding  plane) 
shields  detector  from  stray 
light 


Signal  A  =  InLens 
Mag  =  1.84  KX 


Date  29  Feb  2016 
Time  :1 1:36:31 


80-nm-wide  2-SNAP  Provide  71  ps 
High  efficiency  and  delay,  enough  to 
SNR,  and  faster  reset  separate  two 

adjacent  detectors 


Taper 


Match  1.5  kO 
nanowire  to  50  Q 
readout  circuitry 


16  detectors/chain 


Multiple  photon  detections  in  the  time-miSSlSeSSisl  Mlii 
SNSPD  arrays 


SPICE  simulation  results 


Details  of  the  pulse  shape  could  give  a  full  information  of 

photon  numbers,  arrival  times  and  locations 


voltage  (V) 


Preliminary  results  of  a  single  20-mm 
long  nanowire 
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50,000  photon  detection  events 
(flood  illumination  over  the  entire  area) 
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16  two  photon 
detection  events 


time  (ns) 


For  64-channel  off-chip  detector  system,  we 
need  a  compact  fiber-coupling  method. 
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Optical 


Zircon  ia 
sleeve 


Packaging  used  at  NIST  and  JPL  is  as  big  as  penny. 


a 


iber 


SNSPD 

(inside 

sleeve) 


RF  pins 


A.  J.  Miller  et  a!., 


Optics  Express, 


2011 


We  designed  a  new  packaging  to  couple 
eight  fibers  to  a  row  of  eight  SNSPDs. 
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Fiber  array 

DC  rknnnArktnro 


Raytheon 

We  tested  an  eight  SNSPD  chip  using  thi!SN  Technologies 
approach _ 


30  K  stage 


RF  cables 


3  K  stage 


Optical 

fibers 


Sample 

mount 


RF 

cables 


vve  measured  me  sysiem  aeiecn®wieon  m- 
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efficiency  and  the  reset  time  of  the  eight 


SNSPDs. 


5.5 
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We  are  replacing  RF  coaxial  cables  with  coplanar 
waveguides  of  tape. 


B.  Mazin,  LTD  ‘16 


S12  (dB) 


we  expeci  more  auenuaiion  TrorTram«on  mi-: 
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coplanar  waveguides,  but  we  can  use 
nTrons 


2.5-mm-long  coplanar  waveguide 


Frequency  (GHz) 


Drain 


CftOlte 


Gate" : 


Source 


'VrOrwn 


A.  N.  McCaughan  and 
K.  K.  Berggren, 
Nanoletters,  2014 


Future  steps 
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Nir 


64  SNSPDs  on  PIC  with  >  90  %  efficiency  & 
on-chip  schematic  for  64  detectors  parallel 
read-out . 

•  Fabricate  the  64  detectors  on  AIN  waveguides  and 
microstripline  structure. 

•  Test  the  64-SNSPD  PIC. 

64-channel  fiber-coupled  off-chip  detector 
system 

•  Start  using  fibers  with  higher  NA 

•  Integrate  more  detectors  in  a  cryostat  that  can  cool  down 
to  1  K. 


Time-tagged  readout 


Photon  arrival  time  given  by  (t  +t  )/2 
Photon  arrival  position  given  by  (t  -t  ) 
v 

g 

!b 
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